SEARCH · Search NASA

Results for “crustal structure”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 55 records · Page 3

All sources NASA NTRS DOE OSTI

DOE OSTI2021

Variability of near-fault seismic risk to reinforced concrete buildings based on high-resolution physics-based ground motion simulations

Broadband physics-based simulated earthquake ground motions are utilized to characterize the regional-scale seismic risk to modern reinforced concrete (RC) structures. A highly dense dataset of ground motions covering a 100-km urn:x-wiley:00988847:media:eqe3413:eqe3413-math-0001 40-km domain was generated using kinematic fault rupture models with varying rupture characteristics to represent shallow crustal earthquakes and resolved up to frequencies of 5 Hz. Over 40,000 nonlinear response history simulations of short- and mid-rise RC special moment frame buildings were conducted using simulation models that are capable of representing nonlinear behavior and component deterioration effects. The spatial variability of structural risk within a single earthquake scenario and between different rupture scenarios is examined, and the regions of strongest directivity effects and highest structural demands are identified. The structural demands may vary by factor of up to 8.0 at very short distances from the fault, and the large dispersion in the demands decreases significantly beyond a distance of 15 km. The interstory drift and member rotation demands are substantially impacted by important features of the geological structure and the characteristics of the rupture scenarios, particularly the presence of localized high-slip regions. The frequency characteristics of the structures are found to play an important role in determining the effects of near-fault ground motions on the structural response and expected damage. The results of this study suggest that the simulated ground motions, particularly those generated using the hybrid rupture approach, may offer reasonable structural risk estimates for low-frequency structures and conservative estimates for high-frequency structures.

58 GEOSCIENCES↗

DOE OSTI2022

Aerosol iron speciation and seasonal variation of iron oxidation state over the western Antarctic Peninsula

The iron (Fe) speciation and oxidation state have been considered critical factors affecting Fe solubility in the atmosphere and bioavailability in the surface ocean. In this study, elemental composition and Fe speciation in aerosol samples collected at the Palmer Station in the West Antarctic Peninsula were determined using synchrotron-based X-ray fluorescence (XRF) and X-ray Absorption Near-Edge Structure (XANES) spectroscopy. The elemental composition of coarse-mode (>1 mu m) Fe-containing particles suggests that the region's crustal emission is the primary source of aerosol Fe. The Fe minerals in these aerosol particles were predominantly hematite and biotite, but minor fractions of pyrite and ilmenite were observed as well. The Fe oxidation state showed an evident seasonal variation. Further, the Fe(II) content accounted for 71% of the total Fe in the austral summer, while this fraction dropped to 60% in the austral winter. Multivariate linear models involving meteorological parameters suggested that the wind speed, relative humidity, and solar irradiance were the factors that significantly controlled the percentage of Fe(II) in the austral summer. On the contrary, no relationship was found between these factors and the Fe(II) percentage in the austral winter, suggesting that atmospheric photoreduction and regional dust emission were limited. Moreover, the snow depth was significantly (p < 0.05) correlated with the aerosol Fe concentration, confirming the limiting effect of snow/ice cover on the regional dust emission. Given that the Antarctic Peninsula has experienced rapid warming during recent decades, the ice-free areas in the Antarctic Peninsula may act as potential dust sources.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2023

Does Magnetic Reconnection Occur in the Near Lunar Surface Environment?

The near lunar surface contains small-scale magnetic field structures that provide a natural test bed for observing plasmas with a non-zero Hall electric field, as well as potentially facilitating electron-only reconnection. This study presents observational evidence of magnetized electrons as well as demagnetized ions when THEMIS-ARTEMIS probe B reached an altitude of ~15 km above the lunar surface. Additionally, observations suggest the presence of a field line topology change and traversal of a closed magnetic field structure containing solar wind electrons, suggestive of magnetic reconnection having occurred at some point between the solar wind interplanetary magnetic field and a lunar crustal magnetic field. Thus, the observations presented here are consistent with previous studies that predict prominent Hall electric fields near lunar crustal magnetic fields and further suggest that the solar wind interplanetary magnetic field may reconnect with lunar crustal magnetic fields, most likely via electron-only reconnection.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2021

Thermal conductivity of CaSiO 3 perovskite at lower mantle conditions

