Search NASA⌕ Search

SEARCH · Search NASA

Results for “seismic 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 433 records · Page 24

Connecting Risk and Resilience for a Power System Using the Portland Hills Fault Case Study

Active seismic faults in the Pacific Northwest area have encouraged electric utilities in the region to deeply contemplate and proactively intervene to support grid resilience. To further this effort this research introduces Monte Carlo (MC)-based power system modeling as a means to inform the Performance Based Earthquake Engineering method and simulates 100,000 sample earthquakes of a 6.8 magnitude (M6.8) Portland Hills Fault (PHF) scenario in the Portland General Electric (PGE) service territory as a proof of concept. This paper also proposes the resilience metric Seismic Load Recovery Factor (SLRF) to quantify the recovery of a downed power system and thus can be used to quantify earthquake economic risk. Using MC results, the SLRF was evaluated to be 19.7 h and the expected economic consequence cost of a M6.8 PHF event was found to be $\$180$ million with an annualized risk of $\$90,000$ given the event’s 1 in 2000 year probability of occurrence. The MC results also identified the eight most consequential substations in the PGE system—i.e., those that contributed to maximum load loss. This paper concludes that retrofitting these substations reduced the expected consequence cost of a M6.8 PHF event to $117 million.

24 POWER TRANSMISSION AND DISTRIBUTION↗

The Cretaceous-Tertiary Impact Crater and the Cosmic Projectile that Produced it

