Search NASA⌕ Search

SEARCH · Search NASA

Results for “heterogeneous spatial distributions”

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 109 records · Page 6
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Uncovering Uranium Isotopic Heterogeneity of Fuel Pellets from the Fifth Collaborative Materials Exercise of the Nuclear Forensics International Technical Working Group

In 2017, the Nuclear Forensics International Technical Working Group (ITWG) organized their fifth 37 Collaborative Materials Exercise (CMX-5). The exercise samples were two uranium dioxide fuel pellets 38 manufactured from the same starting materials by different processes to have similar bulk isotopic 39 composition, but different spatial uranium isotopic distributions. Sets of identical materials were sent to 40 all participating laboratories, who then utilized their existing nuclear forensic capabilities to 41 independently analyse fuel pellets and identify similarities and differences of the materials’ 42 characteristics. The analytical methods used to probe the fuel pellets included ex situ, such as sectioning 43 or breaking up the pellets and analyzing dissolved pieces using inductively coupled plasma mass 44 spectrometry (ICP-MS), analyzing particles collected from intact or fragmented pellets by secondary ion 45 mass spectrometry (SIMS), as well as in situ methods, such as laser ablation coupled with ICP-MS, 46 autoradiography and nanoSIMS. In this paper we present the results of these independent analyses and 47 compare the capabilities of those nuclear forensic analytical methods to uncover details of the isotopic 48 heterogeneity of uranium fuel pellets.

Nuclear Forensic Analysis of Uranium Fuel Pellets,↗
DOE OSTI2020

Predictive Inverse Model for Advective Heat Transfer in a Short–Circuited Fracture: Dimensional Analysis, Machine Learning, and Field Demonstration

Identifying fluid flow maldistribution in planar geometries is a well–established problem in subsurface science/engineering. Of particular importance to the thermal performance of enhanced (or “engineered”) geothermal systems is identifying the existence of nonuniform (i.e., heterogeneous) permeability and subsequently predicting advective heat transfer. Here, machine learning via a genetic algorithm (GA) identifies the spatial distribution of an unknown permeability field in a two–dimensional Hele–Shaw geometry (i.e., parallel plates). The inverse problem is solved by minimizing the L2 norm between simulated residence time distribution (RTD) and measurements of an inert tracer breakthrough curve (BTC) (C–Dot nanoparticle). Principal component analysis (PCA) of spatially correlated permeability fields enabled reduction of the parameter space by more than a factor of 10 and restricted the inverse search to reservoir–scale permeability variations. Thermal experiments and tracer tests conducted at the mesoscale Altona Field Laboratory (AFL) demonstrate that the method accurately predicts the effects of extreme flow channeling on heat transfer in a single bedding–plane rock fracture. However, this is only true when the permeability distributions provide adequate matches to both tracer RTD and frictional pressure loss. Without good agreement to frictional pressure loss, it is still possible to match a simulated RTD to measurements, but subsequent predictions of heat transfer are grossly inaccurate. Here, the results of this study suggest that it is possible to anticipate the thermal effects of flow maldistribution, but only if both simulated RTDs and frictional pressure loss between fluid inlets and outlets are in good agreement with measurements.

42 ENGINEERING↗
NASA NTRS2018

Overview of Primitive Object Volatile Explorer (PrOVE) CubeSat or Smallsat Concept

Here we describe the Primitive Object Volatile Explorer (PrOVE), a smallsat mission concept to study the surface structure and volatile inventory of comets in their perihelion passage phase when volatile activity is near peak. CubeSat infrastructure imposes limits on propulsion systems, which are compounded by sensitivity to the spacecraft disposal state from the launch platform and potential launch delays. We propose circumventing launch platform complications by using waypoints in space to park a deep space SmallSat or CubeSat while awaiting the opportunity to enter a trajectory to flyby a suitable target. In our Planetary Science Deep Space SmallSat Studies (PSDS3) project, we investigated scientific goals, waypoint options, potential concept of operations (ConOps) for periodic and new comets, spacecraft bus infrastructure requirements, launch platforms, and mission operations and phases. Our payload would include two low-risk instruments: a visible image (VisCAM) for 5-10 m resolution surface maps; and a highly versatile multispectral Comet CAMera (ComCAM) will measure 1) H2O, CO2, CO, and organics non-thermal fluorescence signatures in the 2-5 μm MWIR, and 2) 7-10 and 8-14 μm thermal (LWIR) emission. This payload would return unique data not obtainable from ground-based telescopes and complement data from Earth-orbiting observatories. Thus, the PrOVE mission would (1) acquire visible surface maps, (2) investigate chemical heterogeneity of a comet nucleus by quantifying volatile species abundance and changes with solar insolation, (3) map the spatial distribution of volatiles and determine any variations, and (4) determine the frequency and distribution of outbursts.

