Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “seismic sources”

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 253 records · Page 14
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Relocating microseismicity from downhole monitoring of the Decatur CCS site using a modified double-difference algorithm

SUMMARY The injection of CO2 at the Decatur carbon capture and storage site has generated significant microseismic activity, which occurs in distinct spatial clusters up to approximately 2.2 km from the primary injection well. Accurate and precise event locations are vital for the characterization of the microseismicity to help understand the reservoir response to the CO2 injection, whilst enabling the identification of minor faults and fractures below the resolution of conventional active seismic imaging. However, microseismic monitoring of fluid injection sites, such as Decatur, is often performed using a network of borehole sensors often from a single well. While these downhole sensors have excellent detection capabilities, their poor azimuthal coverage limits the ability to precisely determine event locations. We have developed a modified double-difference relocation algorithm suitable for both 1-D and 3-D velocity models, and which incorporates differential back azimuth observations to allow the benefits of the original double-difference algorithm to be applicable to a downhole microseismic monitoring setting. Applying the modified double-difference algorithm to the microseismicity at Decatur, we have successfully relocated 4293 events. The relocation included over 59 million observations for 757 285 event pairs, split across seven geographic regions. Despite the majority of observations being recorded in only two boreholes, with an almost identical azimuthal coverage, the results have shown to be reliable with significantly reduced residuals and low uncertainties associated with the final locations. We have analysed the residuals in terms of their association with each geographic region, data type, station and individual events, to fully appreciate their influence in the inversion and the fit of the data to the final set of event locations. For each region, the relocated seismicity has become less diffuse with improved clustering, and with newly visible linear features often orientated in a NE–SW direction. These results show the potential improvements that can be made to microseismic event locations recorded by a borehole network with a limited and variable azimuthal distribution.

58 GEOSCIENCES↗
DOE OSTI2022

Analysis of the Source Physics Experiment Dry Alluvium Geology Ball Drop Sources

Two balls were dropped approximately 150 m down an open borehole as part of an investigation into the aftershocks from a 51 ton TNT-equivalent chemical explosion. The ball drop sources, a concrete-filled basketball and a bowling ball, were intended to mimic falling debris within an underground explosion borehole and cavity. Falling debris or sloughing of the borehole and cavity walls is thought to be one of the sources of explosion aftershocks. The signals from the ball drop sources were recorded by the same seismic network that monitored the aftershocks of the chemical explosion, so the ball drop recordings were used as templates to search through the continuous stream of aftershocks and detect correlated signals as a proportion of all aftershocks. The search resulted in only four similar event recordings out of the thousands of detected aftershocks. A synthetic test showed that the small, high-frequency ball drop templates were limited to producing highly correlated detections for events occurring in a small region around the impact location at 150 m depth. The high-frequency templates were necessary since the small ball drop signals were only above the noise at frequencies greater than 10 Hz. Unfortunately, the impact location is where the borehole was filled with sand in preparation for the subsequent chemical explosion in the test series, making it 150 m above the explosion created cavity centered at 300 m depth. This puts the ball drop templates 150 m away from the cavity, where many more events from falling debris would be expected. Further synthetic tests confirmed events emanating from the cavity 150 m below the impact site would have correlated with the lower frequency template and been detected had the template been usable below 10 Hz. Therefore, we conclude that the four detected events were most probably falling debris within the open borehole that impacted the floor at 150 m depth just as the ball drop sources had, and that this type of explosion aftershock mechanism comprises less than 1% of all explosion aftershocks.

58 GEOSCIENCES↗
DOE OSTI2024

Precise relative magnitude measurement improves fracture characterization during hydraulic fracturing

SUMMARY Microseismic monitoring is an important technique to obtain detailed knowledge of in-situ fracture size and orientation during stimulation to maximize fluid flow throughout the rock volume and optimize production. Furthermore, considering that the frequency of earthquake magnitudes empirically follows a power law (i.e. Gutenberg–Richter), the accuracy of microseismic event magnitude distributions is potentially crucial for seismic risk management. In this study, we analyse microseismicity observed during four hydraulic fracture treatments of the legacy Cotton Valley experiment in 1997 at the Carthage gas field of East Texas, where fractures were activated at the base of the sand-shale Upper Cotton Valley formation. We perform waveform cross-correlation to detect similar event clusters, measure relative amplitude from aligned waveform pairs with a principal component analysis, then measure precise relative magnitudes. The new magnitudes significantly reduce the deviations between magnitude differences and relative amplitudes of event pairs. This subsequently reduces the magnitude differences between clusters located at different depths. Reduction in magnitude differences between clusters suggests that some attenuation-related biases could be effectively mitigated with relative magnitude measurements. The maximum likelihood method is applied to understand the magnitude frequency distributions and quantify the seismogenic index of the clusters. Statistical analyses with new magnitudes suggest that fractures that are more favourably oriented for shear failure have lower b-value and higher seismogenic index, suggesting higher potential for relatively larger earthquakes, rather than fractures subparallel to maximum horizontal principal stress orientation.

