Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “sampling variability”

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 505 records · Page 28

Detecting Long-period Variability in the SDSS Stripe 82 Standards Catalog

We report the results of a search for long-period (100 < P < 600 days) periodic variability in the SDSS Stripe 82 standards catalog. The SDSS coverage of Stripe 82 enables such a search because there are on average 20 observations per band in ugriz bands for about one million sources, collected over about 6 yr, with a faint limit of r ~ 22 mag and precisely calibrated 1%–2% photometry. We calculated the periods of variable source candidates in this sample using the Lomb–Scargle periodogram and considered the three highest periodogram peaks in each of the gri filters as relevant. Only those sources with gri periods consistent within 0.1% were later studied. We use the Kuiper statistic to ensure uniform distribution of data points in phased light curves. We present five sources with the spectra consistent with quasar spectra and plausible periodic variability. This SDSS-based search bodes well for future sensitive large-area surveys, such as the Rubin Observatory Legacy Survey of Space and Time, which, due to its larger sky coverage (about a factor of 60) and improved sensitivity (~2 mag), will be more powerful for finding such sources.

79 ASTRONOMY AND ASTROPHYSICS↗

Computational methods for efficient structural reliability and reliability sensitivity analysis

This paper presents recent developments in efficient structural reliability analysis methods. The paper proposes an efficient, adaptive importance sampling (AIS) method that can be used to compute reliability and reliability sensitivities. The AIS approach uses a sampling density that is proportional to the joint PDF of the random variables. Starting from an initial approximate failure domain, sampling proceeds adaptively and incrementally with the goal of reaching a sampling domain that is slightly greater than the failure domain to minimize over-sampling in the safe region. Several reliability sensitivity coefficients are proposed that can be computed directly and easily from the above AIS-based failure points. These probability sensitivities can be used for identifying key random variables and for adjusting design to achieve reliability-based objectives. The proposed AIS methodology is demonstrated using a turbine blade reliability analysis problem.

Wu, Y.-T.↗

Biogeochemical cycling at the aquatic–terrestrial interface is linked to parafluvial hyporheic zone inundation history

Abstract. The parafluvial hyporheic zone combines the heightened biogeochemical and microbial interactions indicative of a hyporheic region with direct atmospheric/terrestrial inputs and the effects of wet–dry cycles. Therefore, understanding biogeochemical cycling and microbial interactions in this ecotone is fundamental to understanding biogeochemical cycling at the aquatic–terrestrial interface and to creating robust hydrobiogeochemical models of dynamic river corridors. We aimed to (i) characterize biogeochemical and microbial differences in the parafluvial hyporheic zone across a small spatial domain (6 lateral meters) that spans a breadth of inundation histories and (ii) examine how parafluvial hyporheic sediments respond to laboratory-simulated re-inundation. Surface sediment was collected at four elevations along transects perpendicular to flow of the Columbia River, eastern WA, USA. The sediments were inundated by the river 0, 13, 127, and 398 days prior to sampling. Spatial variation in environmental variables (organic matter, moisture, nitrate, glucose, % C, % N) and microbial communities (16S and internal transcribed spacer (ITS) rRNA gene sequencing, qPCR) were driven by differences in inundation history. Microbial respiration did not differ significantly across inundation histories prior to forced inundation in laboratory incubations. Forced inundation suppressed microbial respiration across all histories, but the degree of suppression was dramatically different between the sediments saturated and unsaturated at the time of sample collection, indicating a binary threshold response to re-inundation. We present a conceptual model in which irregular hydrologic fluctuations facilitate microbial communities adapted to local conditions and a relatively high flux of CO2. Upon rewetting, microbial communities are initially suppressed metabolically, which results in lower CO2 flux rates primarily due to suppression of fungal respiration. Following prolonged inundation, the microbial community adapts to saturation by shifting composition, and the CO2 flux rebounds to prior levels due to the subsequent change in respiration. Our results indicate that the time between inundation events can push the system into alternate states: we suggest (i) that, above some threshold of inundation interval, re-inundation suppresses respiration to a consistent, low rate and (ii) that, below some inundation interval, re-inundation has a minor effect on respiration. Extending reactive transport models to capture processes that govern such dynamics will provide more robust predictions of river corridor biogeochemical function under altered surface water flow regimes in both managed and natural watersheds.

