Data-Driven Seismic Waveform Inversion: A Study on the Robustness and Generalization
Not Available
58 GEOSCIENCES↗
SEARCH · Search NASA
Results for “Geosciences (General)”
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.
Convection-Permitting Simulations With the E3SM Global Atmosphere Model
This paper describes the first implementation of the Δx = 3.25 km version of the Energy Exascale Earth System Model (E3SM) global atmosphere model and its behavior in a 40-day prescribed-sea-surface-temperature simulation (January 20 through February 28, 2020). This simulation was performed as part of the DYnamics of the Atmospheric general circulation Modeled On Non-hydrostatic Domains (DYAMOND) Phase 2 model intercomparison. Effective resolution is found to be ~ 6x the horizontal dynamics grid resolution despite using a coarser grid for physical parameterizations. Despite this new model being in an immature and untuned state, moving to 3.25 km grid spacing solves several long-standing problems with the E3SM model. In particular, Amazon precipitation is much more realistic, the frequency of light and heavy precipitation is improved, agreement between the simulated and observed diurnal cycle of tropical precipitation is excellent, and the vertical structure of tropical convection and coastal stratocumulus look good. In addition, the new model is able to capture the frequency and structure of important weather events (e.g., tropical cyclones, extratropical cyclones including atmospheric rivers, and cold air outbreaks). Interestingly, this model does not get rid of the erroneous southern branch of the intertropical convergence zone nor the tendency for strongest convection to occur over the Maritime Continent rather than the West Pacific, both of which are classic climate model biases. Several other problems with the simulation are identified, underscoring the fact that this model is a work in progress.
54 ENVIRONMENTAL SCIENCES↗
Development, verification, and validation of comprehensive acoustic fluid-structure interaction capabilities in an open-source computational platform
The acoustic fluid-structure interaction (FSI) formulation is a practical numerical approach for the seismic analysis of fluid-filled tanks. However, there are no verification and validation studies reported in the literature that demonstrate the ability of an acoustic FSI numerical model to predict responses important to structural and mechanical design for intense translational and rotational earthquake inputs. Herein, an acoustic FSI formulation is implemented in the open-source Multiphysics Object-Oriented Simulation Environment (MOOSE), and is formally verified and validated using analytical solutions and code-to-code verification, and experimental data, respectively. The analytical solutions are for small amplitude, unidirectional seismic inputs. The code-to-code verification utilizes a previously verified and validated Arbitrary Lagrangian-Eulerian (ALE) numerical model in the commercial finite element code LS-DYNA. The validation studies utilize a comprehensive data set assembled from results of 3D earthquake-simulator tests of a fluid-filled vessel. The acoustic numerical model in MOOSE is verified and validated for hydrodynamic pressures and support reactions except for cases that involve significant convective response. For small amplitude inputs, numerically predicted wave heights match those of the analytical solutions. The numerical model is not verified and validated for wave height calculations under intense 3D seismic inputs. The run times for the acoustic FSI simulations in MOOSE are an order of magnitude, or more, shorter than for the corresponding ALE simulations in LS-DYNA. The utility of the MOOSE acoustic FSI implementation is demonstrated by seismic analysis of a building equipped with a fluid-filled, advanced nuclear reactor.
22 GENERAL STUDIES OF NUCLEAR REACTORS↗
Characterization of oceanic post-cold frontal clouds and their model representation. Final report
Low-level clouds are ubiquitous over the oceans, and their impact on the Earth’ radiative budget is important. However, general circulation models (GCMs) that are used to represent our climate and its evolution still experience problems to represent their correct amount and radiative impact. In the extratropical latitudes (~30-60 N/S), cloud properties are strongly modulated by the occurrence of anti-cyclones (high pressure systems) and cyclones (low pressure systems). While deep clouds dominate in cyclones near the center and along the cold and warm fronts, low-level clouds populate the colder regions of the cyclones, and more specifically the region in the wake of the cold fronts, the post-cold frontal region (PCF). The GCMs mentioned above were found to underestimate the amount of these PCF clouds, and this caused errors in their representation of the amount of solar radiation at the surface of the southern hemisphere oceans, causing an excess in absorption and errors in long-term predictions. To better understand which processes are ill-represented in the GCMs to cause such issues, we used ground-based observations from the Eastern North Atlantic (ENA) ARM site and a high resolution regional model (called WRF) to examine the properties of PCF clouds. With the WRF model, we simulated the passage of a cold front at the ENA site using theoretically distinct representations of convection and the physics of the boundary layer. While these physical representations had little impact on the timing, structure and circulation of the cold front, we found a high sensitivity to how the large scale information was fed into the model. Errors rapidly developed if the incoming flow information was injected too far from the site, but 1000 km was found to be an optimal distance between the outer boundary and the ENA site for a realistic cold front passage. With this model, we then conducted a series of experiment where convection and boundary layer schemes were changed, and revealed that 1) the clouds were more sensitive to the convection than boundary layer schemes and 2) the interaction between the two has a large impact. Furthermore, changing the representation of the physics in the model had greater impacts on PCF clouds than changing the initial conditions. As large-scale models increase their spatial resolution, convection will eventually be resolved but it remains that their representation of the physics of the boundary layer will matter for the correct representation of shallow (low-level) clouds. Observations were used to find which large-scale information had a clear relationship with PCF cloud and precipitation properties, in view of helping with GCM parameterizations and evaluation. Both lower troposphere stability and surface forcing were found to be important for PCF clouds that experience more vigorous dynamics than their tropical or subtropical counterparts. For cloud base and top heights, a simple measure of the change in potential temperature between the surface and the 800 hPa was found to be a good predictor. This relationship was found to hold as well when using data from campaigns recently conducted in the southern ocean south of Tasmania. A recent GCM was found to succeed in representing similar relations. For cloud depth, water content, and different characteristics of precipitation both the temperature contrast between the surface and near-air and the near-surface wind speed had a strong relationship. To summarize, through observational constraints and regional model experiments, the project has provided guidance to evaluate GCMs and tools to improve their representation of low-level clouds in the extratropics. It has also highlighted differences in behavior between these clouds and those typically found in the subtropics, thus opening new avenues of research for shallow clouds in general.
54 ENVIRONMENTAL SCIENCES↗
An aggregation-based nonlinear multigrid solver for two-phase flow and transport in porous media
A nonlinear multigrid solver for two-phase flow and transport in a mixed fractional-flow velocity-pressure-saturation formulation is proposed. The solver, which is under the framework of the full approximation scheme (FAS), extends our previous work on nonlinear multigrid for heterogeneous diffusion problems. The coarse spaces in the multigrid hierarchy are constructed by first aggregating degrees of freedom, and then solving some local flow problems. The mixed formulation and the choice of coarse spaces allow us to assemble the coarse problems without visiting finer levels during the solving phase, which is crucial for the scalability of multigrid methods. Specifically, a natural generalization of the upwind flux can be evaluated directly on coarse levels using the precomputed coarse flux basis vectors. We find that the resulting solver is applicable to problems discretized on general unstructured grids. The performance of the proposed nonlinear multigrid solver in comparison with the standard single level Newton's method is demonstrated through challenging numerical examples. It is observed that the proposed solver is robust for highly nonlinear problems and clearly outperforms Newton's method in the case of high Courant-Friedrichs-Lewy (CFL) numbers.
58 GEOSCIENCES↗
The Large-Scale Vorticity Balance of the Antarctic Continental Margin in a Fine-Resolution Global Simulation
The depth-integrated vorticity budget of a global, eddy-permitting ocean/sea ice simulation over the Antarctic continental margin (ACM) is diagnosed to understand the physical mechanisms implicated in meridional transport. The leading-order balance is between the torques due to lateral friction, nonlinear effects, and bottom vortex stretching, although details vary regionally. Maps of the time-averaged depth-integrated vorticity budget terms and time series of the spatially averaged, depth-integrated vorticity budget terms reveal that the flow in the Amundsen, Bellingshausen, and Weddell Seas and, to a lesser extent, in the western portion of East Antarctica, is closer to an approximate topographic Sverdrup balance (TSB) compared to other segments of the ACM. Correlation and coherence analyses further support these findings, and also show that inclusion of the vorticity tendency term in the response (the planetary vorticity advection and the bottom vortex stretching term) increases the correlation with the forcing (the vertical net stress curl), and also increases the coherence between forcing and response at high frequencies across the ACM, except for the West Antarctic Peninsula. Finally, these findings suggest that the surface stress curl, imparted by the wind and the sea ice, has the potential to contribute to the meridional, approximately cross-slope, transport to a greater extent in the Amundsen, Bellingshausen, Weddell, and part of the East Antarctic continental margin than elsewhere in the ACM.
58 GEOSCIENCES↗
Planetary Defense: Hazards Modeling [Slides]
Hydrocodes were quickly adopted for planetary impact simulations and are now being combined with other codes (e.g. tsunami, atmospheric transport, general circulation, etc.) for hazard assessment and quantification.
58 GEOSCIENCES↗
Bryan Mound Collar Compression Study
This report presents a study of collar compression for the wells at the Bryan Mound site. The study includes a comprehensive review of all available Multi-Arm Caliper (MAC) logs from the site, which has been used to identify specific wells that are experiencing casing restrictions at collar joints. Additional site-wide analysis of Bryan Mound well temperature data, proximity to sulphur wells, pulsed eddy current (PEC) logs, and a general spatial analysis has been included.
02 PETROLEUM↗
Systemic racial disparities in funding rates at the National Science Foundation
Concerns about systemic racism at academic and research institutions have increased over the past decade. Here, we investigate data from the National Science Foundation (NSF), a major funder of research in the United States, and find evidence for pervasive racial disparities. In particular, white principal investigators (PIs) are consistently funded at higher rates than most non-white PIs. Funding rates for white PIs have also been increasing relative to annual overall rates with time. Moreover, disparities occur across all disciplinary directorates within the NSF and are greater for research proposals. The distributions of average external review scores also exhibit systematic offsets based on PI race. Similar patterns have been described in other research funding bodies, suggesting that racial disparities are widespread. The prevalence and persistence of these racial disparities in funding have cascading impacts that perpetuate a cumulative advantage to white PIs across all of science, technology, engineering, and mathematics.
99 GENERAL AND MISCELLANEOUS↗
The Stability Transition From Stable to Unstable Frictional Slip With Finite Pore Pressure
Pore fluids are ubiquitous throughout the lithosphere and are commonly invoked as the cause of induced seismicity and slow earthquakes. We perform lab experiments to address these questions for drained fault conditions and low pore pressure. We shear simulated faults at effective normal stress (σ' n ) of 20 MPa and pore pressures Pp from 1 to 4 MPa. We document the full range of lab earthquake behaviors from slow slip to elasto-dynamic rupture and show that slow slip can be explained by the slip rate dependence of the critical rheologic stiffness without dilatancy hardening or other fluid effects. Our fault permeabilities ranges from 10 –18 to 10 –17 m 2 with an initial porosity of 0.1 and estimated fluid diffusion time ≈1 s. Slow slip and quasi-dynamic fault motion may arise from high P p at higher pressures but dilatancy strengthening is not a general requirement.
58 GEOSCIENCES↗
In-Situ Bioleaching of Manganese by Dissimilatory Reduction
This report describes the development of a biological leaching process for recovery of manganese from low-grade ores, with a high degree of selectivity against contaminants such as iron. This results in the production of manganese that is suitable for battery manufacture and other electrical applications. The leaching process makes use of a community of metal reducing organisms that solubilize manganese at a pH of approximately 4.5. These organisms are nourished by simple organic molecules such as acetate that are generated by decomposition of biomass. A series of long-term laboratory experiments were carried out to determine the necessary operating parameters, followed by construction of a small pilot plant processing approximately 100 kg of ore. This pilot unit was operated for two years, demonstrating the ability to consistently produce high-grade manganese at a commercially viable rate.
25 ENERGY STORAGE↗
Optimal Experimental Design for Inverse Problems in the Presence of Observation Correlations
Optimal experimental design (OED) is the general formalism of sensor placement and decisions on the data collection strategy for engineered or natural experiments. This problem is prevalent in many critical fields such as battery design, numerical weather prediction, geosciences, environmental and urban studies. State-of-the-art computational methods for experimental design do not accommodate correlation structure in observational errors produced by many expensive-to-operate devices such as X-ray machines, radars, and satellites. Discarding evident data correlations leads to biased results, higher expenses, and waste of valuable resources. We present a general formulation of the OED formalism for model-constrained large-scale Bayesian linear inverse problems, where measurement errors are generally correlated. The proposed approach utilizes the Hadamard product of matrices to formulate the weighted-likelihood, and is valid for both finite as well as infinite-dimensional Bayesian inverse problems. Extensive numerical experiments are carried out for empirical verification of the proposed approach using an advection-diffusion model, where the objective is to optimally place a small set of sensors, under a limited budget, to predict the concentration of a contaminant in a closed and bounded domain.
Correlated observations↗
Toward more mechanistic representations of biogeochemical processes in river networks: Implementation and demonstration of a multiscale model
A recently introduced multiscale model for representing the effects of small-scale hyporheic-zone biogeochemical processes is extended from the reach scale to river network scales. The model uses advection-dispersion-reaction equations for the channel network and one-dimensional advection-reaction subgrid models for the hyporheic zone. We summarize implementation in the integrated surface/subsurface hydrology modeling system Amanzi-ATS. The extension uses topologically defined meshes to represent stream/river networks and associates a hyporheic subgrid model with each channel grid cell. Biogeochemical reaction modeling capability residing in community software is accessed through an application programming interface. The implementation is verified against independent numerical solutions on a single reach. Mesh convergence studies show that commonly used semi-distributed representations can introduce significant spatial discretization error. Denitrification of farm runoff in a subbasin of the Portage River Basin in Ohio, USA is used to demonstrate the general-purpose reactive transport capability.
54 ENVIRONMENTAL SCIENCES↗
Data Report: Pore Water Geochemistry at Green Canyon 955, deepwater Gulf of Mexico
Two holes were drilled in Green Canyon Block 955 (GC 955) in the deepwater Gulf of Mexico: Hole GC 955 H002 (H002) and Hole GC 955 H005 (H005). Here we present the solute concentrations and oxygen and hydrogen isotope ratios of pore water extracted from sediment whole-round samples collected from these holes. The analytical methods, degree of sample contamination, and the method for contamination correction are discussed as are the general trends in the geochemical data.
03 NATURAL GAS↗
Ice Dusty Plasma (Final Report)
The NSF/DOE funding described in this final report was for constructing a substantial upgrade that would enable operation at variable and much lower cryogenic temperatures than was available in a previous experiment that operated from 2013 to 2019. Components for the upgraded experiment were designed and ordered from August to December 2019 and assembly began in early 2020. However, COVID-19 halted all lab work for much of 2020 and then when lab work resumed, numerous restrictions on lab access and in-lab time further impeded progress. Because of these delays, first plasma was obtained in late 2021 and a normal work pace resumed only in mid-2022. A no-cost extension was requested because the COVID-19 delays prevented work from being done at the anticipated rate. Constructing a water-ice experiment was a significant gamble because the parameter regime differs substantially from conventional laboratory plasmas and injecting water vapor into a plasma goes very much against convention as one normally goes to great lengths to avoid water vapor. This experiment opens a new frontier, namely the study of how charged ice grains form, grow, and interact with a plasma environment. The just-completed upgrade was a further gamble because it combines technology never previously used together (cryocoolers, plasma, and ice).
36 MATERIALS SCIENCE↗
ORION Code Infrastructure Upgrades
Prior to this update, ORION was capable of processing cases and data that spanned a single UTM zone and assumed that all spatial inputs were given in units of meters. During the 2024 SMART bridge funding period, we completed upgrades to the code that permit general spatial projections, basin-scale datasets, and units.
58 GEOSCIENCES↗
Synoptic and Mesoscale Modulation of Dynamic and Thermodynamic Impacts on Central Arctic Sea Ice During MOSAiC
The primary objectives of this proposed project are to determine 1) to what extent dynamic and thermodynamic air‐ice interactions in the Central Arctic are driven by atmospheric synoptic and mesoscale phenomena, and 2) the key processes involved in these interactions. These phenomena include large‐scale features such as fronts, high‐ and low‐pressure systems, low‐level jets, and cloud bands, but may also include smaller‐scale phenomena such as mesoscale instabilities (e.g., symmetric instability or slantwise convection), transient meso‐lows and gravity waves. Dynamic air‐ice interactions include direct atmospheric stress and atmospheric stress divergence/deformation; thermodynamic air‐ice interactions involve the net surface energy budget, especially in association with non‐stratocumulus clouds generally occurring with stronger dynamical forcing.
58 GEOSCIENCES↗
An efficient method to estimate the probability density of seismic Green's functions
We present a computationally efficient method to approximate the probability distribution of seismic Green's functions given the uncertainty of an Earth model. The method is based on the Karhunen-Loève (KL) theorem and an approximation of the Green's function (or seismogram) covariance. Using Monte Carlo (MC) simulations as a control case, we demonstrate that our KL-based method can accurately reproduce a probability distribution of seismograms that results from an uncertain Earth model for a MC-derived seismogram covariance. We then describe a method to estimate the covariance of the seismograms resulting from those Earth models that is not based on MC simulations. We use the estimated Green's function covariance in conjunction with our KL-based method to produce a Green's function probability distribution, and compare that distribution to a Green's function probability distribution produced using a MC finite difference method. We find that the Green's function probability distribution approximated using our KL-based method generally mimics that produced using the MC simulations, especially for direct-arriving body waves. However the accuracy of the KL-based method generally decreases for later times in the simulated Green's function distribution.
58 GEOSCIENCES↗