Search NASA⌕ Search

SEARCH · Search NASA

Results for “Multiscale image analysis”

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 73 records · Page 4

Modeling the Self-organized Critical Behavior of Earth's Plasma Sheet Reconnection Dynamics

Analyses of Polar UVI auroral image data show that bright night-side high-latitude W emissions exhibit so many of the key properties of systems in self-organized criticality that an alternate interpretation has become virtually impossible. These analyses will be reviewed. It is now necessary to find and model the source of this behavior. We note that the most common models of self-organized criticality are numerical sandpiles. These are, at root, models that govern the transport of some quantity from a region where it is loaded to another where it is unloaded. Transport is enabled by the excitation of a local threshold instability; it is intermittent and bursty, and it exhibits a number of scale-free statistical properties. Searching for a system in the magnetosphere that is analogous and that, in addition, is known to produce auroral signatures, we focus on the reconnection dynamics of the magnetotail plasma sheet. In our previous work, a driven reconnection model has been constructed and has been under study. The transport of electromagnetic (primarily magnetic) energy carried by the Poynting flux into the reconnection region of the model has been examined. All of the analysis techniques (and more) that have been applied to the auroral image data have also been applied to this Poynting flux. New results will be presented showing that this model also exhibits so many of the key properties of systems in self-organized criticality that an alternate interpretation is implausible. A strong correlation between these key properties of the model and those of the auroral UV emissions will be demonstrated. We suggest that, in general, the driven reconnection model is an important step toward a realistic plasma physical model of self-organized criticality and we conclude, more specifically, that it is also a step in the right direction toward modeling the multiscale reconnection dynamics of the magnetotail.

Klimas, Alexander J.↗

The Atlanta Urban Heat Island Mitigation and Air Quality Modeling Project: How High-Resoution Remote Sensing Data Can Improve Air Quality Models

The Atlanta Urban Heat Island and Air Quality Project had its genesis in Project ATLANTA (ATlanta Land use Analysis: Temperature and Air quality) that began in 1996. Project ATLANTA examined how high-spatial resolution thermal remote sensing data could be used to derive better measurements of the Urban Heat Island effect over Atlanta. We have explored how these thermal remote sensing, as well as other imaged datasets, can be used to better characterize the urban landscape for improved air quality modeling over the Atlanta area. For the air quality modeling project, the National Land Cover Dataset and the local scale Landpro99 dataset at 30m spatial resolutions have been used to derive land use/land cover characteristics for input into the MM5 mesoscale meteorological model that is one of the foundations for the Community Multiscale Air Quality (CMAQ) model to assess how these data can improve output from CMAQ. Additionally, land use changes to 2030 have been predicted using a Spatial Growth Model (SGM). SGM simulates growth around a region using population, employment and travel demand forecasts. Air quality modeling simulations were conducted using both current and future land cover. Meteorological modeling simulations indicate a 0.5 C increase in daily maximum air temperatures by 2030. Air quality modeling simulations show substantial differences in relative contributions of individual atmospheric pollutant constituents as a result of land cover change. Enhanced boundary layer mixing over the city tends to offset the increase in ozone concentration expected due to higher surface temperatures as a result of urbanization.

Quattrochi, Dale A.↗

Effective Li-Ion Transport Quantification in Composite Cathodes for All-Solid-State Batteries via Multiscale Modeling and Experiments

The tortuosity factor of composite cathodes significantly affects the rate performance of all-solid-state batteries (ASSBs) and has significant differences from systems with liquid electrolytes. Here, in this work, we report a simulation-experiment combined approach that quantifies the effective Li-ion transport in an ASSB composite cathode, which links tortuosity factor on ∼ μm scale to terminal voltage during cycling at the cell level (on ∼ cm scale). Two independent approaches of tortuosity factor quantification are considered: fitting electrochemical cycling data and verifying at different cycling rates and calculating from segmented tomography images, with the tortuosity factor quantified from both methods reaching self-consistency. The simulated terminal voltage using the quantified tortuosity factor has a small relative error of <3% compared to the experimental measurements. We find a significantly reduced value of the Bruggeman exponent of the catholyte phase (1.75), and using shape analysis, we show that rod-shaped catholyte particles play an important role in lowering the tortuosity factor.