58 GEOSCIENCES↗

The extreme ultraviolet variability of quasars

We study the extreme ultraviolet (EUV) variability (rest frame wavelengths 500–920 Å) of high-luminosity quasars using Hubble Space Telescope (HST) (low to intermediate redshift sample) and Sloan Digital sky Survey (SDSS) (high redshift sample) archives. The combined HST and SDSS data indicates a much more pronounced variability when the sampling time between observations in the quasar rest frame is >2×10{sup 7} s compared to <1.5×10{sup 7} s. Based on an excess variance analysis, for time intervals <2×10{sup 7} s in the quasar rest frame, 10% of the quasars (4/40) show evidence of EUV variability. Similarly, for time intervals >2×10{sup 7} s in the quasar rest frame, 55% of the quasars (21/38) show evidence of EUV variability. The propensity for variability does not show any statistically significant change between 2.5×10{sup 7} s and 3.16×10{sup 7} s (1 year). The temporal behavior is one of a threshold time interval for significant variability as opposed to a gradual increase on these timescales. A threshold timescale can indicate a characteristic spatial dimension of the EUV region. We explore this concept in the context of the slim disk models of accretion. We find that for rapidly spinning black holes, the radial infall time to the plunge region of the optically thin surface layer of the slim disk that is responsible for the preponderance of the EUV flux emission (primarily within 0–7 black hole radii from the inner edge of the disk) is consistent with the empirically determined variability timescale.

79 ASTRONOMY AND ASTROPHYSICS↗

An extended polynomial chaos expansion for PDF characterization and variation with aleatory and epistemic uncertainties

This paper presents an extended polynomial chaos formalism for epistemic uncertainties and a new framework for evaluating sensitivities and variations of output probability density functions (PDF) to uncertainty in probabilistic models of input variables. An ”extended” polynomial chaos expansion (PCE) approach is developed that accounts for both aleatory and epistemic uncertainties, modeled as random variables, thus allowing a unified treatment of both types of uncertainty. We explore in particular epistemic uncertainty associated with the choice of prior probabilistic models for input parameters. A PCE-based Kernel Density (KDE) construction provides a composite map from the PCE coefficients and germ to the PDF of quantities of interest (QoI). Here, the sensitivities of these PDF with respect to the input parameters are then evaluated. Input parameters of the probabilistic models are considered. By sampling over the epistemic random variable, a family of PDFs is generated and the failure probability is itself estimated as a random variable with its own PCE. Integrating epistemic uncertainties within the PCE framework results in a computationally efficient paradigm for propagation and sensitivity evaluation. Two typical illustrative examples are used to demonstrate the proposed approach.

Aleatory uncertainty↗

Data for Yield from Iowa’s first commercial miscanthus fields: implications of spatial variability for productivity and sustainability beyond research plots