58 GEOSCIENCES↗
NASA NTRS1974

Lunar seismicity, structure, and tectonics

Natural seismic events have been detected by the long-period seismometers at Apollo stations 16, 14, 15, and 12 at annual rates of 3300, 1700, 800, and 700, respectively, with peak activity at 13- to 14-day intervals. The data are used to describe magnitudes, source characteristics, and periodic features of lunar seismicity. In a present model, the rigid lithosphere overlies an asthenosphere of reduced rigidity in which present-day partial melting is probable. Tidal deformation presumably leads to critical stress concentrations at the base of the lithosphere, where moonquakes are found to occur. The striking tidal periodicities in the pattern of moonquake occurrence and energy release suggest that tidal energy is the dominant source of energy released as moonquakes. Thus, tidal energy is dissipated by moonquakes in the lithosphere and probably by inelastic processes in the asthenosphere.

Lammlein, D. R.↗
DOE OSTI2025

Geometric Interpretation of the Cluster Location Problem Part I: Theory

We present a new framing of the seismic location problem using principles drawn from differential geometry. Our interpretation relies upon the common assumption that travel times observed across a network are continuous, differentiable functions of source location. In consequence, travel‐time functions constitute a differentiable map between the source region and a Riemannian manifold. The manifold is said to be the image of the source region embedded in a generally high‐dimension travel‐time vector space. A cluster of events in the source region has an image of discrete points on the manifold, that, except in the simplest cases, cannot be viewed directly. However, it is possible to project the image of a cluster into a tangent space of the manifold for direct visualization. The projection operator can be computed directly from the data without a velocity model, but produces a distorted rendering of the cluster geometry. With a model we can predict the distortions and correct them to estimate cluster geometry. We develop these points with the simplest possible example, one for which direct visualization of the manifold is possible, using the example as an introduction to the relevant concepts from differential geometry in a familiar setting. The tangent space, a local linearization of the manifold, plays a key role. We develop a metric to estimate the limits of linearization, that is, to determine when the curvature of the manifold invalidates the linear assumption. We also examine the interplay of model error, inadequate network geometry, and pick error. We then generalize our results from the simple case to the general case of 3D source regions observed by general networks. Although we do suggest a new “project and correct” method for location, we do not develop it into a practical algorithm. In conclusion, our intention rather is to highlight new analytical methods grounded in differential geometry.

East Pacific Ocean Islands↗
DOE OSTI2026

Evolution of storage monitoring – update in response to commercial and regulatory drivers

Carbon Capture and Storage (CCS) is in transition from first-of-a kind projects and research-orientated pilots to commercially-motivated applications. Monitoring results from many newly developed and planned large scale commercial projects are limited; however, it is worthwhile to assess their evolution and consider new strategies as part of an effort to assess and document best practices. Commercial monitoring is targeted to activities that comply with regulatory drivers and de-risk investments. Commercial monitoring also supports accounting that storage has occurred and is tied to project financing. It deals with long time frames and large volumes injected into multiple wells and multiple projects in favorable areas. We see developing trends toward reproducible workflows that systematically reduce risks and clarify expectations for oversight and long-term surveillance. Monitoring techniques showing increasing trends include injection zone pressure as a history-matching and compliance tool. To reduce cost and environmental impact of time-lapse seismic data collection, deploying new approaches and tools, such as use of fibre and installed sources are increasingly applied. Concern over the risk of induced seismicity by regulatory bodies and the general public has increased, which has also resulted in increased monitoring. Some techniques used in the early research phases have been sidelined or used only in restricted applications. For example, geochemical analyses in the injection zone as well as the environment are now being deployed less than it was in research-oriented programs, except in the US where it is required by the permitting process. Expectations of frequent area-wide near surface monitoring have also decreased.

25 ENERGY STORAGE↗
DOE OSTI2021

