Deep Spatiotemporal Clustering: A Temporal Clustering Approach for Multi-dimensional Climate Data
Explore the source record for details and available documents.
Geosciences↗
SEARCH · Search NASA
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.
Measurement of the phase between strong and electromagnetic amplitudes in the decay J/ψ → ϕη
The first direct measurement of the relative phase between the strong and electromagnetic amplitudes for a J/ψ decaying into a vector-pseudoscalar final state is performed using 26 energy points of e + e − annihilation data between 3.00 GeV and 3.12 GeV. The data sets were collected by the BESIII detector with a total integrated luminosity of 452 pb −1 . By investigating the interference pattern in the cross section lineshape of e + e − → ϕη, the relative phase between the strong and electromagnetic amplitudes of J/ψ decay is determined to be within [133°, 228°] at 68% confidence level.
QCD↗
Updates to the Regional Seismic Travel Time (RSTT) Model: 1. Tomography
Abstract A function of global monitoring of nuclear explosions is the development of Earth models for predicting seismic travel times for more accurate calculation of event locations. Most monitoring agencies rely on fast, distance-dependent one-dimensional (1D) Earth models to calculate seismic event locations quickly and in near real-time. RSTT (Regional Seismic Travel Time) is a seismic velocity model and computer software package that captures the major effects of three-dimensional crust and upper mantle structure on regional seismic travel times, while still allowing for fast prediction speed (milliseconds). We describe updates to the RSTT model using a refined data set of regional phases (i.e., Pn, Pg, Sn, Lg) using the Bayesloc relative relocation algorithm. The tomographic inversion shown here acts to refine the previous RSTT public model ( rstt201404um ) and displays significant features related to areas of global tectonic complexity as well as further reduction in arrival residual values. Validation of the updated RSTT model demonstrates significant reduction in median epicenter mislocation (15.3 km) using all regional phases compared to the iasp91 1D model (22.1 km) as well as to the current station correction approach used at the Comprehensive Nuclear-Test-Ban Treaty Organization International Data Centre (18.9 km).
58 GEOSCIENCES↗
Updates to the Regional Seismic Travel Time (RSTT) Model: 2. Path-dependent Travel-time Uncertainty
Abstract The regional seismic travel time (RSTT) model and software were developed to improve travel-time prediction accuracy by accounting for three-dimensional crust and upper mantle structure. Travel-time uncertainty estimates are used in the process of associating seismic phases to events and to accurately calculate location uncertainty bounds (i.e. event location error ellipses). We improve on the current distance-dependent uncertainty parameterization for RSTT using a random effects model to estimate slowness (inverse velocity) uncertainty as a mean squared error for each model parameter. The random effects model separates the error between observed slowness and model predicted slowness into bias and random components. The path-specific travel-time uncertainty is calculated by integrating these mean squared errors along a seismic-phase ray path. We demonstrate that event location error ellipses computed for a 90% coverage ellipse metric (used by the Comprehensive Nuclear-Test-Ban Treaty Organization International Data Centre (IDC)), and using the path-specific travel-time uncertainty approach, are more representative (median 82.5% ellipse percentage) of true location error than error ellipses computed using distance-dependent travel-time uncertainties (median 70.1%). We also demonstrate measurable improvement in location uncertainties using the RSTT method compared to the current station correction approach used at the IDC (median 74.3% coverage ellipse).
58 GEOSCIENCES↗
Predicting Fracture Network Development in Crystalline Rocks
The geometric properties of fractures influence whether they propagate, arrest, or coalesce with other fractures. Thus, quantifying the relationship between fracture network characteristics may help predict fracture network development, and perhaps precursors to catastrophic failure. To constrain the relationship and predictability of fracture characteristics, we deform eight one centimeter tall rock cores under triaxial compression while acquiring in situ X-ray tomograms. The tomograms reveal precise measurements of the fracture network characteristics above the spatial resolution of 6.5 µm. We develop machine learning models to predict the value of each characteristic using the other characteristics, and excluding the macroscopic stress or strain imposed on the rock. The models predict fracture development more accurately in the experiments performed on granite and monzonite, than the experiments on marble. Fracture network development may be more predictable in these igneous rocks because their microstructure is more mechanically homogeneous than the marble, producing more systematic fracture development that is not strongly impeded by grain contacts and cleavage planes. The varying performance of the models suggest that fracture volume, length, and aperture are the most predictable of the characteristics, while fracture orientation is the least predictable. Orientation does not correlate with length, as suggested by the idea that the orientation evolves with increasing differential stress and thus fracture length. This difference between the observed and expected relationship between orientation and length highlights the influence of mechanical heterogeneities and local stress perturbations on fracture growth as fractures propagate, link, and coalesce.
58 GEOSCIENCES↗
Developing a Consistent Travel-Time Framework for Comparing Three-Dimensional Velocity Models for Seismic Location Accuracy
Abstract Location algorithms have historically relied on simple, one-dimensional (1D) velocity models for fast seismic event locations. 1D models are generally used as travel-time lookup tables, one for each seismic phase, with travel-times pre-calculated for event distance and depth. These travel-time lookup tables are extremely fast to use and this fast computational speed makes them the preferred type of velocity model for operational needs. Higher-dimensional (i.e., three-dimensional—3D) seismic velocity models are becoming readily available and provide more accurate event locations over 1D models. The computational requirements of these 3D models tend to make their operational use prohibitive. Additionally, comparing location accuracy for 3D seismic velocity models tends to be problematic, as each model is determined using different ray-tracing algorithms. Attempting to use a different algorithm than the one used to develop a model usually results in poor travel-time prediction. We demonstrate and test a framework to create first-P and first-S 3D travel-time correction surfaces using an open-source framework ( PCalc + GeoTess , https://www.sandia.gov/salsa3d/software/geotess ) that easily stores 3D travel-time and uncertainty data. This framework produces fast travel-time and uncertainty predictions and overcomes the ray-tracing algorithm hurdle because the lookup tables can be generated using the exact ray-tracing algorithm that is preferred for a model.
3D velocity models↗
3D Time-Lapse Electrical Resistivity Imaging of Rock Damage Patterns and Gas Flow Paths Resulting from Two Underground Chemical Explosions
Abstract Rock damage from underground nuclear explosions (UNEs) has a strong influence on sub-surface gas movement and on seismic waveform characteristics, both of which are used to detect UNEs. Although advanced numerical simulation capabilities exist to predict rock damage patterns and corresponding detection signals, those predictions are dependent on (generally) unknown properties of the host rock. For example, the effects of in-situ mechanical heterogeneities on the explosively generated damage/fractures that provide gas flow pathways to the surface are not well understood, due largely to the difficulty in accessing and characterizing the near-source region. In this paper we demonstrate the emerging use of electrical resistivity tomography (ERT) for imaging rock damage and gas flow patterns resulting from two relatively small-scale underground chemical explosions. Pre-explosion ERT and crosshole seismic imaging revealed a natural fracture zone within the test bed. Post-explosion imaging revealed that the damage zone was non-symmetric and was focused primarily within the pre-existing fracture zone, located 10 m above the first explosion and 5 m above the second explosion. Time-lapse ERT imaging of heated air injected into the detonation borehole revealed the primary gas flow paths to be within the upper margin of the same primary damage zone. These results point to the utility of ERT imaging for understanding rock damage and gas flow patterns under experimental conditions, and to the importance of understanding the effects of geologic heterogeneity on UNE detection signals, particularly gas surface breakthrough times.
58 GEOSCIENCES↗
The Dynamic Networks Experiments: Virtual Experiments to Quantify Gains in Nuclear Explosion Monitoring
We describe an ongoing series of virtual experiments conducted collaboratively by four United States National Laboratories: Sandia National Laboratories, Los Alamos National Laboratory, Lawrence Livermore National Laboratory, and Pacific Northwest National Laboratory. These Dynamic Network Experiments (DNEs) provide an experimental framework to evaluate the potential impact of new research tools on nuclear explosion monitoring. The second DNE (DNE2), completed in 2024, exploited waveform data (seismic, infrasound, and electromagnetic) that was recorded by multi-modal sensors within and near the Nevada National Security Site and synthetic radionuclide signatures over multiple time periods. During the execution of DNE2, we processed and analyzed data through a multi-stage event processing pipeline that ingested raw data, performed quality control, detected signals, built events from these signals, located these events, and characterized the events’ source types and sizes. For each stage and over the entire event processing pipeline, we evaluated performance changes by comparing the performance of new data processing methods, models, and algorithms against a baseline. We also performed an additional execution phase to assess event processing pipeline function, speed, and efficiency against that of an expert analyst, including computational and manual efforts. Finally, we assessed the impact and effort of modern computing infrastructure on the monitoring pipeline. This paper describes key elements of the DNEs, from formulation through execution, as demonstrated in DNE2. The DNEs introduce several novel concepts to quantitatively measure the potential impact of new methods on explosion monitoring, including the collaborative design of multi-modal datasets, performance and logistical metrics, and integrated analyses.
42 ENGINEERING↗
Settling of two-way momentum and energy coupled particles subject to Boussinesq and non-Boussinesq heating
This work establishes a procedure to accurately compute heat transfer between an Eulerian fluid and Lagrangian point-particles. Recent work has focused on accurately computing momentum transfer between fluid and particles. The coupling term for momentum involves the undisturbed fluid velocity at the particle location which is not directly accessible in the simulation. Analogously, in the context of thermal coupling, the undisturbed fluid temperature at the particle location is not directly accessible in simulations and must be estimated. In this paper, we develop a scheme to accurately estimate the undisturbed fluid temperature of a point-particle exchanging thermal energy with a surrounding fluid. Furthermore, the temperature disturbance is correlated with the enhanced temperature curvature in the vicinity of the particle and is formally valid in the low heating, low convection limit. We conduct extensive verification of the correction procedure for a settling particle subject to radiation. This setup allows the simultaneous testing of thermal and momentum corrections. By considering equations of drag and Nusselt number extended to finite Péclet and Boussinesq numbers, we establish a large range over which the correction procedure can be applied.
42 ENGINEERING↗
Conservative numerical schemes with optimal dispersive wave relations: Part I. Derivation and analysis
An energy-conserving and an energy-and-enstrophy conserving numerical schemes are derived by approximating the Hamiltonian formulation of the inviscid shallow water flows based on the vorticity-divergence variables. These schemes also conserve the first-order moments such as mass and vorticity, as usual. The conservative properties of the schemes stem from the skew-symmetry and singularities of the Poisson brackets, which are carefully retained in the discrete approximations. Here, the schemes operate on unstructured orthogonal dual meshes, over bounded or unbounded domains, and they are also shown to possess the same optimal dispersive wave relations as those of the Z-grid scheme, which is a consequence of the use of the vorticity and divergence variables.
54 ENVIRONMENTAL SCIENCES↗
Lattice dynamics, sound velocities, and atomic environments of szomolnokite at high pressure
Complex mixtures of sulfates, silicates, and ice have been observed in a variety of planetary environments on Earth, Mars and the icy satellites of the solar system. Characterizing the properties of the corresponding compositional endmembers is important for understanding the interiors of a range of planetary bodies in which these phases are observed. To measure the electronic and vibrational properties of the pure ferrous iron endmember of the kieserite group, szomolnokite, (FeSO 4 ∙H 2 O), we have performed synchrotron 57 Fe nuclear resonant inelastic and forward scattering experiments in the diamond-anvil cell up to 14.5 GPa. This pressure range covers depths within Earth’s interior relevant to sulfur cycling in subduction zones and the range of pressures expected within icy satellite interiors. We find evidence of crystal lattice softening, changes in elastic properties, and changes in the electric field gradients of iron atoms associated with two structural transitions occurring within the experimental pressure range. Here we apply these findings to icy satellite interiors, including discussion of elastic properties, modeling of ice-sulfate aggregates, and implications for tidal observations.
58 GEOSCIENCES↗
In situ X-ray and IR probes relevant to Earth science at the Advanced Light Source at Lawrence Berkeley Laboratory
Access to synchrotron X-ray facilities has become an important aspect for many disciplines in experimental Earth science. This is especially important for studies that rely on probing samples in situ under natural conditions different from the ones found at the surface of the Earth. The non-ambient condition Earth science program at the Advanced Light Source (ALS), Lawrence Berkeley National Laboratory, offers a variety of tools utilizing the infra-red and hard X-ray spectrum that allow Earth scientists to probe Earth and environmental materials at variable conditions of pressure, stress, temperature, atmospheric composition, and humidity. These facilities are important tools for the user community in that they offer not only considerable capacity (non-ambient condition diffraction) but also complementary (IR spectroscopy, microtomography), and in some cases unique (Laue microdiffraction) instruments. The availability of the ALS’ in situ probes to the Earth science community grows especially critical during the ongoing dark time of the Advanced Photon Source in Chicago, which massively reduces available in situ synchrotron user time in North America.
58 GEOSCIENCES↗
A data-driven framework for permeability prediction of natural porous rocks via microstructural characterization and pore-scale simulation
Understanding the microstructure–property relationships of porous media is of great practical significance, based on which macroscopic physical properties can be directly derived from measurable microstructural informatics. However, establishing reliable microstructure–property mappings in an explicit manner is difficult, due to the intricacy, stochasticity, and heterogeneity of porous microstructures. In this paper, a data-driven computational framework is presented to investigate the inherent microstructure–permeability linkage for natural porous rocks, where multiple techniques are integrated together, including microscopy imaging, stochastic reconstruction, microstructural characterization, pore-scale simulation, feature selection, and data-driven modeling. A large number of 3D digital rocks with a wide porosity range are acquired from microscopy imaging and stochastic reconstruction techniques. A broad variety of morphological descriptors are used to quantitatively characterize pore microstructures from different perspectives, and they compose the raw feature pool for feature selection. Here high-fidelity lattice Boltzmann simulations are conducted to resolve fluid flow passing through porous media, from which reliable permeability references are obtained. The optimal feature set that best represents permeability is identified through a performance-oriented feature selection process, upon which a cost-effective surrogate model is rapidly fitted to approximate the microstructure-permeability mapping via data-driven modeling. This surrogate model exhibits great advantages over empirical/analytical formulas in terms of prediction accuracy and generalization capacity, which can predict reliable permeability values spanning four orders of magnitude. Besides, feature selection also greatly enhances the interpretability of the data-driven prediction model, from which new insights into the mechanism of how microstructural characteristics determine intrinsic permeability are obtained.
58 GEOSCIENCES↗
Urbanization Impact on Regional Climate and Extreme Weather: Current Understanding, Uncertainties, and Future Research Directions
Abstract Urban environments lie at the confluence of social, cultural, and economic activities and have unique biophysical characteristics due to continued infrastructure development that generally replaces natural landscapes with built-up structures. The vast majority of studies on urban perturbation of local weather and climate have been centered on the urban heat island (UHI) effect, referring to the higher temperature in cities compared to their natural surroundings. Besides the UHI effect and heat waves, urbanization also impacts atmospheric moisture, wind, boundary layer structure, cloud formation, dispersion of air pollutants, precipitation, and storms. In this review article, we first introduce the datasets and methods used in studying urban areas and their impacts through both observation and modeling and then summarize the scientific insights on the impact of urbanization on various aspects of regional climate and extreme weather based on more than 500 studies. We also highlight the major research gaps and challenges in our understanding of the impacts of urbanization and provide our perspective and recommendations for future research priorities and directions.
54 ENVIRONMENTAL SCIENCES↗
Improving Simulations of Vegetation Dynamics over the Tibetan Plateau: Role of Atmospheric Forcing Data and Spatial Resolution
The efficacy of vegetation dynamics simulations in offline land surface models (LSMs) largely depends on the quality and spatial resolution of meteorological forcing data. In this study, the Princeton Global Meteorological Forcing Data (PMFD) and the high spatial resolution and upscaled China Meteorological Forcing Data (CMFD) were used to drive the Simplified Simple Biosphere model version 4/Top-down Representation of Interactive Foliage and Flora Including Dynamics (SSiB4/TRIFFID) and investigate how meteorological forcing datasets with different spatial resolutions affect simulations over the Tibetan Plateau (TP), a region with complex topography and sparse observations. By comparing the monthly Leaf Area Index (LAI) and Gross Primary Production (GPP) against observations, we found that SSiB4/TRIFFID driven by upscaled CMFD improved the performance in simulating the spatial distributions of LAI and GPP over the TP, reducing RMSEs by 24.3% and 20.5%, respectively. The multi-year averaged GPP decreased from 364.68 gC m –2 yr –1 to 241.21 gC m –2 yr –1 with the percentage bias dropping from 50.2% to –1.7%. When using the high spatial resolution CMFD, the RMSEs of the spatial distributions of LAI and GPP simulations were further reduced by 7.5% and 9.5%, respectively. This study highlights the importance of more realistic and high-resolution forcing data in simulating vegetation growth and carbon exchange between the atmosphere and biosphere over the TP.
54 ENVIRONMENTAL SCIENCES↗
The Spatial and Temporal Variability of the Clear Convective Boundary Layer at the ARM SGP Supersite
The convective boundary layer (CBL), also known as the mixing layer, constitutes the critical lower segment of the atmosphere that significantly influences daily human activities. The growing demand for precise weather forecasts is driven by the requirements of agriculture, transportation, and routine societal functions. Here, to enhance understanding of the CBL, this study investigates the spatiotemporal variability in the CBL and its controlling factors using four-year Doppler lidar, surface flux, and profiling measurements at five ARM Southern Great Plains sites within a 100 km radius. This investigation utilizes data collected exclusively under clear-sky conditions or scattered low-cloud conditions. Results reveal significant spatial differences in CBL evolutions. Daily mixing layer heights (MLHs) vary up to 1 km (30% of the mean) in late afternoon. There is a clear east–west contrast: western sites (C1, E32, E37) exhibit higher summer MLH (1.9–2.1 km) and vertical velocity variances (1.0–1.2 m 2 s −2 ) than eastern sites (1.6–1.8 km), reversing in winter. Temporally, the MLH peaks at 70% of the sunrise–sunset interval, the lagging heat flux (HF) peaks at 50%; and the seasonal MLH maxima lag the HF by approximately one month, influenced by nighttime PBL (planetary boundary layer) properties. The HF and lower tropospheric stability are the main factors of influence for the CBL, but site-specific dependencies highlight the critical roles of local factors, underscoring the need for including them in CBL modeling.
ARM SGP site↗
Detailing cloud property feedbacks with a regime-based decomposition
Here, we report diagnosing the root causes of cloud feedback in climate models and reasons for inter-model disagreement is a necessary first step in understanding their wide variation in climate sensitivities. Here we bring together two analysis techniques that illuminate complementary aspects of cloud feedback. The first quantifies feedbacks from changes in cloud amount, altitude, and optical depth, while the second separates feedbacks due to cloud property changes within specific cloud regimes from those due to regime occurrence frequency changes. We find that in the global mean, shortwave cloud feedback averaged across ten models comes solely from a positive within-regime cloud amount feedback countered slightly by a negative within-regime optical depth feedback. These within-regime feedbacks are highly uniform: In nearly all regimes, locations, and models, cloud amount decreases and cloud albedo increases with warming. In contrast, global-mean across-regime components vary widely across models but are very small on average. This component, however, is dominant in setting the geographic structure of the shortwave cloud feedback: Thicker, more extensive cloud types increase at the expense of thinner, less extensive cloud types in the extratropics, and vice versa at low latitudes. The prominent negative extratropical optical depth feedback has contributions from both within- and across-regime components, suggesting that thermodynamic processes affecting cloud properties as well as dynamical processes that favor thicker cloud regimes are important. The feedback breakdown presented herein may provide additional targets for observational constraints by isolating cloud property feedbacks within specific regimes without the obfuscating effects of changing dynamics that may differ across timescales.
58 GEOSCIENCES↗
Diurnal cycle of precipitation over global monsoon systems in CMIP6 simulations
Accurate simulation of the diurnal cycle of precipitation over different climate regimes, including global monsoon systems, is an ongoing challenge for current state-of-the-art climate models. Here, in this study, the diurnal cycles of summertime precipitation from 21 CMIP6 historical simulations are evaluated with the NASA IMERG data via the diurnal harmonic amplitude and phase from the Fourier analysis in seven regional monsoon systems: Northern and Southern Africa monsoons, South Asia monsoon, East Asia monsoon, Australia monsoon, and North and South America monsoons. Particularly, we focus on four subregions including ocean, coastal ocean, land, and coastal land considering the different underlying physical processes. It is found that the CMIP6 historical simulations perform fairly well on reproducing the diurnal cycle of precipitation over ocean and coastal ocean, while the simulated diurnal harmonic phase and amplitude over land and coastal land show larger biases and model spread compared to biases over ocean and coastal ocean. The simulation of precipitation diurnal cycle is substantially improved in the multi-model mean of both CMIP6 historical and CMIP6 AMIP simulations compared with that of CMIP5 historical simulations. But the model bias of precipitation diurnal cycle is in general similar between CMIP6 AMIP and CMIP6 historical simulations, which suggests that the impact of the interactive ocean on the simulated diurnal precipitation averaged over a large domain is generally minor. The influence of model resolution on the diurnal cycle is shown in mixed directions across various subregions, and significant improvement from the low- to high-resolution models is only notable over coastal ocean.
54 ENVIRONMENTAL SCIENCES↗