Clark, Pamela S.↗
NASA NTRS2018

Overview of Primitive Volatile Explorer (PrOVE) Cubesat or Smallsat Concept

Here we describe the Primitive Object Volatile Explorer (PrOVE), a smallsat mission concept to study the surface structure and volatile inventory of comets in their perihelion passage phase when volatile activity is near peak. CubeSat infrastructure imposes limits on propulsion systems, which are compounded by sensitivity to the spacecraft disposal state from the launch platform and potential launch delays. We propose circumventing launch platform complications by using waypoints in space to park a deep space SmallSat or CubeSat while awaiting the opportunity to enter a trajectory to flyby a suitable target. In our Planetary Science Deep Space SmallSat Studies (PSDS3) project, we investigated scientific goals, waypoint options, potential concept of operations (ConOps) for periodic and new comets, spacecraft bus infrastructure requirements, launch platforms, and mission operations and phases. Our payload would include two low-risk instruments: a visible image (VisCAM) for 5-10 m resolution surface maps; and a highly versatile multispectral Comet CAMera (ComCAM) will measure 1) H2O, CO2, CO, and organics non-thermal fluorescence signatures in the 2-5 µm MWIR, and 2) 7-10 and 8-14 µm thermal (LWIR) emission. This payload would return unique data not obtainable from ground-based telescopes and complement data from Earth-orbiting observatories. Thus, the PrOVE mission would (1) acquire visible surface maps, (2) investigate chemical heterogeneity of a comet nucleus by quantifying volatile species abundance and changes with solar insolation, (3) map the spatial distribution of volatiles and determine any variations, and (4) determine the frequency and distribution of outbursts.

Zucherman, Aaron↗
NASA NTRS1995

The lognormal distribution as a model for bio-optical variability in the sea

The lognormal distribution is presented as a useful model for bio-optical variability at a variety of spatial and temporal scales. A parametric statistical framework is presented for using the lognormal model to assess the effects of heterogeneity and scale on closure. Variability at small scales affects but is unresolved by large-scale measurements. An assumed lognormal distribution allows one to integrate over small-scale variability to predict large-scale measurements. Examples are presented to demonstrate how knowledge of the variance can be incorporated into models to relate measurements made at different scales.

Campbell, Janet W.↗
NASA NTRSDate not supplied

Agent-based modeling of microbes in space

Space is tough on organisms. Microorganisms traveling to space experience stress from environmental features such as ionizing radiation and lack of normal microgravity; however, much remains unknown about the mechanisms by which those environmental features affect microbial physiology. Microbes experience changes in gravity not directly but rather through changes in their fluid environment, and deep-space particle radiation causes cell damage that is complex but rare. Computational modeling at the single-cell level (agent-based modeling) can allow us to probe the spatially heterogeneous processes that characterize space stresses, to gain insight into the relationships of microbial cells with their environments and with each other. Here we present two software packages for simulating microbial population dynamics in space conditions: CAMDLES and AMMPER. Microbes growing in liquid culture medium in the microgravity of an orbital space station experience a quiescent, poorly-mixed fluid environment. CAMDLES (CFD-DEM Artificial Microgravity Developments for Living Ecosystem Simulation) simultaneously simulates biological, chemical, and mechanical processes to predict microbial ecological dynamics in microgravity, and in the rotating culture vessels used to create an artificial microgravity environment in the lab. Initial results demonstrate that the growth of a cross-feeding microbial consortium, dependent on the exchange of soluble metabolites, is sensitive to the initial spatial distribution of cells, and grows differently in real versus artificial microgravity. Microbial populations exposed to deep-space radiation experience spatially and temporally heterogeneous damage from the traversal of high-energy particles. AMMPER (Agent-Based Model for Microbial Populations Exposed to Radiation) pairs a 3d model of energy deposition along a radiation particle track with a microbial population growth and damage model to predict the effects of localized radiation damage on population-level responses. It includes a user-friendly graphical interface. AMMPER results agree with experimental data indicating that indirect effects of radiation (reactive oxygen species generation, metabolic impairment) have a greater impact on microorganisms than direct effects (DNA damage).