Evidence gathered to date from topographic data, geophysical data, well logs, and drill-core samples indicates that the buried Chicxulub basin, the source crater for the approximately 65 Ma Cretaceous-Tertiary (K/T) boundary deposits, is approximately 300 km in diameter. A prominent topographic ridge and a ring of gravity anomalies mark the position of the basin rim at approximately 150 km from the center. Wells in this region recovered thick sequences of impact-generated breccias at 200-300 m below present sea level. Inside the rim, which has been severely modified by erosion following impact, the subsurface basin continues to deepen until near the center it is approximately 1 km deep. The best planetary analog for this crater appears to be the 270 km-diameter Mead basin on Venus. Seismic reflection data indicate that the central zone of downward displacement and excavation (the transient crater is approximately 130 km in diameter, consistent with previous studies of gravity anomaly data). Our analysis of projectile characteristics utilizes this information, coupled with conventional scaling relationships, and geochemical constraints on the mass of extraterrestrial material deposited within the K/T boundary layer. Results indicate that the Chicxulub crater would most likely be formed by a long-period comet composed primarily of nonsilicate materials (ice, hydrocarbons, etc.) and subordinate amounts (less than or equal to 50 percent) primitive chondritic material. This collision would have released the energy equivalent to between 4 x 10(exp 8) and 4 x 10(exp 9) megatons of TNT. Studies of terrestrial impact rates suggest that such an event would have a mean production rate of approximately 1.25 x 10(exp -9) y(exp -1). This rate is considerably lower than that of the major mass extinctions over the last 250 million years (approximately 5 x 10(exp -7) y(exp -1). Consequently, while there is substantial circumstantial evidence establishing the cause-effect link between the Chicxulub basin forming event and the K/T biological extinctions, the results of our analysis do not support models of impact as a common or singular causative agent of mass extinctions on Earth.

Sharpton, Virgil L.↗

Time-Lapse Integration at FWU: Fluids, Rock Physics, Numerical Model Integration, and Field Data Comparison

We present the current status of time-lapse seismic integration at the Farnsworth (FWU) CO2 WAG (water-alternating-gas) EOR (Enhanced Oil Recovery) project at Ochiltree County, northwest Texas. As a potential carbon sequestration mechanism, CO2 WAG projects will be subject to some degree of monitoring and verification, either as a regulatory requirement or to qualify for economic incentives. In order to evaluate the viability of time-lapse seismic as a monitoring method the Southwest Partnership (SWP) has conducted time-lapse seismic monitoring at FWU using the 3D Vertical Seismic Profiling (VSP) method. The efficacy of seismic time-lapse depends on a number of key factors, which vary widely from one application to another. Most important among these are the thermophysical properties of the original fluid in place and the displacing fluid, followed by the petrophysical properties of the rock matrix, which together determine the effective elastic properties of the rock fluid system. We present systematic analysis of fluid thermodynamics and resulting thermophysical properties, petrophysics and rock frame elastic properties, and elastic property modeling through fluid substitution using data collected at FWU. These analyses will be framed in realistic scenarios presented by the FWU CO2 WAG development. The resulting fluid/rock physics models will be applied to output from the calibrated FWU compositional reservoir simulation model to forward model the time-lapse seismic response. Modeled results are compared with field time-lapse seismic measurements and strategies for numerical model feedback/update are discussed. While mechanical effects are neglected in the work presented here, complementary parallel studies are underway in which laboratory measurements are introduced to introduce stress dependence of matrix elastic moduli.

4D↗

Pressure-stabilized divalent ozonide CaO 3 and its impact on Earth's oxygen cycles

High pressure can drastically alter chemical bonding and produce exotic compounds that defy conventional wisdom. Especially significant are compounds pertaining to oxygen cycles inside Earth, which hold key to understanding major geological events that impact the environment essential to life on Earth. Here we report the discovery of pressure-stabilized divalent ozonide CaO3 crystal that exhibits intriguing bonding and oxidation states with profound geological implications. Our computational study identifies a crystalline phase of CaO 3 by reaction of CaO and O 2 at high pressure and high temperature conditions; ensuing experiments synthesize this rare compound under compression in a diamond anvil cell with laser heating. High-pressure x-ray diffraction data show that CaO 3 crystal forms at 35 GPa and persists down to 20 GPa on decompression. Analysis of charge states reveals a formal oxidation state of -2 for ozone anions in CaO 3 . These findings unravel the ozonide chemistry at high pressure and offer insights for elucidating prominent seismic anomalies and oxygen cycles in Earth’s interior. We further predict multiple reactions producing CaO 3 by geologically abundant mineral precursors at various depths in Earth’s mantle.

58 GEOSCIENCES↗

Clear Lake Volcanic Field Disasters: Creating a Deformation Record Using InSAR to Assess Hazards and Detect Volcanic Unrest in Clear Lake Volcanic Field

Clear Lake Volcanic Field (CLVF) in northern California is at a high threat potential for volcanic hazards. Eruptions leading to increased seismic activity could result in silicic domes, cinder cones, and flows that would be dangerous to the residential areas near the volcanic field. Remotely sensed Earth observations can reveal volcanic processes in the subsurface, which are essential to the timely monitoring of potential volcanic activity. In particular, Sentinel-1 C-band Synthetic Aperture Radar (C-SAR) and Digital Elevation Model (DEM) data capture relative surface deformation at unprecedented high spatial and temporal resolutions. Leveraging C-SAR and DEMs, we conducted interferometric analysis from January 2016 to December 2023. Our results demonstrate 1) the mean surface displacement velocity of the Clear Lake Volcanic Field is measured to undergo 5 to 10-centimeter scale deformation and shows a strong relationship with the surrounding faults, 2) apparent seasonal differences in rates of surface change, and 3) seismic activity associated with the geyser geothermal field has a strong association with cumulative surface displacement, with active fault zones having 2 to 5 cm of additional displacement. Results indicate that deformation is linked to deep pressure sources causing stresses on the surficial environment that should be considered in hazard mitigation. This study provides a baseline of historic deformation, aiding hazard analysts in communication efforts and streamlining decision-making for potential risks to region residents.

remote sensing↗

Creating a Deformation Time Series Utilizing InSAR to Assess Hazards and Detect Volcanic Unrest in Clear Lake Volcanic Field, California

According to the U.S Geological Survey, the Clear Lake Volcanic Field (CLVF) in northern California is at a high threat potential for volcanic hazards. Eruptions leading to increased seismic activity could result in silicic domes, cinder cones, and flows that would be dangerous to the residential areas near the volcanic field. Remotely sensed Earth observations can reveal volcanic processes in the subsurface, which are essential to the timely monitoring of potential volcanic activity. In particular, Sentinel-1 C-band Synthetic Aperture Radar (C-SAR) and Digital Elevation Model (DEM) data capture relative surface deformation at unprecedented high spatial and temporal resolutions. Leveraging C-SAR and DEMs, we conducted interferometric analysis from January 2016 to December 2023. Our results demonstrate 1) the mean surface displacement velocity of the Clear Lake Volcanic field is measured to undergo 5 to 10-centimeter scale deformation and shows a strong correlation with the surrounding faults, 2) apparent seasonal differences in rates of surface change, and 3) seismic activity associated with the geyser geothermal field has a strong correlation with cumulative surface displacement, with active fault zones having 2 to 5 cm of additional displacement. Results indicate that deformation is linked to deep pressure sources causing stresses on the surficial environment that should be considered in hazard mitigation. This study provides a baseline of historic deformation, aiding hazard analysts in communication efforts and streamlining decision-making for potential risks to region residents.