The First Detection of an Earthquake From a Balloon Using Its Acoustic Signature

Extreme temperature and pressure conditions on the surface of Venus present formidable technological challenges against performing ground-based seismology. Efficient coupling between the Venusian atmosphere and the solid planet theoretically allows the study of seismically generated acoustic waves using balloons in the upper atmosphere, where conditions are far more clement. However, earthquake detection from a balloon has never been demonstrated. We present the first detection of an earthquake from a balloon-borne microbarometer near Ridgecrest, CA in July 2019 and include a detailed analysis of the dependence of seismic infrasound, as measured from a balloon on earthquake source parameters, topography, and crustal and atmospheric structure. Our comprehensive analysis of seismo-acoustic phenomenology demonstrates that seismic activity is detectable from a high-altitude platform on Earth, and that Rayleigh wave-induced infrasound can be used to constrain subsurface velocities, paving the way for the detection and characterization of such signals on Venus.

58 GEOSCIENCES↗
NASA NTRS2001

Ar Atmosphere: Implications for Structure and Composition of Mercury's Crust

We examine the possibilities of sustaining an argon atmosphere by diffusion from the upper 10 km of crust, and alternatively by effusion from a molten or previously molten area at great depth . Ar-40 in the atmospheres of the planets is a measure of potassium abundance in the interiors since Ar-40 is a product of radiogenic decay of K-40 by electron capture with the subsequent emission of a 1.46 eV gamma-ray. Although the Ar-40 in the earth's atmosphere is expected to have accumulated since the late bombardment, Ar-40 in surface-bounded exospheres is eroded quickly by photoionization and electron impact ionization. Thus, the argon content in the exospheres of the Moon, Mercury and probably Europa is representative of current effusion rather than accumulation over the lifetime of the body. Argon content will be a function of K content, temperature, grain size distribution, connected pore volume and possible seismic activity. Although Mercury and the Moon differ in many details, we can train the solutions to diffusion equations to predict the average lunar atmosphere. Then these parameters can be varied for Hermean conditions. Assuming a lunar crustal potassium abundance of 300 ppm, the observed argon atmosphere requires equilibrium between the argon production in the upper 9 Km of the moon (1.135 x 10(exp -3) cm(exp -3) s(exp -1)) and its loss. Hodges et al. conclude that this loss rate and the observed time variability requires argon release through seismic activity, tapping a deep argon source. An important observation is that the extreme surface of the Moon is enhanced in argon rather than depleted, as one would expect from outgassing of radiogenic argon. Manka and Michel concluded that ion implantation explains the surface enhancement of Ar-40. About half of the argon ions produced in the lunar atmosphere would return to the surface, where they would become embedded in the rocks. Similarly, at Mercury we expect the surface rocks to be enhanced in Ar-40 wherever the magnetosphere has been open over time. Thus the measurement of surface composition will reveal the long-term effects of solar wind-magnetosphere interaction. Additional information is contained in the original extended abstract.

Killen, R. M.↗
NASA NTRS1976

Thermal movement of the regolith

High-frequency seismic events observed at the Apollo landing sites indicate that the regolith moves in response to diurnal variations of thermal stresses. While these events (thermal moonquakes) are observed by each of the Passive Seismic Experiment stations, the location of their sources was not possible until they were observed by the Lunar Surface Profiling Experiment (LSPE) array at the Apollo 17 site. Comparison of source locations with topographic features near the LSPE array indicate that thermal moonquakes are associated to some degree with large rocks and to a larger degree with craters. They are not associated with man-made disturbances of the lunar surface. It is suggested that thermal moonquakes are the seismic expression of a phenomenon that is actively degrading slopes on the lunar surface.

Duennebier, F.↗
NASA NTRS2003

Search for Far-Side Deep Moonquakes

A truly unexpected finding of the Apollo missions, 1969-1972, was a discovery of deep moonquakes. Analysis of the data from the seismic network, which operated for eight years from 1969 through 1977, identified more than 100 discrete source regions at depths approximately half way to the center of the moon. Their distribution, however, was not uniform, as all but one of the source regions found were on the front hemisphere of the moon. Thus, a question remains whether the observed one-sided distribution of deep moonquake sources represents their true distribution or instead occurs because all seismic stations are on the near side of the moon. If it is the former, it means that the interior of the moon is truly asymmetric, structurally and dynamically; if it is the latter, it means that we simply did not identify most moonquakes on the far side and that their identification will be a great help in investigating the deep interior of the moon, including existence of a possible metallic core.