Yao, Archie Mingze [Univ. of Michigan, Ann Arbor, ↗

Detecting and Characterizing Fracture Zones Using a Convolutional Neural Network

This project directly supports the Geothermal Technologies Office (GTO) objectives outlined in the Multi-Year Program Plan (MYPP) by advancing two key research areas: “Exploration and Characterization” and “Data, Modeling, and Analysis.” This project has successfully demonstrated a pre-drilling ability to image and characterize the distribution and connectivity of subsurface faults and fractures, key parameters for identifying permeable pathways that enable geothermal fluids to circulate and produce energy. Specifically, we developed and implemented innovative machine learning methodologies to enhance geothermal exploration. Large-scale faults were detected using a Convolutional Neural Network (CNN), while small-scale fractures were characterized using a novel Double-Beam Neural Network (DBNN). These tools have proven both technically effective and cost-efficient by reducing reliance on expensive exploratory drilling. Through collaboration with our geothermal industry partner, this research has significantly advanced techniques for identifying hidden geothermal systems and extending the productive lifespan of existing geothermal fields. We applied our methods to two geothermal fields—Soda Lake (Nevada) and Lightning Dock (New Mexico)—to identify shallow steam-charged fracture zones and characterize deep faults at depths of 1.5-2 km. The steam zone identified at the Soda Lake geothermal field showed excellent agreement with prior drilling data, validating the effectiveness of our approaches. In addition, the analysis revealed three new prospective drilling targets for further development and verification. The outcomes of this project improve our scientific understanding of geothermal reservoir behavior, enhance exploration efficiency, extend the economic life of existing geothermal plants. Ultimately, these advancements contribute to GTO’s goal of achieving more sustainable, affordable, and data-driven geothermal energy development across the United States.

15 GEOTHERMAL ENERGY↗

Influence of microstructure on replacement and porosity generation during experimental dolomitization of limestones

Replacement reactions commonly alter the multiscale pore structures of rocks during fluid-rock interactions. Analysis of these processes in various model fluid-rock systems during controlled laboratory experiments provides insights into the origins of microstructures found in natural materials. This study focused on understanding the effects of initial starting material permeability and resultant differences in transport pathways on porosity and mineralogical changes during limestone dolomitization.A series of replacement experiments (32–317 days in duration) have been conducted in which 1.59 cm (5/8 in.) diameter cores of two different limestones were reacted with saturated MgCl2 solutions at 200 °C. The Texas Cream (Austin Chalk) is a high-porosity, high-permeability limestone, whereas both the porosity and permeability of the Carthage Marble (Burlington Limestone) are relatively low. Altered limestones were imaged using scanning electron microscopy with energy-dispersive X-ray spectroscopy (SEM-EDX), Time-of-Flight Secondary Ion Mass Spectrometry (ToF–SIMS) and electron microprobe analysis (EMPA). A representative grain boundary of the low-porosity limestone was targeted for a focused ion beam (FIB) lift-out and characterized using transmission electron microscopy (TEM). These results were coupled with analyses of radial changes in the porosity distribution of the core derived from X-ray and neutron small- and ultra-small angle scattering ((U)SANS/(U)SAXS).The high-porosity/permeability limestone showed a four times faster bulk replacement rate than the lower-porosity/permeability material, and a different mechanism of porosity development. For the low-porosity limestone, a two-stage replacement occurred, with the reacted region of the core consisting of an inner rim in which the limestone was replaced by two calcite-dolomite solid solutions, and an outer rim in which the dolomite was replaced by magnesite. Elongated pores formed along grain boundaries at the initial limestone/dolomite reaction interface, and additional nanometer-scale porosity was formed at the secondary magnesite replacement rim. Grain boundaries were identified as preferential pathways for transport leading to dolomitization and a grain boundary diffusion rate was calculated based on microstructural characterization. In contrast, replacement in the high-porosity limestone was accompanied by porosity generation through replacement of individual grains by dolomite throughout the sample and, in longer runs, magnesite in outer parts of the core. These observations emphasize that both the mechanisms of the replacement reaction and the microstructure and chemistry of the replaced product are contingent on the initial structure of the starting material.

36 MATERIALS SCIENCE↗

New geomorphic data on the active Taiwan orogen: A multisource approach

A multisource and multiscale approach of Taiwan morphotectonics combines different complementary geomorphic analyses based on a new elevation model (DEM), side-looking airborne radar (SLAR), and satellite (SPOT) imagery, aerial photographs, and control from independent field data. This analysis enables us not only to present an integrated geomorphic description of the Taiwan orogen but also to highlight some new geodynamic aspects. Well-known, major geological structures such as the Longitudinal Valley, Lishan, Pingtung, and the Foothills fault zones are of course clearly recognized, but numerous, previously unrecognized structures appear distributed within different regions of Taiwan. For instance, transfer fault zones within the Western Foothills and the Central Range are identified based on analyses of lineaments and general morphology. In many cases, the existence of geomorphic features identified in general images is supported by the results of geological field analyses carried out independently. In turn, the field analyses of structures and mechanisms at some sites provide a key for interpreting similar geomorphic featues in other areas. Examples are the conjugate pattern of strike-slip faults within the Central Range and the oblique fold-and-thrust pattern of the Coastal Range. Furthermore, neotectonic and morphological analyses (drainage and erosional surfaces) has been combined in order to obtain a more comprehensive description and interpretation of neotectonic features in Taiwan, such as for the Longitudinal Valley Fault. Next, at a more general scale, numerical processing of digital elevation models, resulting in average topography, summit level or base level maps, allows identification of major features related to the dynamics of uplift and erosion and estimates of erosion balance. Finally, a preliminary morphotectonic sketch map of Taiwan, combining information from all the sources listed above, is presented.

Deffontaines, B.↗

Analysis of Ultra High Resolution Sea Surface Temperature Level 4 Datasets

Sea surface temperature (SST) studies are often focused on improving accuracy, or understanding and quantifying uncertainties in the measurement, as SST is a leading indicator of climate change and represents the longest time series of any ocean variable observed from space. Over the past several decades SST has been studied with the use of satellite data. This allows a larger area to be studied with much more frequent measurements being taken than direct measurements collected aboard ship or buoys. The Group for High Resolution Sea Surface Temperature (GHRSST) is an international project that distributes satellite derived sea surface temperatures (SST) data from multiple platforms and sensors. The goal of the project is to distribute these SSTs for operational uses such as ocean model assimilation and decision support applications, as well as support fundamental SST research and climate studies. Examples of near real time applications include hurricane and fisheries studies and numerical weather forecasting. The JPL group has produced a new 1 km daily global Level 4 SST product, the Multiscale Ultrahigh Resolution (MUR), that blends SST data from 3 distinct NASA radiometers: the Moderate Resolution Imaging Spectroradiometer (MODIS), the Advanced Very High Resolution Radiometer (AVHRR), and the Advanced Microwave Scanning Radiometer ? Earth Observing System(AMSRE). This new product requires further validation and accuracy assessment, especially in coastal regions.We examined the accuracy of the new MUR SST product by comparing the high resolution version and a lower resolution version that has been smoothed to 19 km (but still gridded to 1 km). Both versions were compared to the same data set of in situ buoy temperature measurements with a focus on study regions of the oceans surrounding North and Central America as well as two smaller regions around the Gulf Stream and California coast. Ocean fronts exhibit high temperature gradients (Roden, 1976), and thus satellite data of SST can be used in the detection of these fronts. In this case, accuracy is less of a concern because the primary focus is on the spatial derivative of SST. We calculated the gradients for both versions of the MUR data set and did statistical comparisons focusing on the same regions.

sea surface temperature SST)↗