Thermal conductivity (κ) of mantle minerals is key to understanding dynamics in the deep Earth. It controls the style of mantle convection and the timescale of cooling both the mantle and the core. Cubic CaSiO 3 perovskite (CaPv) is the third most abundant mineral in the lower mantle (7 vol %). Despite its importance, no theoretical or experimental estimate of CaPv's κ is available. Theoretical investigations of its properties are challenging because of its strong anharmonicity. Experimental measurements at relevant pressures and temperatures are equally challenging. Here we present ab initio results for CaPv's κ obtained using the phonon quasiparticle approach to address its strong anharmonicity. We also offer experimental measurements of κ up to 67 GPa and 1950 K. Predictions and measurements are in good agreement and reveal a surprisingly large κ for cubic CaPv that can be explained on the basis of its high crystal structure symmetry. Despite its relatively low abundance, CaPv's κ increases the lower mantle κ by ~10 %, if accounted for. κ of mantle regions enriched in subducted crustal materials will be more strongly impacted.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

DOE OSTI2023

Linking impact melt redox with crustal weathering regime

Abstract Impact melting and outgassing of impact melts were important processes shaping the crust and atmosphere of early Earth and Mars. The redox state of impact melts influences the gases they release, but the controls on impact melt redox are not well understood. Here, we present Fe micro–X-ray absorption near edge structure (μ-XANES) measurements of a suite of impact melt glasses from Lonar crater, India. As an impact into continental flood basalts, Lonar crater is an ideal analog for impacts into basaltic crust on early Earth and Mars. The Fe-μ-XANES technique permits robust characterization of the Fe3+/ΣFe ratio and therefore redox conditions in Lonar glasses. We found a range in Fe3+/ΣFe of 0.21–0.49 among Lonar glasses, including substantial intraclast variations. We conclude that a primary factor driving these variations was pre-impact oxidative weathering of target rocks. The coupling between crustal weathering regime and speciation of gases released from impact melts implies that impact melt outgassing could act as a feedback mechanism to reinforce prevailing surface redox conditions.

Geology↗

DOE OSTI2022

Glassy quantum nuclear pasta in neutron star crusts

We conduct a comprehensive survey of the shape parameter space of the nuclear pasta phases in neutron star crusts by conducting three-dimensional Hartree-Fock + BCS calculations. Spaghetti, waffles, lasagna, bicontinuous phases and cylindrical holes occupy local minima in the resulting constant-pressure Gibbs energy surfaces, implying multiple geometries coexist at a given depth. Notably, the bicontinuous phase, in which both the neutron gas and nuclear matter extend continuously in all dimensions appears over a large depth range. Our results support the idea that nuclear pasta is a glassy system. At a characteristic temperature, of order 10 8 – 10 9 K, different phases may become frozen into domains whose sizes we estimate to be 1–50 times the lattice spacing and over which the local density and electron fraction can vary. Above this temperature, very little long-range order exists and matter is an amorphous solid. Electron scattering off domain boundaries may contribute to the disorder resistivity of the pasta phases. Annealing of the domains may occur during cooling; repopulating of local minima during crustal heating might lead to temperature-dependent transport properties in the deep crust layers. We identify four regions distinguished by whether pasta is the true ground state, and whether the pasta structure allows delocalization of protons. The whole pasta region can occupy up to 70% of the crust by mass and 25% by thickness, and the layer in which protons are delocalized could occupy 45% of the crust mass and 15% of its thickness.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2024

Laboratory Study of Magnetic Reconnection in Lunar-relevant Mini-magnetospheres

Abstract Mini-magnetospheres are small ion-scale structures that are well suited to studying kinetic-scale physics of collisionless space plasmas. Such ion-scale magnetospheres can be found on local regions of the Moon, associated with the lunar crustal magnetic field. In this paper, we report on the laboratory experimental study of magnetic reconnection in laser-driven, lunar-like ion-scale magnetospheres on the Large Plasma Device at the University of California, Los Angeles. In the experiment, a high-repetition rate (1 Hz), nanosecond laser is used to drive a fast-moving, collisionless plasma that expands into the field generated by a pulsed magnetic dipole embedded into a background plasma and magnetic field. The high-repetition rate enables the acquisition of time-resolved volumetric data of the magnetic and electric fields to characterize magnetic reconnection and calculate the reconnection rate. We notably observe the formation of Hall fields associated with reconnection. Particle-in-cell simulations reproducing the experimental results were performed to study the microphysics of the interaction. By analyzing the generalized Ohm’s law terms, we find that the electron-only reconnection is driven by kinetic effects through the electron pressure anisotropy. These results are compared to recent satellite measurements that found evidence of magnetic reconnection near the lunar surface.