Ivan Tochimani-Hernandez↗

Preliminary analysis of the JAPE ground vehicle test data with an artificial neural network classifier

Remotely sensing and classifying military vehicles in a battlefield environment have been the source of much research over the past 20 years. The ability to know where threat vehicles are located is an obvious advantage to military personnel. In the past active methods of ground vehicle detection such as radar have been used, but with the advancement of technology to locate these active sensors, passive sensors are preferred. Passive sensors detect acoustic emissions, seismic movement, electromagnetic radiation, etc., produced by the target and use this information to describe it. Deriving the mathematical models to classify vehicles in this manner has been, and is, quite complex and not always reliable. However, with the resurgence of artificial neural network (ANN) research in the past few years, developing models for this work may be a thing of the past. Preliminary results from an ANN analysis to the tank signatures recorded at the Joint Acoustic Propagation Experiment (JAPE) at the US Army White Sands Missile Range, NM, in July 1991, are presented.

Larsen, Nathan F.↗

Utah FORGE: Seismic DAS and Geophone Borehole Data Processing and 3D Imaging of Vp/Vs Ratio in the 2024 Stimulated Reservoir

This dataset includes a final report and a 3D velocity model derived from seismic DAS and geophone borehole data collected during the April 2024 stimulation of the reservoir at Utah FORGE. The report details the processing of over 50,000 P- and S-wave travel times used in a tomographic inversion to estimate the Vp/Vs ratio, revealing anomalies adjacent to wells 16A and 16B that may be associated with injected fracturing fluids. Additionally, Wadati analysis of more than 27,000 pairs of differential P- and S-wave travel times supports this interpretation with a high Vp/Vs estimate of 1.86, averaged over the dimensions of the earthquake cluster. The accompanying data file provides a 3D model of P- and S-wave velocities and Vp/Vs ratios, structured with spatial coordinates and velocity values.

15 GEOTHERMAL ENERGY↗

Seismic Study of The Solar Interior: Inferences from SOI/MDI Observations during Solar Activity

The principal investigator describes several types of solar research conducted during the reporting period and gives a statement of work to be performed in the following year. Research conducted during the reporting period includes: exhaustive analysis of observational and instrumental effects that might cause systematic errors in the characterization of high-degree p-modes; study of the structure, asphericity and dynamics of the solar interior from p-mode frequencies and frequency splittings; characterizing the solar rotation; Time-Distance inversion; and developing and using a new peak-fitting method for very long MDI time series at low degrees.

Korzennik, Sylvain G.↗

Fractal applications to complex crustal problems

Complex scale-invariant problems obey fractal statistics. The basic definition of a fractal distribution is that the number of objects with a characteristic linear dimension greater than r satisfies the relation N = about r exp -D where D is the fractal dimension. Fragmentation often satisfies this relation. The distribution of earthquakes satisfies this relation. The classic relationship between the length of a rocky coast line and the step length can be derived from this relation. Power law relations for spectra can also be related to fractal dimensions. Topography and gravity are examples. Spectral techniques can be used to obtain maps of fractal dimension and roughness amplitude. These provide a quantitative measure of texture analysis. It is argued that the distribution of stress and strength in a complex crustal region, such as the Alps, is fractal. Based on this assumption, the observed frequency-magnitude relation for the seismicity in the region can be derived.

Turcotte, Donald L.↗

Structural Analysis of the US ITER ECH Transmission Line System

The electron cyclotron heating (ECH) and current drive system is one of the main plasma heating systems for ITER. It uses high-power microwave beams with the power deposition location steerable across the plasma cross section. Microwave power is conveyed via transmission lines (TLs) that run from the gyrotrons in the radio frequency building through the assembly hall and tokamak building to the ECH launchers within the tokamak vacuum vessel. The ECH system includes a vast array of interconnected TL waveguides, in-line components, and support structures.Finite element (FE) modeling provides an essential means of simulating the system, applying loads and determining deflections, rotations, forces, moments, and stresses in order to evaluate various structural and microwave transmission performance metrics. A representative FE model of the overall ECH TL system is developed in ANSYS®. This top-level model defines the centerline of the waveguide system. Waveguide segments are represented by line elements (beams and pipes) with equivalent section properties, and support structures are represented by boundary conditions. A systematic approach is used to model each ECH component with lumped masses and structurally equivalent stiffness matrices or ANSYS superelements.The top-level TL FE model is used to evaluate the various loads (thermal, vacuum, seismic, etc.) and operating scenarios. The top-level model directly calculates stresses in the straight aluminum waveguide segments. The model provides the forces and moments acting on the in-line components for detailed submodel assessments. Displacement results from the top-level analysis feed into a separate microwave performance model to help determine operational efficiency. All TL performance and thermal-structural requirements are met, as specified by the applicable codes and standards, and successfully documented in numerous technical reports and demonstrated at the final design review.