Multiscale spatial analysis of fracture arrangement and pattern reconstruction using Ripley's K-function

This work presents novel multiscale spatial data analytics using Ripley's K-function, as a measure of spatial interaction, to study one-dimensional arrangement of fractures. Fracture spatial arrangements are classified into clustered, anticlustered, or indistinguishable from random by testing statistical significance of the calculated Ripley's K-function. Characterizations of fracture arrangements are performed as a function of length scale and position. Analysis of the K-function along the study interval identifies where fracture clustering and anticlustering occur. Here, a simulation technique is also introduced here to statistically reconstruct spatial arrangements and to generate fracture realizations that are spatially similar to the fractures observed in the field. With this simulation technique, one can also fill spatial gaps in fracture measurements where data are absent, unreliable, or unused. Synthetic as well as field-measured 1D fracture datasets are used for testing and demonstration. Methods introduced in this work can be readily applied to fracture datasets observed in outcrops, borehole image logs, and cores.

58 GEOSCIENCES↗

Agile Multi-Scale Decompositions for Automatic Image Registration

In recent works, the first and third authors developed an automatic image registration algorithm based on a multiscale hybrid image decomposition with anisotropic shearlets and isotropic wavelets. This prototype showed strong performance, improving robustness over registration with wavelets alone. However, this method imposed a strict hierarchy on the order in which shearlet and wavelet features were used in the registration process, and also involved an unintegrated mixture of MATLAB and C code. In this paper, we introduce a more agile model for generating features, in which a flexible and user-guided mix of shearlet and wavelet features are computed. Compared to the previous prototype, this method introduces a flexibility to the order in which shearlet and wavelet features are used in the registration process. Moreover, the present algorithm is now fully coded in C, making it more efficient and portable than the MATLAB and C prototype. We demonstrate the versatility and computational efficiency of this approach by performing registration experiments with the fully-integrated C algorithm. In particular, meaningful timing studies can now be performed, to give a concrete analysis of the computational costs of the flexible feature extraction. Examples of synthetically warped and real multi-modal images are analyzed.