Jessica A Lee↗
NASA NTRS2022

Agent-Based Modeling of Microbes in Space

Space is tough on organisms. Microorganisms traveling to space experience stress from environmental features such as ionizing radiation and lack of normal gravity, and much remains unknown about the mechanisms by which those environmental features affect microbial physiology. Microbes experience changes in gravity not directly but rather through changes in their fluid environment, and deep-space particle radiation causes cell damage that is complex but rare. Computational modeling at the single-cell level (agent-based modeling) can allow us to probe the spatially heterogeneous processes that characterize space stresses, to gain insight into the relationships of microbial cells with their environments and with each other. Here we present two software packages for simulating microbial population dynamics in space conditions: CAMDLES and AMMPER. Microbes growing in liquid culture medium in the microgravity of an orbital space station experience a quiescent, poorly-mixed fluid environment. CAMDLES (CFD-DEM Artificial Microgravity Developments for Living Ecosystem Simulation) simultaneously simulates biological, chemical, and mechanical processes to predict microbial ecological dynamics in microgravity, and in the rotating culture vessels used to create an artificial microgravity environment in the lab. Initial results demonstrate that the growth of a cross-feeding microbial consortium, dependent on the exchange of soluble metabolites, is sensitive to the initial spatial distribution of cells, and grows differently in real versus artificial microgravity. Microbial populations exposed to deep-space radiation experience spatially and temporally heterogeneous damage from the traversal of high-energy particles. AMMPER (Agent-Based Model for Microbial Populations Exposed to Radiation) pairs a 3d model of energy deposition along a radiation particle track with a microbial population growth and damage model to predict the effects of localized radiation damage on population-level responses. It includes a user-friendly graphical interface. AMMPER growth curves recapitulate experimental results, and allow comparison between direct effects (DNA damage) and indirect effects (reactive oxygen species generation, metabolic impairment) of radiation.

microbiology↗
DOE OSTI2025

Evaluating post-fire watershed response to varying burn severity and precipitation regimes using fully-distributed and integrated hydrologic models

Wildfires can cause significant changes in vegetation and soil, which may lead to increased surface runoff and soil erosion, thereby affecting water cycling within ecosystems. This study uses the Advanced Terrestrial Simulator (ATS), an integrated and fully distributed hydrologic model at the watershed scale, to examine post-fire hydrologic responses in selected watersheds with varying burn severities in the Pacific Northwest region of the United States. The model integrates surface overland flow, subsurface flow, and canopy biophysical processes. We have developed a new fire module in ATS to account for changes in soil hydraulic properties caused by fire in the topsoil layer. Modeling results indicate that, in the year following a high-severity burn, watershed-averaged evapotranspiration decreases by 25%. Additionally, post-fire peak flows increase by 18-29% in watersheds burned with medium to high severity due to changes in soil properties. Conversely, a low-severity burn results in less than a 1% increase in post-fire peak flow. Furthermore, a high-severity fire causes a 38% reduction in the infiltration rate within the affected watershed during the first post-fire wet season. Hypothetical numerical experiments with varying precipitation regimes after a high-severity fire show that post-fire peak flows can increase by 1-29% due to fire-induced changes in soil hydraulic properties. This study highlights the importance of using fully distributed hydrologic models to quantify disturbance-feedback loops, which are essential for understanding the complexities brought about by spatial heterogeneity in post-fire landscapes.

ATS↗
NASA NTRS1993

Bright soil units on Mars determined from ISM imaging spectrometer data