Astronomy & Astrophysics↗

DOE OSTI2021

Profiling the Quito basin (Ecuador) using seismic ambient noise

SUMMARY Quito, the capital of Ecuador, with more than 2.5 M inhabitants, is exposed to a high seismic hazard due to its proximity to the Pacific subduction zone and active crustal faults, both capable of generating significant earthquakes. Furthermore, the city is located in an intermontane piggy-back basin prone to seismic wave amplification. To understand the basin’s seismic response and characterize its geological structure, 20 broad and medium frequency band seismic stations were deployed in Quito’s urban area between May 2016 and July 2018 that continuously recorded ambient seismic noise. We first compute horizontal-to-vertical spectral ratios to determine the resonant frequency distribution in the entire basin. Secondly, we cross-correlate seismic stations operating simultaneously to retrieve interstations surface-wave Green’s functions in the frequency range of 0.1–2 Hz. We find that Love waves travelling in the basin’s longitudinal direction (NNE–SSW) show much clearer correlograms than those from Rayleigh waves. We then compute Love wave phase-velocity dispersion curves and invert them in conjunction with the HVSR curves to obtain shear-wave velocity profiles throughout the city. The inversions highlight a clear difference in the basin’s structure between its northern and southern parts. In the centre and northern areas, the estimated basin depth and mean shear-wave velocity are about 200 m and 1800 ms−1, respectively, showing resonance frequency values between 0.6 and 0.7 Hz. On the contrary, the basement’s depth and shear-wave velocity in the southern part are about 900 m and 2500 ms−1, having a low resonance frequency value of around 0.3 Hz. This difference in structure between the centre-north and the south of the basin explains the spatial distribution of low-frequency seismic amplifications observed during the Mw 7.8 Pedernales earthquake in April 2016 in Quito.

Pacheco, D. (ORCID:0000000282355096)↗

DOE OSTI2022

Atomic-scale mechanism of carbon nucleation from a deep crustal fluid by replica exchange reactive molecular dynamics simulation

Here we present a mechanistic model of carbon nucleation and growth from a fluid at elevated temperature (T) and pressure conditions, typical of those found in the shallow Earth’s lithosphere. Our model uses a replica exchange reactive molecular dynamics framework in which molecular configurations are swapped between adjacent T replica at regular intervals according to underlying statistical mechanics. This framework allows predicting complex molecular structures and thermodynamics while remaining computationally efficient. Here we simulate the reactivity of an unstable mixture of CO 2 and CH 4 at 1000 K and 1 GPa. We find that the path to thermodynamic equilibrium is initially entropy-driven, producing a diversity of short-lived species, including various alcohols with intermediate carbon oxidation states. Cyclic and polycyclic radicals that are sometimes resonance-stabilized form next and set the stage for carbon nucleation. The carbon exsolution process releases abundant water, is exothermic and starts with the nucleation of a large aggregate of hydrogenated graphene flakes from covalently bonded polycyclic units. The carbon backbone of this nucleus subsequently grows into a hydrogen-depleted fullerene-like structure, before evolving toward a partially bilayered graphene layer. Overall, our results show that the mechanism of graphitic C formation is certainly not bimolecular, and that it may involve a combination of key condensation and radical chain reactions. This will help understand the isotopic, and reactive characteristics of carbon-bearing fluids during their upward transit through the Earth’s mantle and crust. Moreover, the mechanistic insights outlined here present intriguing similarities with the process of soot and interstellar dust formation, which suggests that the widespread distribution of abiotic polyaromatic and graphitic material on Earth and beyond may reflect the prevalence of a fundamental chemical pathway.

58 GEOSCIENCES↗

DOE OSTI2021

Dynamics of Laterally Propagating Flames in X-Ray Bursts. II. Realistic Burning and Rotation