Science Data Processing↗

Agile Multi-Scale Decompositions for Automatic Image Registration

In recent works, the first and third authors developed an automatic image registration algorithm based on a multiscale hybrid image decomposition with anisotropic shearlets and isotropic wavelets. This prototype showed strong performance, improving robustness over registration with wavelets alone. However, this method imposed a strict hierarchy on the order in which shearlet and wavelet features were used in the registration process, and also involved an unintegrated mixture of MATLAB and C code. In this paper, we introduce a more agile model for generating features, in which a flexible and user-guided mix of shearlet and wavelet features are computed. Compared to the previous prototype, this method introduces a flexibility to the order in which shearlet and wavelet features are used in the registration process. Moreover, the present algorithm is now fully coded in C, making it more efficient and portable than the MATLAB and C prototype. We demonstrate the versatility and computational efficiency of this approach by performing registration experiments with the fully-integrated C algorithm. In particular, meaningful timing studies can now be performed, to give a concrete analysis of the computational costs of the flexible feature extraction. Examples of synthetically warped and real multi-modal images are analyzed.

Science Data Processing↗

SQMS Quantum R&D in Machine Learning, Optimization and Sensing beyond Fundamental Physics Applications

This newly formed team at SQMS under the Ecosystem Thrust is looking to develop capabilities impacting societal advances outside the core domain of HEP and condensed matter physics. We explicitly leverage the experimental and algorithmic innovations developed across all groups as well as connect to broad-scope external projects of the diverse team of PIs. As the inaugural set of projects, we are studying numerically quantum machine learning models inspired by efficiently trainable echo-state and orthogonal neural networks and developing designs for related experiments to be performed on quantum processors based on SQMS SRF cQED technology and Rigetti s transmon arrays. Investigated models exploit ideas and lessons learned from multiple prior work by SQMS team members in a variety of internal and external activities [R1]. Target initial applications include noisy signal processing, potentially captured by quantum sensors or noisy QPUs, as well as simulation and classification of healthcare data. For instance, image reconstruction of the brain s electrical properties by solving the inverse Maxwell equation problem with uncertainty [R2] through a hybrid quantum-classical physics-informed architecture for time-dependent processes [R3]. The group is also investigating the application and development of novel quantum sensors based on magnetic levitation of a superconducting sphere coupled to a superconducting qubit. This coupling enables high-precision measurements of the position of the sphere, which can be used for sensitive detection of forces, enabling practical applications such as gravimetry for geophysics analysis, or accelerometry for GPS-denied navigation [R4] [R1] Rieffel, Eleanor G., Ata Akbari Asanjan, M. Sohaib Alam, Namit Anand, David E. Bernal Neira, Sophie Block, Lucas T. Brady et al. "Assessing and advancing the potential of quantum computing: A NASA case study." Future Generation Computer Systems (2024). [R2] Yu, X., Serrall s, J.E., Giannakopoulos, I.I., Liu, Z., Daniel, L., Lattanzi, R. and Zhang, Z., 2023. Pifon-ept: Mr-based electrical property tomography using physics-informed fourier networks. IEEE Journal on Multiscale and Multiphysics Computational Techniques. [R3] Wudarski, Filip, Daniel OConnor, Shaun Geaney, Ata Akbari Asanjan, Max Wilson, Elena Strbac, P. Aaron Lott, and Davide Venturelli. "Hybrid quantum-classical reservoir computing for simulating chaotic systems." arXiv preprint arXiv:2311.14105 (2023). [R4] Higgins, Gerard, Saarik Kalia, and Zhen Liu. "Maglev for dark matter: Dark-photon and axion dark matter sensing with levitated superconductors." Physical Review D 109.5 (2024): 055024.