The lithology of bright Martian soil provides evidence for chemical and physical processes that have modified the planet's surface. Data from the ISM imaging spectrometer, which observed much of the equatorial region at a spatial resolution of approximately 22 km, cover the NIR wavelength range critical to ascertaining the presence and abundance of Fe-containing phases, hydroxylated silicates, and H2O in the bright soil. ISM data previously have revealed spatial variations in depth of the 3.0-microns H2O absorption suggesting differences in water content, a weak absorption at 2.2 microns indicative of metal-OH in phyllosilicate, and variations in the 1-micron Fe absorption indicative of differences in Fe mineralogy. This paper summarizes first results of a systematic investigation of spectral heterogeneity in bright soils observed by ISM. At least seven 'units' with distinctive properties were discriminated. Comparison of their spatial distributions with Viking data shows that they generally correspond with previously recognized morphologic, color, and thermal features. These correspondences and the units' spectral attributes provide evidence for lithologic differences between the soils in different geologic settings.

Murchie, Scott↗
NASA NTRS2004

Lunar and Planetary Science XXXV: Early Solar System Chronology

The following topics were presented in this report: Iron Isotopic Fractionation During Vacuum Evaporation of Molten W?stite and Solar Compositions; Mg Isotope Ratio Zonation in CAIs - New Constraints on CAI Evolution; Sm-Nd Systematics of Chondrites; AMS Measurement of Mg-24(He-3,p)Al-26 Cross Section, Implications for the Al-26 Production in the Early Solar System; On Early Solar System Chronology: Implications of an Initially Heterogeneous Distribution of Short-lived Radionuclides; Revisiting Extraterrestrial U Isotope Ratios; Helium-Shell Nucleosynthesis and Extinct Radioactivities; High Spatial Resolution Ion Microprobe Measurements Refine Chronology of Orgueil Carbonate Formation; and Calibration of the Galactic Cosmic Ray Flux.

Source record↗
DOE OSTI2023

Amino mapping: possibility to visualize amino-N compounds in the rhizosphere of Zea Mays L.

Abstract Understanding N uptake by plants, the N cycle, and their relationship to soil heterogeneity has generated a great deal of interest in the distribution of amino-N compounds in soil. Visualization of the spatial distribution of amino-N in soil can provide insights into the role of labile N in plant-microbial mechanisms of N acquisition and plant N uptake, but until now, it has remained technically challenging. Here, we describe a novel technique to visualize the amino-N distribution at the root-soil interface. The technique is based on time-lapse amino mapping (TLAM) using membranes saturated with the fluorogenic OPAME reagent (O-phthalaldehyde and β-mercaptoethanol). OPAME in the membrane reacts with organic compounds containing a NH 2 functional group at the membrane-soil interface, generating a fluorescent product visible under UV light and detectable by a digital camera. The TLAM amino-mapping technique was applied to visualize and quantify the concentration of amino-N compounds in the rhizosphere of maize (Zea MaysL.). A ten times greater amino-N concentration was detected in the rhizosphere compared to non-rhizosphere soil. The high content of amino-N was mainly associated with the root tips and was 3 times larger than the average amino-N content at seminal roots. The amino-N rhizosphere was 2 times broader around the root tips than around other parts of the roots. We concluded that TLAM is a promising approach for monitoring the fate of labile N in soils. However, the technique needs to be standardized for different soil types, plant species, and climate conditions to allow wider application.

Agriculture↗
NASA NTRS1982

Spatial variability of soil hydraulics and remotely sensed soil parameters

The development of methods to correctly interpret remotely sensed information about soil moisture and soil temperature requires an understanding of water and energy flow in soil, because the signals originate from the surface, or from a shallow surface layer, but reflect processes in the entire profile. One formidable difficulty in this application of soil physics is the spatial heterogeneity of natural soils. Earlier work has suggested that the heterogeneity of soil hydraulic properties may be described by the frequency distribution of a single scale factor. The sensitivity of hydraulic and energetic processes to the variation of this scale factor is explored with a suitable numerical model. It is believed that such an analysis can help in deciding how accurately and extensively basic physical properties of field soils need to be known in order to interpret thermal or radar waveband signals. It appears that the saturated hydraulic conductivity needs to be known only to its order of magnitude, and that the required accuracy of the soil water retention function is about 0.02 volume fraction. Furthermore, the results may be helpful in deciding how the total scene or view field, as perceived through a sensor, is composed from the actual mosaic of transient soil properties, such as surface temperature or surface soil moisture. However, the latter proposition presupposes a random distribution of permanent properties, a condition that may not be met in many instances, and no solution of the problem is apparent.