In this work, we continue to investigate two-dimensional laterally propagating flames in type I X-ray bursts using fully compressible hydrodynamics simulations. In the current study we relax previous approximations where we artificially boosted the flames. We now use more physically realistic reaction rates, thermal conductivities, and rotation rates, exploring the effects of neutron star rotation rate and thermal structure on the flame. We find that at lower rotation rates the flame becomes harder to ignite, whereas at higher rotation rates the nuclear burning is enhanced by increased confinement from the Coriolis force and the flame propagates steadily. At higher crustal temperatures, the flame moves more quickly and accelerates as it propagates through the atmosphere. If the temperature is too high, instead of a flame propagating across the surface the entire atmosphere burns uniformly. Our findings could have implications for the relationship between observed burst rise times and neutron star rotation and accretion rates. All of the software used for these simulations is freely available.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2021

Interplay Between Crustal‐Scale Thrusting, High Metamorphic Heating Rates, and the Development of Inverted Thermal‐Metamorphic Gradients: Numerical Models and Examples From the Caledonides of Northern Scotland

Abstract Understanding the processes that produce high prograde metamorphic heating rates and the development of inverted metamorphic sequences in collisional thrust belts remains a fundamental challenge for tectonics and metamorphic petrology. New 2D finite element models of crustal‐scale thrusts with variable slip rates (10, 20, 35, 50 km Myr −1 ) are used to examine how thrust sheet emplacement contributes to these processes. In the models, average prograde heating rates of 31–118 °C Myr −1 are observed in the footwall, with maximum transient heating rates of ∼167 °C Myr −1 occurring at the highest slip rate. Also, thrust sheet emplacement produces an inverted thermal gradient in the model footwall. At slip rates of 10 km Myr −1 , heating magnitudes >200 °C are observed >7 km structurally beneath the thrust plane. At slip rates of 20–50 km Myr −1 , models produce thermal penetration depths of 4–5 km for similar heating magnitudes. Petrologically determined heating rates for thrust footwalls in the Scandian orogenic wedge of northern Scotland mostly yield rates consistent with model results (10–230 °C Myr −1 ). The model‐derived magnitude of inverted thermal gradients is also similar to those indicated in texturally determined deformation temperature transects across the Scandian orogenic wedge (100 °C–180 °C). Combined, these results indicate that high heating rates can be produced by fault slip at typical plate velocities. Additionally, this implies that crustal‐scale thrusts are likely to produce inverted metamorphic sequences in most systems, provided that P‐T conditions during slip allow for metamorphic or recrystallization processes that would preserve evidence for such features.

Thigpen, J. Ryan↗

DOE OSTI2023

High-resolution lithospheric structure of continental China from joint inversion of surface wave and gravity data

The lithospheric structure of continental China has been previously determined by seismic travel time tomography, surface wave tomography, and joint inversion of body wave and surface wave data. However, due to the inherent limitations of seismic data, the lithospheric structure of continental China is still not well resolved in the shallow part and in some regions where the station coverage is relatively sparse. In this study, we aim at improving the lithospheric structure by joint inversion of seismic surface wave data and satellite gravity data to take advantage of the uniform distribution and complementary strength of the gravity data. Here, the empirical relationship between velocity and density is used as a bridge for joint inversion of surface wave and gravity data. The joint inversion shear-wave velocity Vs (density) model, named as USTClitho1.0g, can fit both surface wave and gravity data well. This high-resolution Vs model can better fit the active airgun source seismic arrival times and better delineate some features in continental China, such as the velocity contrast across the north–south gravity lineament (NSGL), the lithosphere thinning in eastern China, the crustal footprint of the Hainan mantle plume, the likely magma chamber beneath volcanos in northeast China, the middle-lower crust low velocity layer beneath the Tibetan plateau, and the tearing of subducted Indian mantle lithosphere. Our joint inversion Vs model can provide a reference model for geosciences in continental China and surrounding areas.

58 GEOSCIENCES↗

DOE OSTI2022

Seismic Tomographic Modeling of the Crust and Upper Mantle beneath Israel and the Middle East: Improved Resolution through Optimized Model Parameterization