Nakamura, Y.↗
DOE OSTI2020

Play Fairway Analysis: Structurally Controlled Geothermal Systems in the Eastern Great Basin Extensional Regime, Utah

A research team with membership from the University of Utah/Energy & Geoscience Institute, the University of Utah/Dept. of Geology & Geophysics, and the Utah Geological Survey, undertook a play fairway analysis (PFA) for geothermal resources in the Eastern Great Basin (EGB) extensional tectonic regime of western Utah. This is a high-priority region for geothermal exploration because active Basin and Range (B&R) extension with volcanism having a N-S strike is superimposed upon pre-existing E-W belts of plutonic rocks and large-scale structural lineaments. Cumulative heat flow along the N-S strike of the state totals approximately 5 GWt above background stable interior. Three electricity producing power plants currently exist with substantial potential for increase. Succinctly, our PFA approach aims to resolve potential sources of heat and permeability in the region, which are the two principal criteria for establishing a geothermal resource. An initial Phase 1 was carried out using only existing geoscientific data in the area. Criteria selected for focusing heat potential include direct heat flow measurements in boreholes, magnetotelluric (MT) low resistivity anomalies, fluid/gas geochemistry, and proximity to recent volcanic eruptions. Permeability is established through geological structures (fault density, critically stressed areas, seismicity, gravity), and MT low resistivity anomalies. In Phase II of this PFA project, additional geological, geophysical and geochemical data were acquired and analysis carried out primarily over promising composite common risk (CCR) areas initially identified in Phase I in order to focus prospectivity and prepare for drilling recommendations. These prospects are near the Twin Peaks rhyolite field, high heat flow areas north of the producing Cove Fort system, and geophysical structure beneath the Crater Knoll area off the northeast flank of the Mineral Mountains. Additional data included MT site fill-in, structural mapping and analysis using high-resolution imagery, gravity and on-ground mapping, Nodal 3C passive seismic collection, and passive 3He surveying. Heat source and permeability potential are again expressed in terms of their individual common risk segment (CRS) maps, with a color scheme using green for most favorable (low risk) and red for least favorable (high risk). Diverse data types are united through the technique of probability kriging, which establishes prospectivity thresholds for each data type and then computes probability of exceeding that threshold over the PFA area. Modified CRS and CCRS maps are compared to those of Phase I to highlight tighter prospectivity focus. In doing so, the prospectivity threshold for heat was increased significantly to narrow the targeting. In the final Phase III of this project where a recommended deep thermal gradient hole was sited, additional geophysical, geochemical and geological field collection and analysis was carried out to refine drill hole targeting. This includes prospect-scale MT, gravity, structure, passive seismic deployment (Cove Fort area), and detailed 3He isotope profiling. It was the recommendation of the DOE Technical Monitoring Team (TMT) that one or more holes be sited in the north-ern Cove Fort area where legacy TG gradient holes showed high cumulative heat flow. These were to be of moderate depth, 2000-3000 feet, to reach the geothermal fluid table expected to start in excess of 1000 feet depth. The drilling organization stipulated by the DOE/GTO was that of the USGS Research Drilling Program (RDP) centered in Las Vegas, NV. A detailed well plan, appended to this report, was developed principally by Dr. Ben Barker consulting to University of Utah, Dr. Steve Pye on the DOE TMT, Mr. Steven Crawford of the USGS-RDP, and the project PI Phil Wannamaker. However, temperature and possible H 2 S at the systems lead to cancellation of the drilling last-minute as this appeared outside the experience base of the USGS-RDP. We hope to have the opportunity to revisit the test drilling and expand the Play Fairway Analysis of this region at some point in the future.

15 GEOTHERMAL ENERGY↗
DOE OSTI2021

Wastewater Disposal Has Not Significantly Altered the Regional Stress State in Southern Kansas