Lascano, R. J.↗
DOE OSTI2021

Spatial heterogeneity analyses of pore structure and mineral composition of Barnett Shale using X-ray scattering techniques

The spatial heterogeneity of pore structure and mineral composition at mu m to mm scales in mm to cm-sized shale samples complicates the laboratory measurements of porosity & pore size distribution and consequently reservoir quality evaluation. Traditional fluid invasion characterization methods typically take hours to days to get data and can only test relatively small samples (usually 2.54 cm diameter core plugs and mm-sized rock chips). Furthermore, the Xray scattering techniques at the Advanced Photon Source of the Argonne National Laboratory have the capability to non-destructively determine the pore structure and mineral composition of mm x mm area in a short time of seconds and to hold several to hundreds of mm 2 sized sample to conduct large area heterogeneity analysis in the short duration of hours. With the use of ultra-small angle X-ray scattering (USAXS), small-angle X-ray scattering (SAXS), and wide-angle X-ray scattering (WAXS), the porosity, surface area and pore size distribution of two Barnett Shale samples have been mapped out and the mineral compositions are determined qualitatively. Overall, the porosity and pore surface area of the Barnett Shale Sample A over a 36 cm2 area ranges from 1.70 % to 8.13 % and 8.87 m 2 /g to 40.2 m 2 /g respectively. The porosity and surface area of Sample B over a 24 cm 2 area varies from 3.15 % to 11.7 % and 18.3 m 2 /g to 60.2 m 2 /g. In Sample A, the mineral composition shows a shift from carbonate-rich to siliciclastic-rich over a distance of several millimeters. For Sample B, the mineral compositions are mainly siliciclastic and do not show an obvious lithological change. The direct observations of high spatial heterogeneity of pore structure and mineral composition in shale illustrate the utility of integrated X-ray scattering techniques to provide valuable insights into shale studies and reservoir evaluation.

58 GEOSCIENCES↗
DOE OSTI2022

Probing charge density in materials with atomic resolution in real space

The charge distribution in materials at the nanoscale can often explain the origin of macroscopic properties such as localized conductivity or the plasmonic response and illuminate more fundamental changes in the microscopic structure such as changes in chemical bonding characteristics. Previously, direct visualization of the charge density with high spatial resolution was often a missing link in the formation of structure–property relationships, especially in heterogeneous materials systems. Furthermore, recent advancements in microscopy technology have enabled researchers to visualize the charge distribution in materials down to subatomic length scales. In this Technical Review, we discuss the developments in high-resolution real-space charge distribution imaging using diffraction techniques and electron microscopy, with a focus on the recent advancement of four-dimensional scanning transmission electron microscopy, electron holography, and applications to materials interfaces.

Electronic properties and materials↗
DOE OSTI2022

Parameterizing Convective Organization Effects With a Moisture-PDF Approach in Climate Models: Concept and a Regional Case Simulation

We propose a parameterization scheme of convective organization effects based on a moisture-distribution approach, which can reflect aggregation of convective cells within a model grid as well as the interaction between convection and spatial heterogeneity in free-troposphere moisture. With this concept, convective cells in an aggregated state are surrounded by air that is moister than the grid-mean condition, which provides a “shielding” effect favorable for deeper convective updrafts. Such effects are represented quantitatively via utilizing a schematic diagram about the geometric interpretation of convective cluster, dry area, and their contact area, in which both the convective cluster size and dry area fraction increase during the development of convection. Our sensitivity analyses indicate that the new scheme performs well in capturing the mean precipitation features. Excluding the convective organization effects in the model leads to a considerable reduction in the simulated precipitation magnitude. Notably, there exists an inverse relationship between background mean moisture condition and sub-grid moisture variability for a given rain rate, suggesting that with a high spatial heterogeneity in free-troposphere moisture, convection can still develop and maintain its strength under relatively dry background conditions, consistent with the convection-permitting model simulation. Overall, our new parameterization of convective organization effects can successfully reproduce the relationship between precipitation and sub-grid moisture variability, which is a missing element in traditional convection parameterization schemes but important for the simulations of precipitation variability at various scales in climate models.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Transit Time Estimation in Catchments: Recent Developments and Future Directions

