Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Covariance Generation”

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 181 records · Page 10

Online Control of Process Variance Using Feedback

For many stochastic control problem formulations, design objectives typically include constraining the covariance of the outputs of the system. Recent engineering challenges have demonstrated an interest in controlling system output covariances to a setpoint. However, open-loop control of system output covariance can be exceptionally challenging when the system has significant plant uncertainties. This paper explores closed-loop control of the system’s output covariance using a windowed variance function. The stochastic nature of the windowed variance feedback means that root locus and other classical control design techniques cannot be immediately used. Expected value analysis is used to generate transfer functions for the average response of such systems, and the results are used to design controllers. Stochastic simulations are presented to verify the analytical results and make other observations concerning closed-loop variance control.

Bieniek, Mitchell↗

Detection of supernova magnitude fluctuations induced by large-scale structure

The peculiar velocities of supernovae and their host galaxies are correlated with the large-scale structure of the Universe, and can be used to constrain the growth rate of structure and test the cosmological model. In this work, we measure the correlation statistics of the large-scale structure traced by the Dark Energy Spectroscopic Instrument Bright Galaxy Survey Data Release 1 sample, and magnitude fluctuations of type Ia supernova from the Pantheon+ compilation across redshifts z < 0.1. We find a detection of the cross-correlation signal between galaxies and type Ia supernova magnitudes. Fitting the normalised growth rate of structure f sigma_8 to the auto- and cross-correlation function measurements we find f sigma_8 = 0.384 +0.094 -0.157, which is consistent with the Planck LambdaCDM model prediction, and indicates that the supernova magnitude fluctuations are induced by peculiar velocities. Using a large ensemble of N-body simulations, we validate our methodology, calibrate the covariance of the measurements, and demonstrate that our results are insensitive to supernova selection effects. We highlight the potential of this methodology for measuring the growth rate of structure, and forecast that the next generation of type Ia supernova surveys will improve f sigma_8 constraints by a further order of magnitude.