This dataset contains biomass yield measurements and associated vegetation index data collected from commercial Miscanthus × giganteus fields in eastern Iowa during the 2022–2023 growing seasons. The data support the analyses presented in the article: “Yield From Iowa's First Commercial Miscanthus Fields: Implications of Spatial Variability for Productivity and Sustainability Beyond Research Plots.” We collected 105 ground-truth biomass samples from four mature commercial fields (>4 years old) covering 92.81 ha. Samples were taken from 3 m² quadrats that were hand-harvested in alignment with commercial harvest timing. Stem biomass (excluding leaves) was weighed, moisture-corrected, and converted to dry-matter yield expressed in Mg DM ha⁻¹. Sampling locations were selected to capture spatial variability visible in aerial imagery and were recorded using RTK GPS. Each biomass observation was paired with vegetation indices derived from high-resolution PlanetScope satellite imagery (3 m resolution). Images were acquired throughout the growing season, and indices were calculated to evaluate their ability to predict end-of-season biomass yield. Statistical and machine learning approaches were used to identify key predictors, and a linear regression model based on end-of-July Green Normalized Difference Vegetation Index (GNDVI) was developed and evaluated. This repository includes the data used in that modeling workflow. Management practices, economic data, full imagery time series, and additional methodological details are described in the associated publication and are not included here. The dataset consists of three comma-separated value (CSV) files: 1. Combine_Groundtruth_Yield_VI_22_23.csv This file contains ground-truth biomass yield measurements and associated key vegetation index values collected during the 2022 and 2023 growing seasons. Rows: 105 observations Columns: Year — Year of observation (2022 or 2023) Field — Field location identifier Sample_number — Unique sample identifier GNDVI_End_Jul — Green Normalized Difference Vegetation Index calculated at end of July GNDVI_End_Aug — Green Normalized Difference Vegetation Index calculated at end of August NDRE_End_Aug — Normalized Difference Red Edge index calculated at end of August Biomass_Stem_Yield_MgDM/ha — Measured stem biomass yield (megagrams dry matter per hectare) 2. trainData_GNDVI.csv This file contains the subset of observations used to train the predictive relationship between July GNDVI and biomass yield. Rows: 76 observations Columns: Unnamed: 0 — Row index retained from the original data processing workflow GNDVI_End_Jul — GNDVI at end of July Stem_Yield_MgDM/ha — Observed stem biomass yield (Mg DM ha⁻¹) 3. testData_GNDVI.csv This file contains the test dataset used to evaluate model performance. Rows: 29 observations Columns: Unnamed: 0 — Row index retained from the original data processing workflow GNDVI_End_Jul — GNDVI at end of July Predicted_Yield_MgDM/ha — Model-predicted stem biomass yield (Mg DM ha⁻¹) Observed_Yield_MgDM/ha — Measured stem biomass yield (Mg DM ha⁻¹)

Potential yield, yield gap, in-field management, y↗

X-Ray Emission of Ultraviolet Variable Active Galactic Nucleus Candidates

We analyze the X-ray properties of a sample of 23 high-probability active galactic nucleus (AGN) candidates with ultraviolet variability identified in Wasleske et al. Using data from the Chandra X-ray Observatory and XMM-Newton Observatory, we that find 11/23 nuclei are X-ray detected. We use spectral energy distribution modeling to compute star formation rates and show that the X-ray luminosities are typically in excess of the X-ray emission expected from star formation by at least an order of magnitude. Interestingly, this sample shows a diversity of optical spectroscopic properties. We explore possible reasons for why some objects lack optical spectroscopic signatures of black hole activity while still being UV-variable and X-ray bright. We find that host galaxy stellar emission and obscuration from gas and dust are potential factors. We study where this sample falls on relationships such as α OX - L 2500 and L X - L IR and find that some of the sample falls outside the typical scatter for these relations, indicating that they differ from the standard quasar population. With the diversity of optical spectroscopic signatures and varying impacts of dust and stellar emissions on our sample, these results emphasize the strength of variability in selecting the most complete set of AGN, regardless of other host galaxy properties.

79 ASTRONOMY AND ASTROPHYSICS↗

An omnipresent diversity and variability in the chemical composition of atmospheric functionalized organic aerosol