Abstract Water transit time is now a standard measure in catchment hydrological and ecohydrological research. The last comprehensive review of transit time modeling approaches was published 15+ years ago. But since then the field has largely expanded with new data, theory and applications. Here, we review these new developments with focus on water‐age‐balance approaches and data‐based approaches. We discuss and compare methods including StorAge‐Selection functions, well/partially mixed compartments, water age tracking through spatially distributed models, direct transit time estimates from controlled experiments, young water fractions, and ensemble hydrograph separation. We unify some of the heterogeneity in the literature that has crept in with these many new approaches, in an attempt to clarify the key differences and similarities among them. Finally, we point to open questions in transit time research, including what we still need from theory, models, field work, and community practice.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Envisioning urban environments resilient to vector-borne diseases: a protocol to study dengue in Vietnam

Transmitted primarily by Aedes aegypti (Ae. aegypti) and Aedes albopictus (Ae. albopictus), arboviral diseases pose a major global public health threat. Dengue, chikungunya, and zika are increasingly prevalent in Southeast Asia. Among other arboviruses, dengue and zika are becoming more common in Central and South America. Given human encroachment into previously uninhabited, often deforested areas, to provide new housing in regions of population expansion, conceptualizing built urban environments in a novel way is urgently needed to safeguard against the growing climate change-driven threat of vector-borne diseases. By understanding the spread from a One Health perspective, enhanced control and prevention can be achieved. This is particularly important considering that climate change is likely to significantly impact the persistence of ponded water where mosquitoes breed due to increasing temperature and shifting rainfall patterns with regard to magnitude, duration, frequency, and season. Models can incorporate aquatic mosquito stages and adult spatial dynamics when habitats are heterogeneously available, thereby including dispersal and susceptible-exposed-infected-recovered (SEIR) epidemiology. Coupled with human population distribution (density, locations), atmospheric conditions (air temperature, precipitation), and hydrological conditions (soil moisture distribution, ponding persistence in topographic depressions), modeling has improved predictive ability for infection rates. However, it has not informed interventional approaches from an urban environment perspective which considers the role of ponds/lakes that support green spaces, the density of population that enables rapid spread of disease, and varying micro-habitats for various mosquito stages under climate change. Here, for an example of dengue in Vietnam, a preventive and predictive approach to design resilient urban environments is proposed, which uses data from rapidly expanding metropolitan communities to learn continually. This protocol deploys computational approaches including simulation and machine learning/artificial intelligence, underpinned by surveillance and medical data for validation and adaptive learning. Its application may best inform urban planning in low-middle income countries in tropical zones where arboviral pathogens are prevalent.

60 APPLIED LIFE SCIENCES↗
NASA NTRS2005

Planetary Geophysics and Tectonics

The broad objective of this work is to improve understanding of the internal structures and thermal and stress histories of the solid planets by combining results from analytical and computational modeling, and geophysical data analysis of gravity, topography and tectonic surface structures. During the past year we performed two quite independent studies in the attempt to explain the Mariner 10 magnetic observations of Mercury. In the first we revisited the possibility of crustal remanence by studying the conditions under which one could break symmetry inherent in Runcorn's model of a uniformly magnetized shell to produce a remanent signal with a dipolar form. In the second we applied a thin shell dynamo model to evaluate the range of intensity/structure for which such a planetary configuration can produce a dipole field consistent with Mariner 10 results. In the next full proposal cycle we will: (1) develop numerical and analytical and models of thin shell dynamos to address the possible nature of Mercury s present-day magnetic field and the demise of Mars magnetic field; (2) study the effect of degree-1 mantle convection on a core dynamo as relevant to the early magnetic field of Mars; (3) develop models of how the deep mantles of terrestrial planets are perturbed by large impacts and address the consequences for mantle evolution; (4) study the structure, compensation, state of stress, and viscous relaxation of lunar basins, and address implications for the Moon s state of stress and thermal history by modeling and gravity/topography analysis; and (5) use a three-dimensional viscous relaxation model for a planet with generalized vertical viscosity distribution to study the degree-two components of the Moon's topography and gravity fields to constrain the primordial stress state and spatial heterogeneity of the crust and mantle.

Zuber, Maria↗