Nguyen, A. [Swinburne U., Ctr. Astrophys. Supercom↗

NGEE Arctic CO2, CH4 and Energy Eddy-Covariance (EC) Flux Tower Auxiliary Measurements, Utqiagvik (Barrow), Alaska, 2012 - 2019

Measurements began in September 2012 and are generally made throughout the snow-free growing season starting from April (pre-thawing) until November (freeze-up) when the instrumentation is taken down for maintenance and instrument calibration. CO2, CH4 and energy fluxes using the Eddy Covariance (EC) technique (Baldocchi 2003) and meteorological measurements are reported as 30-minute averages. The tower site is registered with AmeriFlux as US-NGB. See http://ameriflux.lbl.gov/sites/siteinfo/US-NGB for more information. Reported data include: (1) additional/extra meteorological and eddy covariance CO2, CH4 and energy flux data, (2) binned spectra for the three wind components, the sonic temperature and gas densities together with the binned cospectra for covariances of w (vertical wind component) and gas densities, (3) binned ogives (cumulative (co)spectra) for the three wind components, the sonic temperature and gas densities, and (4) tower site footprint matrices. Meteorological and eddy covariance CO2, CH4 and energy flux data (AmeriFlux BASE data product) are available from AmeriFlux https://doi.org/10.17190/AMF/1436326. The tower site is adjacent to the NGEE Arctic Tram system with several complementary time-series measurements https://doi.org/10.5440/1558782.This dataset supersedes the previously released dataset "Eddy-Covariance and auxiliary measurements, NGEE-Barrow, 2012-2013" DOI:10.5440/1124200 as this new product includes data for the longer time period processed using consistent methods for the entire record. The Next-Generation Ecosystem Experiments: Arctic (NGEE Arctic), was a research effort to reduce uncertainty in Earth System Models by developing a predictive understanding of carbon-rich Arctic ecosystems and feedbacks to climate. NGEE Arctic was supported by the Department of Energy's Office of Biological and Environmental Research. The NGEE Arctic project had two field research sites: 1) located within the Arctic polygonal tundra coastal region on the Barrow Environmental Observatory (BEO) near Utqiagvik (Barrow), Alaska and 2) multiple areas in the discontinuous permafrost region of the Seward Peninsula north of Nome, Alaska. Through observations, experiments, and synthesis with existing datasets, NGEE Arctic provided an enhanced knowledge base for multi-scale modeling and contributed to improved process representation at global pan-Arctic scales within the Department of Energy's Earth system Model (the Energy Exascale Earth System Model, or E3SM), and specifically within the E3SM Land Model component (ELM).

54 ENVIRONMENTAL SCIENCES↗

EOS attitude determination and next generation star tracker enhancements

The pointing knowledge required for the Earth Observing System (EOS) AM mission is at the limit of the current generation of star trackers, with little margin. Techniques for improving the performance of existing star trackers are explored, with performance sensitivities developed for each alternative. These are extended to define the most significant performance enhancements for a next generation star tracker. Since attitude determination studies tend to be computationally intensive, an approach for using a simpler one degree of freedom formulation is contrasted with a full three degree of freedom formulation. Additionally, covariance analysis results are compared with time domain simulation performance results.

Kudva, P.↗

Parametric Analysis of a Hover Test Vehicle using Advanced Test Generation and Data Analysis

Large complex aerospace systems are generally validated in regions local to anticipated operating points rather than through characterization of the entire feasible operational envelope of the system. This is due to the large parameter space, and complex, highly coupled nonlinear nature of the different systems that contribute to the performance of the aerospace system. We have addressed the factors deterring such an analysis by applying a combination of technologies to the area of flight envelop assessment. We utilize n-factor (2,3) combinatorial parameter variations to limit the number of cases, but still explore important interactions in the parameter space in a systematic fashion. The data generated is automatically analyzed through a combination of unsupervised learning using a Bayesian multivariate clustering technique (AutoBayes) and supervised learning of critical parameter ranges using the machine-learning tool TAR3, a treatment learner. Covariance analysis with scatter plots and likelihood contours are used to visualize correlations between simulation parameters and simulation results, a task that requires tool support, especially for large and complex models. We present results of simulation experiments for a cold-gas-powered hover test vehicle.

Gundy-Burlet, Karen↗

Performance Assessment of Different Pulse Reconstruction Algorithms for the ATHENA X-Ray Integral Field Unit

The X-ray Integral Field Unit (X-IFU) microcalorimeter, on-board Athena, with its focal plane comprising 3840 Transition Edge Sensors (TESs) operating at 90 mK, will provide unprecedented spectral-imaging capability in the 0.2-12 keV energy range. It will rely on the on-board digital processing of current pulses induced by the heat deposited in the TES absorber, as to recover the energy of each individual events. Assessing the capabilities of the pulse reconstruction is required to understand the overall scientific performance of the X-IFU, notably in terms of energy resolution degradation with both increasing energies and count rates. Using synthetic data streams generated by the X-IFU End-to-End simulator, we present here a comprehensive benchmark of various pulse reconstruction techniques, ranging from standard optimal filtering to more advanced algorithms based on noise covariance matrices. Beside deriving the spectral resolution achieved by the different algorithms, a first assessment of the computing power and ground calibration needs is presented. Overall, all methods show similar performances, with the reconstruction based on noise covariance matrices showing the best improvement with respect to the standard optimal filtering technique. Due to prohibitive calibration needs, this method might however not be applicable to the X-IFU and the best compromise currently appears to be the so-called resistance space analysis which also features very promising high count rate capabilities.

microcalorimeters↗

Tractable Learning in Underexcited Power Grids

Estimating the structure of physical flow networks, such as power grids, is critical to secure delivery of energy. This article discusses statistical structure estimation in power grids in the “underexcited” regime, where a subset of internal nodes has zero injection fluctuations. Prior estimation algorithms based on nodal voltages fail for such grids as the voltage covariance matrix is not invertible. Here, we propose a novel topology learning algorithm for learning underexcited general networks. Our algorithm uses physics-informed conservation laws to first identify the zero-injection buses and their neighbors, and then estimates the remaining edges in the grid. We prove the asymptotic correctness of our algorithm for grids with nonadjacent internal zero-injection nodes. More important, we theoretically analyze our algorithm’s efficacy under noisy measurements, and determine bounds on maximum noise under which asymptotically correct recovery is guaranteed. Our approach is validated through simulations with voltage samples generated on test distribution grids with real injection data and nonlinear power flow models.

24 POWER TRANSMISSION AND DISTRIBUTION↗

Covariant spectator theory of n p scattering: Deuteron form factors

The deuteron form factors are calculated using two model wave functions obtained from the 2007 CST high precision fits to $np$ scattering data. Included in the calculation are a new class of isoscalar $np$ interaction currents which are automatically generated by the nuclear force model used in these fits. If the nuclear model WJC2 is used, a precision fit ($\chi^2$/datum $\eqsim1$) to the Sick Global Analysis (GA) of all $ed$ elastic scattering data can be obtained by adjusting the unknown off-shell nucleon form factors $F_3(Q^2)$ (discussed before) and $F_4(Q^2)$ (introduced in this paper), and predicting the high $Q^2$ behavior of the neutron charge form factor $G_{En}(Q^2)$ well beyond the region where it has been measured directly. Relativistic corrections, isoscalar interaction currents, and off-shell effects are defined, discussed, and their size displayed. A rationale for extending $ed$ elastic scattering measurements to higher $Q^2$ is presented.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

X-BASE: the first terrestrial carbon and water flux products from an extended data-driven scaling framework, FLUXCOM-X

Mapping in situ eddy covariance measurements of terrestrial land–atmosphere fluxes to the globe is a key method for diagnosing the Earth system from a data-driven perspective. We describe the first global products (called X-BASE) from a newly implemented upscaling framework, FLUXCOM-X, representing an advancement from the previous generation of FLUXCOM products in terms of flexibility and technical capabilities. The X-BASE products are comprised of estimates of CO 2 net ecosystem exchange (NEE), gross primary productivity (GPP), evapotranspiration (ET), and for the first time a novel, fully data-driven global transpiration product (ETT), at high spatial (0.05°) and temporal (hourly) resolution. X-BASE estimates the global NEE at −5.75 ± 0.33 Pg C yr −1 for the period 2001–2020, showing a much higher consistency with independent atmospheric carbon cycle constraints compared to the previous versions of FLUXCOM. The improvement of global NEE was likely only possible thanks to the international effort to increase the precision and consistency of eddy covariance collection and processing pipelines, as well as to the extension of the measurements to more site years resulting in a wider coverage of bioclimatic conditions. However, X-BASE global net ecosystem exchange shows a very low interannual variability, which is common to state-of-the-art data-driven flux products and remains a scientific challenge. With 125 ± 2.1 Pg C yr −1 for the same period, X-BASE GPP is slightly higher than previous FLUXCOM estimates, mostly in temperate and boreal areas. X-BASE evapotranspiration amounts to 74.7×10 3 ± 0.9×10 3 km 3 globally for the years 2001–2020 but exceeds precipitation in many dry areas, likely indicating overestimation in these regions. On average 57 % of evapotranspiration is estimated to be transpiration, in good agreement with isotope-based approaches, but higher than estimates from many land surface models. Despite considerable improvements to the previous upscaling products, many further opportunities for development exist. Pathways of exploration include methodological choices in the selection and processing of eddy covariance and satellite observations, their ingestion into the framework, and the configuration of machine learning methods. For this, the new FLUXCOM-X framework was specifically designed to have the necessary flexibility to experiment, diagnose, and converge to more accurate global flux estimates.

Nelson, Jacob A.↗

Kalman filter estimation of attitude and gyro bias with the QUEST observation model

The loss function for the Wahba attitude estimation problem employs unit vector observations with scalar weights. It is usually associated with the QUEST observation model, where the actual sensor noise is assumed the same for all components of the observed vector, regardless of the position in the field of view. The QUEST model has the great advantage of being sensor-independent except for the scalar parameter that characterizes the sensor errors. Although efficient algorithms for solving the Wahba problem exist, extending these algorithms to estimate gyro biases or sensor alignments has had mixed success. However, it is straightforward to estimate bias and alignment parameters with a Kalman filter. This paper investigates the use of an extended Kalman filter for the attitude and gyro bias that incorporates the QUEST observation model, to be referred to as the Unit Vector Filter (UVF). The UVF results are compared with those from a more conventional filter, the Real-Time Sequential Filter (RTSF), for which the residual is the two-dimensional projection of the unit vector onto the plane perpendicular to the sensor boresight. The RTSF is similar in design to that used by Multimission Modular Spacecraft for onboard attitude determination. An apparent obstacle to the use of unit vectors as measurements is their singular noise covariance matrix. Shuster has shown that this problem should not affect filter performance. The UVF and RTSF are tested using actual flight data from the Extreme Ultraviolet Explorer (EUVE). It is found that these filters generate nearly identical attitude and gyro bias estimates, thus validating the use of unit vectors and the QUEST noise model.

Sedlak, J.↗

GeoNEX-SUBSETS: Cutouts of geostationary satellite data over long-term monitoring sites

Satellite remote sensing data are important tool for extrapolating the knowledge obtained at the Fluxnet towers. We introduce our GEO-NEX geostationary satellite subset products, which make it easy to compare between ground observation and Geo-NEX products. The MODIS subset has been frequently used for validating ecosystem models developed at flux towers and upscaling the observed flux data to regional scales. However, MODIS on the polar orbiting satellites can observe target regions only once a day, while the Fluxnet eddy-covariance data are compiled as sub-hourly datasets. As a results, summarizing the sub-hourly flux data into daily statistics is necessary for the comparison between MODIS and Flux data. Here, the new generation geostationary satellite sensors (GOES-16/17 ABI and Himawari-8/9 AHI) has the high-frequent observation feature (10 minutes) in addition to similar spectral band and spatial resolution with MODIS. The high frequent observation allows us to compare the flux diurnal cycle with geostationary satellite sensor data. Some studies have already shown the effective utilization of time series of geostationary satellite data for ecosystem modeling. We are producing NEX Level-1G products, which is the gridded Top-of-Atmosphere reflectance and brightness temperature data of geostationary satellite sensors. We cut out the NEX Level-1G data using the same file format with the MODIS subset except for the projection. The other data products (e.g., surface reflectance, land surface temperature, vegetation indices and climate data) will be also added upon their availability. We selected the ground observation sites from Fluxnet, Phenocam, and AERONET networks. The NEX subset data will be provided through NASA NEX data portal.

GeoNEX↗

A Graphical Model for Fusing Diverse Microbiome Data

This paper develops a Bayesian graphical model for fusing disparate types of count data. The motivating application is the study of bacterial communities from diverse high-dimensional features, in this case, transcripts, collected from different treatments. In such datasets, there are no explicit correspondences between the communities and each corresponds to different factors, making data fusion challenging. We introduce a flexible multinomial-Gaussian generative model for jointly modeling such count data. This latent variable model jointly characterizes the observed data through a common multivariate Gaussian latent space that parameterizes the set of multinomial probabilities of the transcriptome counts. The covariance matrix of the latent variables induces a covariance matrix of co-dependencies between all the transcripts, effectively fusing multiple data sources. We present a computationally scalable variational Expectation-Maximization (EM) algorithm for inferring the latent variables and the parameters of the model. Here, the inferred latent variables provide a common dimensionality reduction for visualizing the data and the inferred parameters provide a predictive posterior distribution. In addition to simulation studies that demonstrate the variational EM procedure, we apply our model to a bacterial microbiome dataset.

59 BASIC BIOLOGICAL SCIENCES↗

NGEE Arctic Soil Micro-warming Experiment Temperature Profiles, Council Road Mile Marker 71, Seward Peninsula, Alaska, 2017-2019

The reported soil temperature profile measurements (24 locations - 2 probes with 5 and 6 thermistors, respectively - Fig. 1 and Fig. 2) were initiated as part of a soil micro-warming experiment at the Council Road Mile Marker 71 Site (CN_MM71) in September 2017. From 2017 through August 2019, these were measurements of ambient pre-treatment plot conditions. This tussock tundra site (with underlying permafrost) experiences annual frost heaving that causes a vertical movement of the ground surface, hence causing some of the upper most temperature thermistors to be at surface level or above ground and recorded therefore air temperatures near the surface and not below ground soil temperature (see section on Quality Assurance). After August 2019, these measurements ended in preparation for transitioning to the experimental warming of individual plots. The reported temperature data are nominally 3-hour averages of measurements made at varying frequencies (3-hour maximum) with frequencies that depended upon power (solar) availability to operate the sensors (see section in documentation on Methods). The number of values and the standard deviation for each 3-hour average (where appropriate) are also provided. These measurements are located near the NGEE-Arctic CO2 and CH4 eddy covariance tower at the Council Road Site (US-NGC: NGEE Arctic Council, https://ameriflux.lbl.gov/sites/siteinfo/US-NGC#overview) and auxiliary data at https://doi.org/10.5440/1526749. This dataset contains four *.csv files and one *.pdf file. The Next-Generation Ecosystem Experiments: Arctic (NGEE Arctic), was a research effort to reduce uncertainty in Earth System Models by developing a predictive understanding of carbon-rich Arctic ecosystems and feedbacks to climate. NGEE Arctic was supported by the Department of Energy's Office of Biological and Environmental Research. The NGEE Arctic project had two field research sites: 1) located within the Arctic polygonal tundra coastal region on the Barrow Environmental Observatory (BEO) and the North Slope near Utqiagvik (Barrow), Alaska and 2) multiple areas on the discontinuous permafrost region of the Seward Peninsula north of Nome, Alaska. Through observations, experiments, and synthesis with existing datasets, NGEE Arctic provided an enhanced knowledge base for multi-scale modeling and contributed to improved process representation at global pan-Arctic scales within the Department of Energy's Earth system Model (the Energy Exascale Earth System Model, or E3SM), and specifically within the E3SM Land Model component (ELM).

54 ENVIRONMENTAL SCIENCES↗

Structure in the event-by-event energy-dependent neutron-$γ$ multiplicity correlations in 252 Cf(sf)

The emission of neutrons and γ rays by fission fragments reveal important information about the properties of fragments immediately following scission. The initial fragment properties, correlations between fragments, and emission competition give rise to correlations in neutron-γ emission. Recent theoretical and experimental advances have been proposed to explain the mechanism of angular momentum generation in fission, which would result in observable signature in neutron-γ emission correlations. In this paper, we present a novel analysis method of neutrons and γ rays emitted by fission fragments that allows us to discern structure in the observed correlations. We have analyzed data collected on 252 Cf (sf) at the Chi-Nu array at the Los Alamos Neutron Science Center. Through our analysis of the energy-differential neutron-γ multiplicity covariance, we have observed enhanced neutron-γ correlations, corresponding to rotational band γ-ray transitions, at γ-ray energies of 0.7 and 1.2 MeV. To shed light on the origin of this structure, we compare the experimental data with the predictions of three model calculations. The origin of the observed correlation structure is understood in terms of a positive spin-energy correlation in the generation of angular momentum in fission.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Coherent manipulation of graph states composed of finite-energy Gottesman-Kitaev-Preskill-encoded qubits

Graph states are a central resource in measurement-based quantum information processing. In the photonic qubit architecture based on Gottesman-Kitaev-Preskill (GKP) encoding, the generation of high-fidelity graph states composed of realistic, finite-energy approximate GKP-encoded qubits thus constitutes a key task. We consider the finite-energy approximation of GKP-qubit states given by a coherent superposition of shifted finite-squeezed vacuum states, where the displacements are Gaussian distributed. We present an exact description of graph states composed of such approximate GKP qubits as a coherent superposition of a Gaussian ensemble of randomly displaced ideal GKP-qubit graph states. Using standard Gaussian dynamics, we track the transformation of the covariance matrix and the mean-displacement vector elements of the Gaussian distribution of the ensemble under tools such as GKP-Steane error-correction and fusion operations that can be used to grow large high-fidelity GKP-qubit graph states. The covariance matrix elements capture the noise in the graph state due to the finite-energy approximation of GKP qubits, while the mean displacements relate to the possible absolute shift errors on the individual qubits arising conditionally from the homodyne measurements that are a part of these tools. Our work thus pins down an exact coherent error model for graph states generated from truly finite-energy GKP qubits, which can shed light on their error-correction properties.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Selection of a Pair of Experiments to Optimally Reduce Uncertainty in Targeted Nuclear Data

We propose a novel process to select a pair of differential and integral experiments that best reduce uncertainties in targeted 239 ⁢Pu nuclear data while compressing the current nuclear data pipeline from 20 to 3 years. 239⁢ Pu nuclear data are poorly understood for neutrons in the intermediate energy range due to sparsity and uncertainty in historical experiments. New experiments targeting this range will enable better understanding of these nuclear data, but choosing the ideal experiments to conduct is challenging. Beginning with a prior distribution represented by samples of nuclear data generated from theory, generalized least squares adjustments are made to incorporate data from historical experiments. To quantify potential uncertainty reduction obtainable from a pair of candidate experiments, we compute the D-optimality criterion of the posterior covariance of intermediate energy range nuclear data compared to the equivalent covariance after additional adjustment to the pair of candidate experiments. Repeating the process for each of many candidate pairs facilitates the final selection. Results support 63⁢ Cu total cross section measurements for differential experiments and alumina and alumina/graphite configurations for integral experiments. This analysis enables choosing differential and integral experiments to be executed concurrently while shortening decision times relative to the current nuclear data pipeline.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

The DESI-Lensing Mock Challenge: large-scale cosmological analysis of 3x2-pt statistics

The current generation of large galaxy surveys will test the cosmological model by combining multiple types of observational probes. Realising the statistical promise of these new datasets requires rigorous attention to all aspects of analysis including cosmological measurements, modelling, covariance and parameter likelihood. In this paper we present the results of an end-to-end simulation study designed to test the analysis pipeline for the combination of the Dark Energy Spectroscopic Instrument (DESI) Year 1 galaxy redshift dataset and separate weak gravitational lensing information from the Kilo-Degree Survey, Dark Energy Survey and Hyper-Suprime-Cam Survey. Our analysis employs the 3x2-pt correlation functions including cosmic shear and galaxy-galaxy lensing, together with the projected correlation function of the spectroscopic DESI lenses. We build realistic simulations of these datasets including galaxy halo occupation distributions, photometric redshift errors, weights, multiplicative shear calibration biases and magnification. We calculate the analytical covariance of these correlation functions including the Gaussian, noise and super-sample contributions, and show that our covariance determination agrees with estimates based on the ensemble of simulations. We use a Bayesian inference platform to demonstrate that we can recover the fiducial cosmological parameters of the simulation within the statistical error margin of the experiment, investigating the sensitivity to scale cuts. This study is the first in a sequence of papers in which we present and validate the large-scale 3x2-pt cosmological analysis of DESI-Y1.

79 ASTRONOMY AND ASTROPHYSICS↗

SNC Meteorites and Martian Reservoirs

Jones first suggested that the inverse covariation of initial epsilon (Nd-143) and Sr-87/Sr-86 of the shergottites could be explained by interaction between mantle-derived magmas with another isotopic reservoir(s) (i.e., assimilation or contamination). In that model, magmas were generated in a source region that was isotopically similar to the Nakhla source and the second reservoir(s) was presumed to be crust. The text also permitted the second reservoir to be another type of mantle, but I can confirm that a second mantle reservoir was never seriously considered by that author. Other features of this model were that (i) it occurred at a particular time, 180 m.y. ago, and (ii) the interacting reservoirs had been separated at approximately 4.5 b.y. In a later paper Jones explored this mixing model more quantitatively and concluded that magmas from a Nakhla-like source region at 180 m.y. would fall on or near an isotopic Nd-Sr-Pb hyperplane defined by the shergottites. This criterion was a necessary prerequisite for the parent magma(s) of the shergottites to have initially been Nakhla-like isotopically. At this juncture, it is perhaps worthwhile to note that this mixing model was not presented to explain geochemical variations but as a justification for a 180 m.y. crystallization age for the shergottites and a 1.3 b.y. crystallization age for the nakhlites. In the mid-1980's crystallization ages estimated for Nakhla ranged from approximately 1.3 b.y to 4.5 b.y. Similarly, preferred crystallization ages for the shergottites ranged from 360 m.y., to 1.3 b.y., to 4.5 b.y. In all these models, the 180 m.y. event seen in the shergottites was deemed to be metamorphic. The fit between the Nakhla-like source region and the shergottite hyperplane was a validation both of the 1.3 b.y. igneous age of Nakhla and the 180 m.y. igneous age of the shergottites.

Jones, J. H.↗