Abstract The atmospheric evolution of organic compounds encompasses many thousands of compounds with varying volatility, polarity, and water solubility. The molecular-level chemical composition of this mixture plays a major, yet uncertain, role in its transformations and impacts. Here we perform a non-targeted molecular-level intercomparison of functionalized organic aerosol from three diverse field sites and a chamber. Despite similar bulk composition, we report large molecular-level variability between multi-hour organic aerosol samples at each site, with 66 ± 13% of functionalized compounds differing between consecutive samples. Single precursor environmental laboratory chamber experiments and fully chemically-explicit modeling confirm this variability is due to changes in emitted precursors, chemical age, and/or oxidation conditions. These molecular-level results demonstrate greater compositional variability than is typically observed in less-speciated measurements, such as bulk elemental composition, which tend to show less daily variability. These observations should inform future field and laboratory studies, including assessments of the effects of variability on aerosol properties and ultimately the development of strategic organic aerosol parameterizations for air quality and climate models.

Ditto, Jenna C.↗

An Examination of the Space Weathering Patina of Lunar Rock 76015

Space weathering discussions have generally centered around soils but exposed rocks will also incur the effects of weathering. Rocks have much longer surface lifetimes than an individual soil grain and thus record a longer history of exposure. By studying the weathering products which have built up on a rock surface, we can gain a deeper perspective on the weathering process and better assess the relative importance of various weathering components. The weathered coating, or patina, of the lunar rock 76015 has been previously studied under SEM and also by TEM using ultramicrotome sample preparation methods. However, to really understand the products involved in creating these coatings, it is helpful to examine the patina in cross section, something which is now possible though the use of Focused Ion Beam (FIB) sample prep techniques, which allows us to preserve intact the delicate stratigraphy of the patina coating and provides a unique cross-sectional view of the space weathering process. Several samples have been prepared from the rock and the coatings are found to be quite variable in thickness and composition from one sample to the next.

Noble, S.↗

Homogenization of Satellite Based Hyperspectral Infrared Sounder Data To Build Long Term Climate Record

Building long term climate record using data from multiple hyperspectral Infrared (IR)sounders requires the homogenization of different data record to ensure the consistency andcontinuity. Such a requirement comes from two perspectives: 1) the need to adjust theoverlapping measurements of different sounders to ensure the radiometric consistency in thespectral radiance domain; 2) the need for a rigorously defined scheme to ensure the radiometricconsistency being transferred to the essential climate variables derived from the radiance record.We develop a solution that uses a spectral fingerprinting scheme to derive anomalies of keyclimate variables from long term spectral radiance data record constructed using both AIRS andCrIS observations aboard AQUA, SNPP and JPSS satellites. The fingerprinting scheme usescommon radiative kernels for all sounder measurements and therefore effectively avoids thealgorithm introduced inconsistency in retrieved geophysical variables. Our approach uses aunified sampling scheme to match AIRS and CrIS observations in both spectral and spatial-temporal domain, facilitating the intercomparison of spectral radiances from different sensors(platforms). The optimized liner inversion scheme allows the direct quantification and thereforethe adjustment for the impact on the derived climate anomalies imposed by potential radiometricinconsistency between the overlapping measurements. Such a scheme also enables the low-latency data processing of long term hyperspectral sounder data records. This paper provides a detailed introduction of the spectral fingerprinting methodology.Also introduced here is the climate fingerprinting Sounder Product (ClimFiSP) developed basedon the fingerprinting methodology. ClimFiSP products include the space-time averagedproperties of key climate variables that are derived from the long-term, space-time averagedradiances from AQUA-AIRS, SNPP-CrIS, and JPSS1-CrIS. ClimFiSP will be available to usersthrough NASA'sGoddard Earth Sciences Data and Information Services Center (GES DISC).

Wan Wu↗

Efficiently sampling conformations and pathways using the concurrent adaptive sampling (CAS) algorithm