Shanmugasundaram, Aravind↗

Analysis of Data Fusion Between Waveform Events and Radionuclide Detections Reported in 2021 by the International Data Centre

During this project, we assessed the percentage of events reported by the International Data Centre (IDC) that are candidates for evaluation using data fusion methods. We calculated the percentage of waveform (seismic, hydroacoustic, infrasound, SHI) events that may be associated with radionuclide (RN) detections and vice versa. For this task, data fusion means identifying RN detections that may be associated with SHI events. The task of producing a fused SHI and RN bulletin is included in the IDC operational manual. As of 2022, the production of a fused event bulletin remains incomplete, largely because methods for producing the bulletin are under development. Nonetheless, simple methods for data fusion have been developed and preliminary software tools are available at the IDC to associate SHI events and RN detections. This report is a first look at the percentage of SHI events that may be candidates for fusion with RN detections or bulletins. One year of IDC bulletin data (2021) is considered in this report. It is important to note that there were no known nuclear explosions during this time period. Hence, we expect that events during this time period would generate SHI and RN signals. If many SHI events are candidates for fusion with RN data, greater priority should be placed on establishing screening procedures before and/or during the fusion process. By contrast, if the number of candidate events is diminishingly small, then data fusion is not a driver for further development of event or detection screening at the IDC. LLNL reviewed the events reported in the Standard Event Bulletin (SEB) and determined which events are candidates for fusion with the Reviewed Radionuclide Reports (RRR). Much of the work involves establishing a metric or metrics for determining whether a Reviewed Event Bulletin (REB) event is a candidate for data fusion.

58 GEOSCIENCES↗

Analysis of passive acoustic ranging of helicopters from the joint acoustic propagation experiment

For more than twenty years, personnel of the U.S.A.E. Waterways Experiment Station (WES) have been performing research dealing with the application of sensors for detection of military targets. The WES research has included the use of seismic, acoustic, magnetic, and other sensors to detect, track, and classify military ground targets. Most of the WES research has been oriented toward the employment of such sensors in a passive mode. Techniques for passive detection are of particular interest in the Army because of the advantages over active detection. Passive detection methods are not susceptible to interception, detection, jamming, or location of the source by the threat. A decided advantage for using acoustic and seismic sensors for detection in tactical situations is the non-line-of-sight capability; i.e., detection of low flying helicopters at long distances without visual contact. This study was conducted to analyze the passive acoustic ranging (PAR) concept using a more extensive data set from the Joint Acoustic Propagation Experiment (JAPE).

Carnes, Benny L.↗

Distributed Acoustic Sensing (DAS) on Opportunistic Networks_ A Feasibility Study Utilizing Fiber Optic Infrastructure at NEON Sites for DAS

This study demonstrates the feasibility of using Distributed Acoustic Sensing (DAS) on the opportunistic optical fiber infrastructure presented by NEON (National Ecological Observatory Network) sites. The test took place at the NEON Terrestrial Site designated D10 CPER. Unlit portions of the local fiber network were successfully spliced together and optimized for length and a Silixa iDAS sensor was installed. The resulting array (composed of 12 spliced fiber runs) totaled 875m in length. Over a period of 2 weeks, we gathered a high spatial and temporal resolution data set to characterize discernable signals and understand background noise. The data set was reviewed to extract wave signatures of identifiable origin, both anthropogenic and otherwise. Continent wide distribution of these sites makes them an interesting candidate for DAS and identifying measurable signal sources informs future DAS application on opportunistic fiber networks. These observations validate the utility of both the NEON and similarly small fiber networks (<1km) for geoscience, environmental, and security DAS applications and indicate possibilities for further inquiry in methods for DAS data collection and analysis.

42 ENGINEERING↗

Stress analysis using direct-S wavelets produced by a vertical vibrator