Venturelli, Davide↗

Understanding Solution State Conformation and Aggregate Structure of Conjugated Polymers via Small Angle X-ray Scattering

Donor-acceptor (D-A) conjugated polymers are high-performance organic electronic materials that exhibit complex aggregation behavior. Understanding the solution state conformation and aggregation of conjugated polymers is crucial for controlling morphology during thin-film deposition and the subsequent electronic performance. However, a precise multiscale structure of solution state aggregates is lacking. Here, we present an in-depth small-angle X-ray scattering (SAXS) analysis of the solution state structure of an isoindigo-bithiophene-based D-A polymer (PII-2T) in chlorobenzene and decane as our primary system. Modeling the system as a combination of hierarchical fibrillar aggregates mixed with dispersed polymers, we extract information about conformation and multiscale aggregation and also clarify the physical origin of features often observed but unaddressed or misinterpreted in small-angle scattering patterns of conjugated polymers. The persistence length of the D-A polymer extracted from SAXS agrees well with a theoretical model based on the dihedral potentials. Additionally, we show that the broad high q structure factor peak seen in scattering profiles can be attributed to lamellar stacking occurring within the fibril aggregates and that the low q aggregate scattering is strongly influenced by the polymer molecular weight. Overall, the SAXS profiles of D-A polymers in general exhibit a sensitive dependence on the co-existence of fibrillar aggregate and dispersed polymer chain populations. We corroborate our findings from SAXS with electron microscopy of freeze-dried samples for direct imaging of fibrillar aggregates. Finally, we demonstrate the generality of our approach by fitting the scattering profiles of a variety of D-A polymers based on thieno-isoindigo (PTII-2T), diketopyrrolopyrrole (DPP2T-TT, DPP-BTZ, PDPP2FT-C- 16 ), naphthalenediimide (P(NDI2OD-T2)), and a conjugated block copolymer P3HT-b-DPPT-T. The results presented here establish a picture of the D-A polymer solution state structure and provide a general method of interpreting and analyzing their scattering profiles.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Multiscale Data-Driven Seismic Full-Waveform Inversion With Field Data Study

Seismic full-waveform inversion (FWI), which uses iterative methods to estimate high-resolution subsurface models from seismograms, is a powerful imaging technique in exploration geophysics. In recent years, the computational cost of FWI has grown exponentially due to the increasing size and resolution of seismic data. Moreover, it is a nonconvex problem and can encounter local minima due to the limited accuracy of the initial velocity models or the absence of low frequencies in the measurements. To overcome these computational issues, we develop a multiscale data-driven FWI method based on fully convolutional networks (FCNs). In preparing the training data, we first develop a real-time style transform method to create a large set of synthetic subsurface velocity models from natural images. We then develop two convolutional neural networks with encoder-decoder structures to reconstruct the low- and high-frequency components of the subsurface velocity models, separately. To validate the performance of our data-driven inversion method and the effectiveness of the synthesized training set, we compare it with conventional physics-based waveform inversion approaches using both synthetic and field data. Finally, these numerical results demonstrate that, once our model is fully trained, it can significantly reduce the computation time and yield more accurate subsurface velocity models in comparison with conventional FWI.

58 GEOSCIENCES↗

Fractal Characterization of Multimodal, Multiscale Images of Shale Rock Fracture Networks

An array of multimodal and multiscale images of fractured shale rock samples and analogs was collected with the aim of improving the numerical representation of fracture networks. 2D and 3D reconstructions of fractures and matrices span from 10 –6 to 100 m. The origin of the fracture networks ranged from natural to thermal maturation to hydraulically induced maturation. Images were segmented to improve fracture identification. Then, the fractal dimension and length distribution of the fracture networks were calculated for each image dataset. The resulting network connectivity and scaling associations are discussed at length on the basis of scale, sample and order of magnitude. Fracture network origin plays an important role in the resulting fracture systems and their scaling. Rock analogs are also evaluated using these descriptive tools and are found to be faithful depictions of maturation-induced fracture networks.

58 GEOSCIENCES↗

Proceedings of the 20th International Symposium on Space Flight Dynamics