Molecular dynamics (MD) simulations are useful in obtaining thermodynamic and kinetic properties of bio-molecules but are limited by the timescale barrier, i.e., we may be unable to efficiently obtain properties because we need to run microseconds or longer simulations using femtoseconds time steps. While there are several existing methods to overcome this timescale barrier and efficiently sample thermodynamic and/or kinetic properties, problems remain in regard to being able to sample un- known systems, deal with high-dimensional space of collective variables, and focus the computational effort on slow timescales. Hence, in this work, a new sampling method, called the “Concurrent Adaptive Sampling (CAS) algorithm,” has been developed to tackle these three issues and efficiently obtain conformations and pathways. The method is not constrained to use only one or two collective variables, unlike most reaction coordinate-dependent methods. Instead, it can use a large number of collective vari- ables and uses macrostates (a partition of the collective variable space) to enhance the sampling. The exploration is done by running a large number of short simula- tions, and a clustering technique is used to accelerate the sampling. In this paper, we introduce the new methodology and show results from two-dimensional models and bio-molecules, such as penta-alanine and triazine polymer

59 BASIC BIOLOGICAL SCIENCES↗

Whole genome sequencing of Mycobacterium bovis directly from clinical tissue samples without culture

Advancement in next generation sequencing offers the possibility of routine use of whole genome sequencing (WGS) for Mycobacterium bovis (M. bovis) genomes in clinical reference laboratories. To date, the M. bovis genome could only be sequenced if the mycobacteria were cultured from tissue. This requirement for culture has been due to the overwhelmingly large amount of host DNA present when DNA is prepared directly from a granuloma. To overcome this formidable hurdle, we evaluated the usefulness of an RNA-based targeted enrichment method to sequence M. bovis DNA directly from tissue samples without culture. Initial spiking experiments for method development were established by spiking DNA extracted from tissue samples with serially diluted M. bovis BCG DNA at the following concentration range: 0.1 ng/μl to 0.1 pg/μl (10 –1 to 10 –4 ). Library preparation, hybridization and enrichment was performed using SureSelect custom capture library RNA baits and the SureSelect XT HS2 target enrichment system for Illumina paired-end sequencing. The method validation was then assessed using direct WGS of M. bovis DNA extracted from tissue samples from naturally (n = 6) and experimentally (n = 6) infected animals with variable Ct values. Direct WGS of spiked DNA samples achieved 99.1% mean genome coverage (mean depth of coverage: 108×) and 98.8% mean genome coverage (mean depth of coverage: 26.4×) for tissue samples spiked with BCG DNA at 10 –1 (mean Ct value: 20.3) and 10 –2 (mean Ct value: 23.4), respectively. The M. bovis genome from the experimentally and naturally infected tissue samples was successfully sequenced with a mean genome coverage of 99.56% and depth of genome coverage ranging from 9.2× to 72.1×. The spoligoyping and M. bovis group assignment derived from sequencing DNA directly from the infected tissue samples matched that of the cultured isolates from the same sample. Our results show that direct sequencing of M. bovis DNA from tissue samples has the potential to provide accurate sequencing of M. bovis genomes significantly faster than WGS from cultures in research and diagnostic settings.

59 BASIC BIOLOGICAL SCIENCES↗

Using low volume eDNA methods to sample pelagic marine animal assemblages

Environmental DNA (eDNA) is an increasingly useful method for detecting pelagic animals in the ocean but typically requires large water volumes to sample diverse assemblages. Ship-based pelagic sampling programs that could implement eDNA methods generally have restrictive water budgets. Studies that quantify how eDNA methods perform on low water volumes in the ocean are limited, especially in deep-sea habitats with low animal biomass and poorly described species assemblages. Using 12S rRNA and COI gene primers, we quantified assemblages comprised of micronekton, coastal forage fishes, and zooplankton from low volume eDNA seawater samples (n = 436, 380–1800 mL) collected at depths of 0–2200 m in the southern California Current. We compared diversity in eDNA samples to concurrently collected pelagic trawl samples (n = 27), detecting a higher diversity of vertebrate and invertebrate groups in the eDNA samples. Differences in assemblage composition could be explained by variability in size-selectivity among methods and DNA primer suitability across taxonomic groups. The number of reads and amplicon sequences variants (ASVs) did not vary substantially among shallow (<200 m) and deep samples (>600 m), but the proportion of invertebrate ASVs that could be assigned a species-level identification decreased with sampling depth. Using hierarchical clustering, we resolved horizontal and vertical variability in marine animal assemblages from samples characterized by a relatively low diversity of ecologically important species. Low volume eDNA samples will quantify greater taxonomic diversity as reference libraries, especially for deep-dwelling invertebrate species, continue to expand.