We have determined how to extract the azimuth of maximum horizontal stress (SHmax) in deep rocks by doing a simple 360° mathematical rotation of a downgoing direct-S wavelet generated at the baseplate of a surface-based vertical vibrator. We worked with direct-S wavelets that travel through stressed rocks to a deep horizontal vertical seismic profile (VSP) geophone. We find that the azimuth where a polarity reversal occurs in mathematical rotations of this downgoing direct-S wavelet defines the azimuth of the SHmax in the overlying rocks. We tested this direct-S wavelet rotation method for determining the SHmax azimuth at a site in the Illinois Basin using legacy VSP data acquired in 2013. SHmax azimuths indicated by this simple wavelet-rotation method were determined when vertical vibrators were stationed at zero offset, at far offset, and at different azimuths around a VSP receiver well. These VSP-based SHmax azimuths agreed with the azimuth of the SHmax found by traditional minifrac tests in the VSP receiver well. This simple VSP data analysis procedure for detecting the azimuth of maximum horizontal stress has never, to our knowledge, been reported or discussed in the geophysical literature. This technical finding should be of interest to the worldwide geophysical community, especially to people who need to monitor how stress fields shift when fluids are injected into, or extracted from, deep porous reservoirs.

Geochemistry & Geophysics↗

Geothermal Play-Fairway Analysis of Washington State Prospects: Final Report

The Washington State Geothermal Play-Fairway Analysis overcomes the exploration challenges posed by dense vegetation, glacial deposits, and extreme precipitation. The geothermal play-fairways we target are locations where heat, permeability, and saturated porosity are present in sufficient volume to provide adequate heat exchange at depths accessible by modern drilling technology. The three study areas lie along the Cascade Range magmatic arc and are near Mount Baker, Mount St. Helens, and the Wind River Valley. The seven-year project is divided into three phases. In Phase 1 we build on a previous statewide assessment of geothermal resources and develop an initial modeling approach. The results are a series of favorability, uncertainty, and risk maps for three targeted study areas. Based on these initial results, we collect new geologic and geophysical data to further refine our modeling and reduce exploration uncertainty in Phase 2. We improve the modeling method to handle the new data and update the favorability, uncertainty, and risk maps. We also update the conceptual geothermal resource models. In Phase 3 we validate our modeling approach by drilling two temperature-gradient holes and collecting and analyzing core, image logs, and new geochemistry. Our modeling approach improves on an earlier statewide method through a more-rigorous and detailed assessment of heat and permeability. Permeability potential is assessed through geomechanical modeling of the deformation that can generate and maintain reservoir porosity and permeability. Metrics to inform heat potential include temperature-gradient wells, which are sparse in Washington; proximity of Quaternary volcanic vents and young intrusive rock; spring temperature; and reservoir temperature inferred from geothermometry. We weight the individual components using an expert-guided approach known as the Analytical Hierarchy Process. During Phase 2 we also develop a fluid-filled fracture model, and an infrastructure model that helps to delineate areas which are more favorable for geothermal development based on proximity to transmission lines, elevation, land ownership and use restrictions, and availability of process water. New geologic and geophysical data is collected during Phase 2 in each of our three main study areas. At Mount Baker and north of Mount St. Helens we conduct 1:24,000-scale geologic mapping and lidar analysis to better constrain the location and character of surface faults; detailed mapping in the Wind River Valley was completed just prior to the start of this project. Ages of intrusive rocks are determined with 40 Ar/ 39 Ar geochronology, though all of our samples are Miocene or older. We collect ground based gravity observations (a total of 1,580 new stations) in all of our study areas and ground-based magnetic lines (a total of 93 km) at Mount Baker. These data are combined with existing gravity and aeromagnetic data and used to constrain fault locations and geometry. Two to three cross sections are constructed at each study area using the mapped surface geology and forward-modeling of the gravity and magnetic data; these cross sections form the basis for our updated conceptual models. We collect magnetotelluric surveys at Mount Baker and Mount St. Helens and these data are inverted to form a resistivity model from the surface to about 10 km depth; each model shows conductive zones that can be interpreted as upwelling geothermal fluids. At Mount St. Helens we deploy a passive seismic array and use the newly detected events to refine the location of the Saint Helens seismic zone. We also employ ambient-noise tomography to develop a detailed seismic-velocity model for the study area and use this model to help constrain our cross sections and conceptual model. Based on the new data collected during Phase 2—and our updated models—we develop a campaign of temperature-gradient holes and core analysis to validate our modeling in Phase 3. Drill hole MB76-31 is located near Little Park Creek, 11 km west-southwest of the summit of Mount Baker, and is 1,471 ft deep. About 410 ft of core from the lower portion of the hole—and image logs from ~175 ft below ground surface to the bottom—are collected and analyzed. Water samples are collected and processed for geothermometry. Drill hole MSH17-24 is located along upper Schultz Creek, 16 km north-northeast of Mount St. Helens and has core from 470 ft to the bottom at 1,053 ft. We did not collect image logs due to borehole stability concerns, but water samples are collected and analyzed for geothermometry. Repeat temperature-gradient measurements are made at both sites and thermal conductivity is measured from core samples. At MB76-31, the equilibrated temperature gradient of 64°C/km and calculated heat flow of 141–159 mW/m 2 is more than twice the regional average. Detailed mapping and analysis of the core, coupled with correlation to the image logs, indicates a history of permeability generation consistent with our predictions of high permeability. Because the site has high favorability in the Phase 2 model, we consider the results a positive validation of the modeling. At site MSH17-24, the equilibrated temperature gradient of ~15°C/km and calculated heat flow of 41–43 mW/m 2 are similar to regional. Geochemical analysis of the water samples indicates a meteoric source without any geothermal component. Detailed outcrop-based mapping of fault exposures near the drill site and analysis of image logs from nearby boreholes indicates a history of permeability generation consistent with our predictions. Because the site has low favorability in the Phase 2 model, we consider the results a positive validation of the modeling. Together, the two sites provide a reasonably positive validation of the Phase 2 modeling and should encourage future use of this modeling approach.