Wastewater disposal is primarily responsible for the increased seismicity rate since ~2013 in southern Kansas. Previous work that used shear-wave splitting (SWS) in southern Kansas interpreted an ~90° temporal rotation in the fast polarization direction and attributed it to increased pore pressures resulting from fluid injection. However, this interpreted rotation coincided with a change in the stations used to make the SWS measurements. We investigate the temporal variability of fast azimuths in southern Kansas by making SWS measurements on earthquake families with similar source–receiver paths recorded on a stable local seismic network. We select high-quality SWS measurements by investigating the stability of results across 65 different frequency bands between 0.5 and 15 Hz. We find that the fast polarization direction in southern Kansas is relatively constant with an average east-northeast (~N79°E) orientation between 2014 and 2017. Our fast polarization measurements are primarily a reflection of the maximum principal horizontal stress direction (SHmax). We observe a slight spatial change in SHmax to the northeast (~N55°E) near the Nemaha ridge in Oklahoma. However, we do not observe any significant temporal rotation of SHmax or variation in delay time (i.e., crack density) in southern Kansas, contrary to the earlier study. The previously interpreted ~90° rotation may either be a reflection of a very local stress change or a misinterpretation of SWS results potentially due to the use of inconsistent source–receiver paths. Our SWS measurements cover the period of peak wastewater disposal and seismicity rates and suggest an absence of significant temporal rotations in the local anisotropy and stress orientations associated with wastewater disposal.

Geochemistry & Geophysics↗
DOE OSTI2025

LATTE: open-source, high-performance traveltime computation, tomography and source location in acoustic and elastic media

Traveltime-based tomography and source location are fundamental approaches for imaging subsurface structures and understanding the spatiotemporal distribution of seismicity from local to global scales. We present an open-source, high-performance framework integrating eikonal equation solvers and adjoint-state theory for traveltime computation, velocity tomography, source location and joint tomography-location in 2-D/3-D acoustic and elastic media. We introduce novel regularization schemes based on total generalized p-variation, structural similarity and multitask machine learning to enhance the fidelity and interpretability of inverted models and source locations. Key features of our implementation also include the ability to leverage both absolute-difference and double-difference traveltime misfits for high-fidelity velocity tomography and source parameter estimation; support for traveltime computation and inversion in diverse 2-D/3-D scenarios with arbitrary source and receiver distributions; and a perturbation-based optimal step-size estimation method to reduce computational costs. In addition, our implementation employs shared-memory and distributed-memory parallelization to provide an efficient solution for traveltime computation, tomography, and source location. In conclusion, we validate the efficacy and accuracy of our approach through multiple synthetic data examples.

58 GEOSCIENCES↗
DOE OSTI2022

Report on the LLNL Global Seismic Waveform Tomography Modeling Project

Earth models have important applications for seismic event monitoring including the location and characterization of potential underground nuclear tests. Seismic waves generated by earthquakes and man-made events are altered by variations in Earth’s subsurface properties; and these distortions effect the arrival times of energy packets (phases) and obscure the characteristics of the original energy source (event mechanism). Three-dimensional models of Earth’s seismic properties can help explain and predict the distortions to the seismic wavefield and reveal the properties of a source that generated the waves, including the event location using the model-predicted timing of the phases and source mechanism using the modelpredicted waveform characteristics.

58 GEOSCIENCES↗
DOE OSTI2024

AmpSuite

Seismic amplitudes offer vital information about explosion source characteristics, including discrimination and yield estimation. To take advantage of this, we developed an interactive Python package to measure, control data quality, generate broad area propagation models and perform discrimination and estimate yield. Propagation models are essential in support of transportable yield and broad area discrimination. The key benefit of this package will be its ability to continuously integrate data and new techniques. The AmpSuite framework will provide standardized, repeatable, and accurate model generation and characterization routines. The capability is crucial for monitoring agencies tasked with rapid and high-quality seismic event characterization. The AmpSuite software includes a series of independent modules to perform: • Direct Phase Amplitude Measurement and Storage • Coda Envelope Measurement and Storage • Data Quality Control • New Propagation Model Developments • Seismic Discrimination and Analysis • Yield Estimation and supporting utility software. The AmpSuite software provides comprehensive solutions for monitoring agencies seeking to optimize model generation and event analysis within a contemporary Python framework. Stakeholders (AFTAC) have begun to move towards the Python language for scientific analysis as a new workforce emerges.

Alfaro, Richard↗
DOE OSTI2025

CLEAP Project: OR-SAGE Analysis for MT, UT, and CO States