59 BASIC BIOLOGICAL SCIENCES↗

CARBON AND OXYGEN ISOTOPIC RATIOS FOR NEARBY MIRAS

Carbon and oxygen isotopic ratios are reported for a sample of 46 Mira and SRa-type variable asymptotic giant branch (AGB) stars. Vibration–rotation first and second-overtone CO lines in 1.5–2.5 μ m spectra were measured to derive isotopic ratios for {sup 12}C/{sup 13}C, {sup 16}O/{sup 17}O, and {sup 16}O/{sup 18}O. Comparisons with previous measurements for individual stars and with various samples of evolved stars, as available in the extant literature, are discussed. Models for solar composition AGB stars of different initial masses are used to interpret our results. We find that the majority of M-stars have main sequence masses ≤2 M {sub ⊙} and have not experienced sizable third dredge-up (TDU) episodes. The progenitors of the four S-type stars in our sample are slightly more massive. Of the six C-stars in the sample three have clear evidence relating their origin to the occurrence of TDU. Comparisons with O-rich presolar grains from AGB stars that lived before the formation of the solar system reveal variations in the interstellar medium chemical composition. The present generation of low-mass AGB stars, as represented by our sample of long period variables (LPVs), shows a large spread of {sup 16}O/{sup 17}O ratios, similar to that of group 1 presolar grains and in agreement with theoretical expectations for the composition of mass 1.2–2 M {sub ⊙} stars after the first dredge-up. In contrast, the {sup 16}O/{sup 18}O ratios of present-day LPVs are definitely smaller than those of group 1 grains. This is most probably a consequence of the the decrease with time of the {sup 16}O/{sup 18}O ratio in the interstellar medium due to the chemical evolution of the Milky Way. One star in our sample has an O composition similar to that of group 2 presolar grains originating in an AGB star undergoing extra-mixing. This may indicate that the extra-mixing process is hampered at high metallicity, or, equivalently, favored at low metallicity. Similarly to O-rich grains, no star in our sample shows evidence of hot bottom burning, which is expected for massive AGB stars.

79 ASTRONOMY AND ASTROPHYSICS↗

Cloud, Aerosol and Boundary Layer Structure across the Northeast Pacific Stratocumulus-Cumulus Transition as observed during CSET

During the Cloud System Evolution in the Trades (CSET) field study, 14 research flights of the National Science Foundation G-V sampled the stratocumulus-cumulus transition between Northern California and Hawaii and its synoptic variability. The G-V made vertically resolved measurements of turbulence, cloud microphysics, aerosol characteristics, and trace gases. It also carried dropsondes and a vertically pointing W-band radar and lidar. This paper summarizes these observations with the goals of fostering novel comparisons with theory, models and reanalyses, and satellite-derived products. A longitude-height binning and compositing strategy mitigates limitations of sparse sampling and spatiotemporal variability. Typically, a 1-km-deep decoupled stratocumulus-capped boundary layer near California evolved into 2-km-deep precipitating cumulus clusters surrounded by patches of thin stratus that dissipated toward Hawaii. Low cloud cover was correlated with estimated inversion strength more than with cloud droplet number, even though the thickest clouds were generally precipitating and ultraclean layers indicative of aerosol-cloud-precipitation interaction were common west of 140 degrees W. Accumulation-mode aerosol concentration correlated well with collocated cloud droplet number concentration and was typically largest near the surface. Aitken mode aerosol concentration was typically larger in the free troposphere. Wildfire smoke produced spikes of aerosol and trace gases on some flights. CSET data are compared with space-time collocated output from MERRA-2 reanalysis and from the CAM6 climate model run with winds and temperature nudged toward this reanalysis. The reanalysis compares better with the observed relative humidity than does nudged CAM6. Both vertically diffuse the stratocumulus cloud layer versus observations. MERRA-2 slightly underestimates in situ carbon monoxide measurements and underestimates ozone depletion within the boundary layer.