Accurate regional seismic travel-time (RSTT) predictions rely on regional phases (e.g., Pg, Lg, Pn, Sn ) to account for 3D effects in the crust and upper mantle that are not captured by 1D models traditionally used for real-time location. The RSTT prediction model accounts for regional-scale crust and upper mantle structure globally by incorporating regional seismic phases into its travel-time calculations. Previous versions of the RSTT model have used a constant grid cell size of 1°. To improve the tomographic accuracy of recovering velocity structure at regional scales, we, in this study, perform data-driven grid refinement on the RSTT model down to a 0.125° grid (~14 km) in pursuit of two main goals: (1) to test the limits of RSTT capability and accuracy of determined velocity structure through variable grid refinement and (2) to image smaller structures in Israel and the Middle East and illuminate upper mantle dynamics operating in this complex tectonic area. We investigate the effects of model parameterization as grid cell size decreases and the trade-offs between recovered velocity structures. Our final dataset includes 4751 events and 499 stations that recorded 79,344 Pn and 7489 Pg . The variable grid refinement method allows recovery of finer-scale velocity structures and reduces travel-time residuals in areas with the highest data coverage. At smaller grid cell sizes, longer paths need to be upweighted to stabilize the inversion. Results illuminate tectonic features undefined in coarser grid-size models; in particular, we observe mantle perturbations related to the subduction zone around the Cyprian arc and crustal anomalies near the Dead Sea fault and throughout the Anatolian plate.

58 GEOSCIENCES↗

DOE OSTI2023

Internal Structure of the Central Garlock Fault Zone From Ridgecrest Aftershocks Recorded by Dense Linear Seismic Arrays

Abstract We provide high‐resolution seismic imaging of the central Garlock fault using data recorded by two dense seismic arrays that cross the Ridgecrest rupture zone (B4) and the Garlock fault (A5). Analyses of fault zone head waves and P ‐wave delay times at array A5 show that the Garlock fault is a sharp bimaterial interface with P waves traveling ∼5% faster in the northern crustal block. The across‐fault velocity contrast agrees with regional tomography models and generates clear P ‐wave reflections in waveforms recorded by array B4. Kirchhoff migration of the reflected waves indicates a near‐vertical fault between 2 and 6 km depth. The P ‐wave delay times imply a ∼300‐m‐wide transition zone near the Garlock fault surface trace beneath array A5, offset to the side with faster velocities. The results provide important constraints for derivations of earthquake properties, simulations of ruptures and ground motion, and future imaging studies associated with the Garlock fault.

58 GEOSCIENCES↗

DOE OSTI2023

Relocation of the Seismicity of the Caucasus Region

Our objective is to improve the view of the seismicity in the Caucasus region using instrumental data between 1951 and 2019. To create a comprehensive catalog, we combine the bulletins of local agencies and the International Seismological Centre, and use an advanced single-event location algorithm, iLoc, to obtain better locations. We show that relocations with iLoc, using travel-time predictions from the 3D upper mantle velocity model, Regional Seismic Travel Time, improve the locations. Then, using the iLoc results as initial locations and the ground-truth events identified in the iLoc results as fix points, we apply Bayesloc, a multiple-event location algorithm, to simultaneously relocate the entire seismicity of the Caucasus region. Here we demonstrate that the simultaneous relocation of the seismicity with Bayesloc clarifies the location and geometry of major active structures accommodating ongoing convergence between the Arabian and Eurasian continents between the Black and Caspian Seas. Among our major findings is the confirmation of widespread seismicity in the mantle beneath the northern flank of the Greater Caucasus and central Caspian, resulting from north-dipping subduction of the Kura and South Caspian basins and the identification of a discrete band of crustal seismicity beneath the southern flank of the Greater Caucasus.

58 GEOSCIENCES↗

DOE OSTI2025

Eocene-Oligocene Metamorphism, Fluid Flow and Deformation in the Ruby Mountains-East Humboldt Range Metamorphic Core Complex