The OR-SAGE tool is designed to use industry-accepted practices in screening sites and then employ the proper array of data sources through the considerable computational capabilities of GIS technology available at ORNL. The tool was developed to screen the potential for NPP siting on a national and regional basis. However, because of the tool granularity, it is often focused specifically on the immediate area around user sites of interest. If data center siting parameters can be added to OR-SAGE, the ability to evaluate data center siting on a localized scale will be beneficial.1 More than 60 data sets have been collected and processed by ORNL to develop exclusionary, avoidance, and suitability criteria for screening sites for a variety of power generation types, including nuclear power plants. Available site evaluation parameters include population density, slope, seismic activity, proximity to cooling-water sources, proximity to hazard facilities, avoidance of protected lands and floodplains, susceptibility to landslide hazards, and many others. All siting parameters should be considered as flags to inform siting decisions and should not be used to rule in or rule out any NPP site. Once data center siting parameters are identified, appropriate data sets will be collected and processed. The OR-SAGE process is very versatile. Essentially, OR-SAGE is a visual, relational database. The database partitions the contiguous United States, a total of 720 million hectares (~1.8 billion acres), into 100-m by 100-m (1 hectare or ~2.5 acre) cells. The database is tracking just under 700 million individual land cells. Successive suitability criterion is applied to each cell in the database. User-specified thresholds can be applied to each siting parameter data layer. In this manner, a variety of scenarios can be quickly and thoroughly evaluated. Data can be added and/or revised within OR-SAGE to address user interests. Siting security assessment capability is currently being added to OR-SAGE. Security is expected to be of concern at data centers whether it is collocated with a nuclear power generating technology or not. If data center is collocated with a nuclear power generating source, the security threat attractiveness level of both will likely increase. It will be of additional benefit if a potential data center site is also assessed for security vulnerability.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2026

Comprehensive earthquake catalogue update and spatiotemporal distribution analysis for Iraq and surrounding regions, northeastern Arabian Plate

The updated earthquake catalogue for Iraq covers the period from 1900 to the end of 2021 and includes over 37 000 recorded earthquakes. To create this comprehensive catalogue, five key steps were taken: compiling bulletins, calculating moment magnitudes, harmonizing magnitudes, establishing empirical conversion relations and evaluating the completeness of the catalogue. A notable enhancement in this update is the direct calculation of moment magnitudes for approximately 2800 earthquakes, achieved through the coda envelope technique and waveform data from the Mesopotamian Seismological Network (MPSN) in Iraq. This updated catalogue serves as a valuable resource for examining the spatiotemporal distribution of earthquakes, with respect to earthquake density, maximum moment magnitude and seismogenic depths. Additionally, the Gutenberg–Richter relationship was applied to calculate the a- and b-values specific to Iraq. The findings show that the Zagros Fold-Thrust Belt has a seismogenic layer (source) that ranges from 2 to 33 km deep and experiences high seismic activity. In contrast, the Mesopotamian Foredeep has a seismogenic layer ranging from 1 to 25 km deep and has lower seismic activity. The greatest seismic activity is concentrated around the Mandili-Badra-Teeb fault, which has experienced significant ruptures over time. The Outer Arabian Platform is identified as the main area of seismic activity, while additional activity occurs on the Inner Arabian Platform. Three major tectonic boundaries define the distribution of earthquakes in the northeastern Arabian Plate. These boundaries are defined by the Main Zagros Reverse Fault, the Zagros Foredeep Fault and the Anah Graben and Abu Jir-Euphrates Fault Zone. These boundaries highlight variations in seismicity levels and the spatial distribution of deformation in the region. The updated earthquake catalogue presented in this study is expected to play a vital role in regional seismicity assessments and seismic hazard analyses for Iraq and its surrounding areas.

58 GEOSCIENCES↗
NASA NTRS1979

The Moho as a magnetic boundary

Magnetism in the crust and the upper mantle and magnetic results indicating that the seismic Moho is a magnetic boundary are considered. Mantle derived rocks - peridotites from St. Pauls rocks, dunite xenoliths from the Kaupulehu flow, and peridotite, dunite, and eclogite xenoliths from Roberts Victor and San Carlos diatremes - are weakly magnetic with saturation magnetization values from 0.013 emu/gm to less than 0.001 emu/gm which is equivalent to 0.01 to 0.001 wt% Fe304. Literature on the minerals in mantle xenoliths shows that metals and primary Fe304 are absent, and that complex Cr, Mg, Al, and Fe spinels are dominant. These spinels are non-magnetic at mantle temperatures, and the crust/mantle boundary can be specified as a magnetic mineralogy discontinuity. The new magnetic results indicate that the seismic Moho is a magnetic boundary, the source of magnetization is in the crust, and the maximum Curie isotherm depends on magnetic mineralogy and is located at depths which vary with the regional geothermal gradient.

Wasilewski, P. J.↗