Topics include: Measuring Image Navigation and Registration Performance at the 3-Sigma Level Using Platinum Quality Landmarks; Flight Dynamics Performances of the MetOp A Satellite during the First Months of Operations; Visual Navigation - SARE Mission; Determining a Method of Enabling and Disabling the Integral Torque in the SDO Science and Inertial Mode Controllers; Guaranteeing Pointing Performance of the SDO Sun-Pointing Controllers in Light of Nonlinear Effects; SDO Delta H Mode Design and Analysis; Observing Mode Attitude Controller for the Lunar Reconnaissance Orbiter; Broken-Plane Maneuver Applications for Earth to Mars Trajectories; ExoMars Mission Analysis and Design - Launch, Cruise and Arrival Analyses; Mars Reconnaissance Orbiter Aerobraking Daily Operations and Collision Avoidance; Mars Reconnaissance Orbiter Interplanetary Cruise Navigation; Motion Parameters Determination of the SC and Phobos in the Project Phobos-Grunt; GRAS NRT Precise Orbit Determination: Operational Experience; Orbit Determination of LEO Satellites for a Single Pass through a Radar: Comparison of Methods; Orbit Determination System for Low Earth Orbit Satellites; Precise Orbit Determination for ALOS; Anti-Collision Function Design and Performances of the CNES Formation Flying Experiment on the PRISMA Mission; CNES Approaching Guidance Experiment within FFIORD; Maneuver Recovery Analysis for the Magnetospheric Multiscale Mission; SIMBOL-X: A Formation Flying Mission on HEO for Exploring the Universe; Spaceborne Autonomous and Ground Based Relative Orbit Control for the TerraSAR-X/TanDEM-X Formation; First In-Orbit Experience of TerraSAR-X Flight Dynamics Operations; Automated Target Planning for FUSE Using the SOVA Algorithm; Space Technology 5 Post-Launch Ground Attitude Estimation Experience; Standardizing Navigation Data: A Status Update; and A Study into the Method of Precise Orbit Determination of a HEO Orbiter by GPS and Accelerometer.

Woodard, Mark↗

An Integrated Multiscale Experimental-Numerical Analysis on Reconsolidation of Salt-Clay Mixture for Disposal of Heat-Generating Waste (Final NEUP Technical Report)

The overall purpose of this research is to improve understanding of THMC coupling effect on the reconsolidation of granular (or crushed) salt-clay mixture used for seal systems of shafts and drifts in salt repositories. This proposed work is partially motivated by the recent work on the Waste Isolation Pilot Plant (WIPP) that shows the promising sealing capability of clay-salt mixture compared to crushed salt. In particular, primary emphasis is to develop a fully integrated multiscale experiment-numerical study to determine and explain what leads to the superior sealing ability of the clay-salt mixture. These research activities are designed to seek further understanding of (1) why clay additives may enhance the fluid trapping and (2) whether this flow barrier effect may prevail under different combinations of temperature, confining pressure, deviatoric stress and other foreseeable environmental factors. If successful, this enhanced flow trapping ability of the seal provides significant improvement to the seal and repository performance and therefore make the repository safer in the long-term. The experiment component includes microstructural investigation and macroscopic tests on a reconsolidated salt-clay mixture. In the microstructural study, the goal is to (1) characterize microscopic distributions of distinct phases (e.g., clay, salt crystal boundaries, trapped brine, and pore) to examine the connectivity of the pore network inside the salt-clay mixture with different amounts of clay additive and moisture content and (2) analyze multiscale imaging data to reconstruct the polycrystalline microstructures for numerical simulations. Meanwhile, macroscopic tests are performed to analyze how clay alters the failure/creep mechanisms in the salt-clay mixture. Microscopic and macroscopic experimental observations will both be used to calibrate and validate a multiscale model that explicitly simulates the capillary and multiphase flow in the connected pores and the deformation due to the presence of intra-crystalline brine at the pore scale via a new polyhedral discrete element–lattice Boltzmann method (DEM-LBM) coupling model. The pore-scale simulations are homogenized via an upscaling procedure that converts pore-scale information (e.g. force exerted on grain boundary, sliding, pressure-solution) to continuum measures (e.g. Cauchy stress, Darcy’s flow) at each integration point in the macroscopic multiphase TMHC model. This multiscale scheme will allow coupling be- tween high-fidelity simulations of brine-salt-clay interaction and the macroscopic TMHC model. The multiscale model helps the understanding of how the trapped brine inclusion affects the pressure-solution mechanism with the presence of clay and moisture. This work brings new insight into the sealing capacity of salt-clay mixture under elevated temperature over a long period of time - a key to evaluating the potential of salt-clay mixture usage for salt repositories.