Bretherton, Christopher S.↗

The x ray variability of NGC6814: Power spectra

Simulation techniques are used to obtain the X-ray variability power spectrum of unevenly sampled GINGA data from NGC6814. A simple power law is not an adequate description of the power spectrum, with the residuals showing excess power on timescales consistent with the periodicity seen in EXOSAT observations of this object. However the shape of the folded lightcurve is very different, with 3 main peaks, two of which are separated by an extremely sharp dip instead of the single peak and small harmonic structure observed by EXOSAT. Using the dip as a fiducial mark, a second GINGA observation of this source taken one year later is found to be consistent with being completely periodic and phase coherent with this first GINGA observation. Thus the period is consistent with being constant over a period of 6 years, but phase coherence is only maintained on timescales of approximately 1 year. Over 75 percent of the total source variability is due to the periodic component (r.m.s. amplitude of 36 percent). The residual variability can be described as the more usual 'flicker noise' f(exp -1.1) powerlaw. This shows no apparent high frequency break on timescales greater than 300 seconds. Subtle differences in the shape of the folded light curve with energy, and the very large amount of power in the periodic component suggest occultation as its origin, though amplification of variability from an X-ray emitting 'hot spot' at the disk inner radius through gravitational lensing is also possible. The former suffers from the very arbitrary nature of the periodic timescale, while the latter is unattractive as it cannot simply explain the lack of high frequency break in the residual power. That these models probably fail to provide an adequate explanation may be due to the added complexity of anisotropy of the X-ray emission, suggested by the discrepancy between the lack of soft photons implied by the flat spectrum and the copious source of soft photons available from reprocessing in the iron line producing material.

Done, C.↗

The Effects of High-Temperature Fractional Crystallization on Calcium Isotopic Composition

We report Ca isotope fractionation among a co-genetic suite of samples from the Mauritanian Guelb el Azib ultramafic-mafic-anorthosite complex (GAC), which represents the fractional crystallization sequence of an igneous magma chamber. We investigate how the composition of an evolving crystallizing liquid and the resultant mineralogy of co-crystallizing assemblages affects inter-phase Ca isotope fractionation. Because the GAC is an Archean metamorphosed complex, we first investigate the resistance of the Ca isotope signature to secondary hydrothermal alteration and to metamorphism. While we cannot conclude here that Ca isotopes are undisturbed at the mineral scale, we show that they are preserved at the bulk rock scale. This is adequate for our study where several samples are almost monomineralic at the bulk scale level. The δ(exp 44/40)Ca (relative to NIST standard SRM 915a) of GAC layers range from −1.53 to 1.61‰, with the earliest-formed ultramafic cumulate layer being isotopically heaviest and the later, more felsic layers being isotopically lighter. For the first time, we show that plagioclase is much more enriched in light Ca isotopes compared to olivine and (Ca-) pyroxene. Monte Carlo simulations suggest that Ca isotopes are fractionated among co-existing silicate minerals during fractional crystallization, possibly in relation with the residual liquid composition. In qualitative agreement with computational models based on first principles lattice dynamics, we observe that Ca isotope fractionation is mineralogically controlled and importantly, that the degree of fractionation can vary according to the CaO composition of the residual liquid. While previous studies have aimed to understand partial melting as a source of Ca isotope fractionation, our results suggest that fractional crystallization is also a source of Ca isotope variability among co-genetic samples.

Maria C. Valdes↗