15 GEOTHERMAL ENERGY↗

On asymmetric dark matter constraints from the asteroseismology of a subgiant star

ABSTRACT The asteroseismic modelling of solar-like stars has proved to be valuable in constraining dark matter. In this work, we study for the first time the influence of asymmetric dark matter (ADM) in the evolution of a subgiant star (KIC 8228742) by direct comparison with observational data. Both spectroscopic and seismic data are analysed with a new approach to the stellar calibration method, in which DM properties can also be considered as free inputs. In another phase of this study, a calibrated standard stellar model (without DM) is used as the benchmark for DM models. We find that the latter models consistently outperform the former for 10−40 ≤ σSD < 10−38 cm2, hinting that the presence of ADM in stars of this type does not go against observations. Moreover, we show that stellar seismology allows us to suggest exclusion limits that complement the constraints set by direct detection experiments. Different seismic observables are proposed to study DM properties and ΔΠℓ is found to be the most reliable, having the potential to build future DM exclusion diagrams. This new methodology can be a powerful tool in the analysis of the data coming from the next generation of asteroseismic missions.

Rato, João (ORCID:0000000320802361)↗

Thermal equation of state for zoisite: Implications for the transportation of water into the upper mantle and the high-velocity anomaly in the Farallon plate

Seismic tomography studies have revealed a high-velocity anomaly at depths between 100 km and 300 km in the Farallon plate. However, the reasons for the high-velocity anomaly continue to be debated. An analysis of the mineral proportions of eclogites exhumed from the Farallon plate shows that the average amount of zoisite in eclogite is ∼16.0 vol%. Therefore, the presence of zoisite eclogite needs to be considered to explain the high-velocity anomaly of the Farallon plate. However, the thermal equation of state and stability of zoisite have not been fully investigated under high pressure−temperature (P-T) conditions. We investigated the high-pressure and high-temperature behavior of natural zoisite utilizing synchrotron single crystal-X-ray diffraction (XRD). The results indicate that zoisite is metastable up to 24.8 GPa and 700 K. We obtained the ambient unit-cell volume V0 = 901.26(3) Å3 by synchrotron single crystal-XRD measurement. We also fitted the pressure-volume-temperature data to a high-temperature Birch-Murnaghan equation of state and obtained the zero-pressure bulk modulus K0 = 134.7(8) GPa, the temperature derivative of the bulk modulus (∂K/∂T)P = −0.011(4) GPa/K, and the thermal expansion coefficients α0V = 1.9(7) × 10−5 K−1 and α1V = 3(2) × 10−8 K−2 when the pressure derivative of bulk modulus (K′0) is fixed at 4. By incorporating the results from previous studies, we calculated the density and bulk sound velocity profiles of zoisite eclogite along the Farallon plate geotherm. Finally, we infer that zoisite could carry water to depths of ∼300 km within cold subducting slabs and that the bulk sound velocity of typical zoisite eclogite with ∼61.0 vol% omphacite, ∼23.0 vol% garnet, and ∼16.0 vol% zoisite could cause the high-velocity anomalies at depths of 100−300 km in the Farallon plate.

Geology↗