42 ENGINEERING↗

Fast Plasma Instrument for MMS: Data Compression Simulation Results

Magnetospheric Multiscale (MMS) mission will study small-scale reconnection structures and their rapid motions from closely spaced platforms using instruments capable of high angular, energy, and time resolution measurements. To meet these requirements, the Fast Plasma Instrument (FPI) consists of eight (8) identical half top-hat electron sensors and eights (8) identical ion sensors and an Instrument Data Processing Unit (IDPU). The sensors (electron or ion) are grouped into pairs whose 6 deg x 180 deg fields-of-view (FOV) are set 90 deg apart. Each sensor is equipped with electrostatic aperture steering to allow the sensor to scan a 45 deg x 180 deg fan about its nominal viewing (0 deg deflection) direction. Each pair of sensors, known as the Dual Electron Spectrometer (DES) and the Dual Ion Spectrometer (DIS), occupies a quadrant on the MMS spacecraft and the combination of the eight electron/ion sensors, employing aperture steering, image the full-sky every 30-ms (electrons) and 150-ms (ions), respectively. To probe the results in the DES complement of a given spacecraft generating 6.5-Mbs(exp -1) of electron data while the DIS generates 1.1-Mbs(exp -1) of ion data yielding an FPI total data rate of 6.6-MBs(exp -1). The FPI electron/ion data is collected by the IDPU then transmitted to the Central Data Instrument Processor (CIDP) on the spacecraft for science interest ranking. Only data sequences that contain the greatest amount of temporal/spatial structure will be intelligently down-linked by the spacecraft. Currently, the FPI data rate allocation to the CIDP is 1.5-Mbs(exp -1). Consequently, the FPI-IDPU must employ data/image compression to meet this CIDP telemetry allocation. Here, we present simulations of the CCSDS 122.0-B-1 algorithm-based compression of the FPI-DES electron data. Compression analysis is based upon a seed of re-processed Cluster/PEACE electron measurements. Topics to be discussed include: review of compression algorithm; data quality; data formatting/organization; and, implications for data/matrix pruning. To conclude a presentation of the base-lined FPI data compression approach is provided.

Barrie, A.↗

Internal structure of the San Jacinto fault zone at the Ramona Reservation, north of Anza, California, from dense array seismic data

SUMMARY We image the internal structure of the San Jacinto fault zone (SJFZ) near Anza, California, with seismic data recorded by two dense arrays (RA and RR) from ∼42 000 local and ∼180 teleseismic events occurring between 2012 and 2017. The RA linear array has short aperture (∼470 m long with 12 strong motion sensors) and recorded for the entire analysed time window, whereas the RR is a large three-component nodal array (97 geophones across a ∼2.4 km × 1.4 km area) that operated for about a month in September–October 2016. The SJFZ at the site contains three near-parallel surface traces F1, F2 and F3 from SW to NE that have accommodated several Mw > 6 earthquakes in the past 15 000 yr. Waveform changes in the fault normal direction indicate structural discontinuities that are consistent with the three fault surface traces. Relative slowness from local events and delay time analysis of teleseismic arrivals in the fault normal direction suggest a slower SW side than the NE with a core damage zone between F1 and F2. This core damage zone causes ∼0.05 s delay at stations RR26–31 in the teleseismic P arrivals compared with the SW-most station, and generates both P- and S-type fault zone trapped waves. Inversion of S trapped waves indicates the core damaged structure is ∼100 m wide, ∼4 km deep with a Q value of ∼20 and 40 per cent S-wave velocity reduction compared with bounding rocks. Fault zone head waves observed at stations SW of F3 indicate a local bimaterial interface that separates the locally faster NE block from the broad damage zone in the SW at shallow depth and merges with a deep interface that separates the regionally faster NE block from rocks to the SW with slower velocities at greater depth. The multiscale structural components observed at the site are related to the geological and earthquake rupture history at the site, and provide important information on the preferred NW propagation of earthquake ruptures on the San Jacinto fault.

Geochemistry & Geophysics↗