Metamorphic core complexes are ubiquitous in collapsed orogens globally and play a primary role in crustal exhumation. Here, we investigate the metamorphic history of the Ruby Mountains-East Humboldt Range metamorphic core complex, Nevada, using petrochronology to understand how magmatism, metamorphism and deformation interact to modulate crustal rheology, and the timing of exhumation within the Sevier orogenic belt. Field study, microstructural analysis, thermobarometry and laser ablation split-stream monazite, titanite and allanite petrochronology of the mylonitic footwall were integrated to elucidate the pressure-temperature-time (P-T-t) evolution. Major-, trace-element and quartz-in-garnet thermobarometry show peak metamorphism occurred at 5.5–6 kbar and 600–650°C across the study area. Monazite and titanite U-(Th)-Pb petrochronology constrain this episode of metamorphism to 88–81 Ma with the structurally deepest samples yielding dates down to ca. 71 Ma. The preservation of peak metamorphic mineral assemblages indicates these rocks remained at near-peak temperature conditions into the Cenozoic. Cenozoic metamorphism occurred during a punctuated episode contemporaneous with the Great Basin ignimbrite flare up. All samples show a distinct 39–37 Ma population of monazite, titanite and allanite dates that correlate with the emplacement of quartz diorite and gabbro intrusions throughout the footwall. Accessory phases contain zoning textures in backscattered electron images that are consistent with fluid-mediated dissolution-reprecipitation reactions. The U-(Th)-Pb analyses from these zones yield 39–27 Ma dates interpreted as the timing of fluid flow. The results of this study elucidate the integrated importance of magmatism, melting and fluid flow in driving the coupled evolution of Late Cretaceous metamorphism within the Sevier hinterland and the subsequent Cenozoic metamorphic core complex development that exhumed the middle crust.

58 GEOSCIENCES↗

DOE OSTI2020

Viscosity of the Mercurian Magma Ocean: Implications of Sulfur-Free and Sulfur-Bearing Magma Oceans for Differentiation and Crustal Petrogenesis

Mercury is believed to have had a magma ocean in its early evolution whose solidification dynamics contributed to the planet’s present day surface composition. The surface is characterized by two main geological provinces; the Northern Smooth Plains (NSP),and the Heavily Cratered Terrain (HCT) and Intercratered Terrian (IcT). These provinces are defined by the differences in morphology, composition, and age. The NSP are located in the northern hemisphere and have low Mg/Si, Ca/Si, and S/Si ratios and high Al/Si ratios, and are compositionally analogous to a magnesian alkali basalt. The HCT/IcT cover the majority of the planet and are older than the NSP, and are characterized by high Mg/Si, Ca/Si, and S/Si ratios and low Al/Si ratios. The dichotomy in surface composition motivates the study into the structure of the mantle as a product of magma ocean solidification. Two important and interrelated parameters that determine the efficiency of crystal fractionation during magma ocean solidification are convection vigor and viscosity of the magma ocean. Crystal suspension in a cooling magma ocean could be efficient or inefficient which suggests two end member possibilities for magma ocean solidification. The first is the formation of a compositionally stratified mantle reflecting a cumulate sequence produced along a fractional liquid line of descent. The second is the formation of a well-mixed, batch crystallized mantle by high viscosity and/or vigorously convecting magma ocean. Determining which of these end members, or where on the spectrum between them is most relevant to Mercury’s magma ocean can be experimentally constrained through high P-T falling-sphere viscometry experiments.

Mouser, M. D.↗

DOE OSTI2020

Mineral reaction kinetics constrain the length scale of rock matrix diffusion

Mass transport by aqueous fluids is a dynamic process in shallow crustal systems, redistributing nutrients as well as contaminants. Rock matrix diffusion into fractures (void space) within crystalline rock has been postulated to play an important role in the transient storage of solutes. The reacted volume of host rock involved, however, will be controlled by fluid-rock reactions. Here we present the results of a study which focusses on defining the length scale over which rock matrix diffusion operates within crystalline rock over timescales that are relevant to safety assessment of radioactive and other long-lived wastes. Through detailed chemical and structural analysis of natural specimens sampled at depth from an active system (Toki Granite, Japan), we show that, contrary to commonly proposed models, the length scale of rock matrix diffusion may be extremely small, on the order of centimetres, even over timescales of millions of years. This implies that in many cases the importance of rock matrix diffusion will be minimal. Additional analyses of a contrasting crystalline rock system (Carnmenellis Granite, UK) corroborate these results.

58 GEOSCIENCES↗