Search NASA⌕ Search

SEARCH · Search NASA

Results for “samples”

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 487 records · Page 27

Achieving Accuracy Requirements for Forest Biomass Mapping: A Spaceborne Data Fusion Method for Estimating Forest Biomass and Lidar Sampling Error

The synergistic use of active and passive remote sensing (i.e., data fusion) demonstrates the ability of spaceborne light detection and ranging (LiDAR), synthetic aperture radar (SAR) and multispectral imagery for achieving the accuracy requirements of a global forest biomass mapping mission (+/-20 Mg/ha or 20%, the greater of the two, for at least 80% of grid cells). A data fusion approach also provides a means to extend 3D information from discrete spaceborne LiDAR measurements of forest structure across scales much larger than that of the LiDAR footprint. For estimating biomass, these measurements mix a number of errors including those associated with LiDAR footprint sampling over regional-global extents. A general framework for mapping above ground live forest biomass density (AGB) with a data fusion approach is presented and verified using data from NASA field campaigns near Howland, ME, USA, to assess AGB and LiDAR sampling errors across a regionally representative landscape. We combined SAR and Landsat-derived optical (passive optical) image data to identify contiguous areas (>0.5 ha) that are relatively homogenous in remote sensing metrics (forest patches). We used this image-derived data with simulated spaceborne LiDAR derived from orbit and cloud cover simulations and airborne data from NASA's Laser Vegetation Imaging Sensor (LVIS) to compute AGB and estimate LiDAR sampling error for forest patches and 100 m, 250 m, 500 m, and 1 km grid cells. At both the patch and grid scales, we evaluated differences in AGB estimation and sampling error from the combined use of LiDAR with both SAR and passive optical and with either SAR or passive optical alone. First, this data fusion approach demonstrates that incorporating forest patches into the AGB mapping framework can provide sub-grid forest information for coarser grid-level AGB reporting. Second, a data fusion approach for estimating AGB using simulated spaceborne LiDAR with SAR and passive optical image combinations reduced forest AGB sampling errors 12%-38% from those where LiDAR is used with SAR or passive optical alone. In absolute terms, sampling errors were reduced from 14-40 Mg/ha to 11-28 Mg/ha across all grid scales and prediction methods, where minimum sampling errors were 11, 15, 18, and 22 Mg/ha for 1 km, 500 m, 250 m, and 100 m grid scales, respectively. Third, spaceborne global scale accuracy requirements were achieved whereby at least 80% of the grid cells at 100 m, 250 m, 500 m, and 1 km grid levels met AGB accuracy requirements using a combination of passive optical and SAR along with machine learning methods to predict vegetation structure metrics for forested areas without LiDAR samples. Finally, using either passive optical or SAR, accuracy requirements were met at the 500 m and 250 m grid level, respectively..

LiDAR↗

Active and Passive Storage Solutions for Low Temperature Lunar Sample Return

Active and passive thermal storage solutions to return lunar biological, physical science and/or geology samples to the Earth are under development. Temperature requirements for biological/physical science sample conditioning range from -100˚C to -153˚C and possibly down to -253˚C for geological samples. The lower limit for geological samples is derived from the temperature of permanently shadowed regions on the lunar surface and the sublimation temperature of specific volatiles of interest (i.e., H2O, NH3 and CH4). Individual samples for the lunar application are expected to be less than 50 kg with a 700 kg allowance for the container and refrigeration. Requiring no electrical power or heat rejection, passive approaches, comprised of high performance insulation and consumable Joule-Thompson cooling, may be preferred for shorter duration missions (< 30 days) to provide significant mass savings. Active storage approaches with cryogenic cooling may be necessary to preserve samples for longer periods of time. A notional passive storage concept with an internal vapor cooled shield is shown below. The vapor cooled shield contains the sample and is isolated from the outer container with concentric reflective rigid shields and conventional multi-layer and/or aerogel insulation on the outer layer. A special removable, insulating end cap to stow or retrieve the sample is included.

Lunar Sample Return↗

The Planetary Protection Strategy of Mars Sample Return’s Earth Return Orbiter Mission

The Mars Sample Return campaign aims to use three flight missions and one ground element to safely bring rock cores, regolith and atmospheric samples from the surface of Mars to Earth to answer key questions about the geologic and climate history of Mars, including the potential for ancient life. Since its landing in Jezero Crater in 2021, the first mission, NASA’s Mars 2020, has collected a number of samples on the crater floor and on the delta using the Perseverance rover. Subsequent missions would recover the sealed sample tubes, launch them into Mars orbit, and transport them back to Earth. The ground element would be a high-containment facility that would isolate and protect the samples during initial sample characterization, which would include sample safety assessments and time-sensitive scientific investigations. These elements are currently in the planning and design stages of development, and represent an international effort of NASA, the European Space Agency (ESA), and many industry partners. The work presented here provides an overview of the planetary protection strategy of the third flight mission, the ESA-led Earth Return Orbiter (ERO), which hosts the NASA-provided Capture, Containment, and Return System (CCRS). ERO-CCRS would capture the container with up to 30 sealed tubes previously put in Martian orbit, contain them in redundant containers to ensure that no potentially hazardous Mars particles are released, and return them to Earth through an entry vehicle. Both NASA and ESA policies comply with the United Nations’ Outer Space Treaty by planning to protect Earth’s biosphere from any potential adverse effects from material returned from solar system bodies beyond the Earth-Moon system. In the conduct of Mars Sample Return, the two agencies have agreed to apply approaches consistent with their own planetary protection standards to the campaign elements each provides.

Mars Sample Return↗

The Planetary Protection Strategy of Mars Sample Return Earth Return Orbiter Mission

The Mars Sample Return campaign aims to use three flight missions and one ground element to safely bring rock cores, regolith and atmospheric samples from the surface of Mars to Earth to answer key questions about the geologic and climate history of Mars, including the potential for ancient life. Since its landing in Jezero Crater in 2021, the first mission, NASA’s Mars 2020, has collected a number of samples on the crater floor and on the delta using the Perseverance rover. Subsequent missions would recover the sealed sample tubes, launch them into Mars orbit, and transport them back to Earth. The ground element would be a high-containment facility that would isolate and protect the samples during initial sample characterization, which would include sample safety assessments and time-sensitive scientific investigations. These elements are currently in the planning and design stages of development, and represent an international effort of NASA, the European Space Agency (ESA), and many industry partners. The work presented here provides an overview of the Planetary Protection strategy of the third flight mission, the ESA-led Earth Return Orbiter (ERO), which hosts the NASA-provided Capture, Containment, and Return System (CCRS). ERO-CCRS would detect and capture the container with up to 30 sealed tubes previously put in Martian orbit, contain them in redundant containers to ensure that no potentially hazardous Mars particles are released, and return them to Earth through an entry vehicle. Both NASA and ESA policies comply with the United Nations’ Outer Space Treaty by planning to protect Earth’s biosphere from any potential adverse effects from material returned from solar system bodies beyond the Earth-Moon system. In the conduct of Mars Sample Return, the two agencies have mutually agreed to apply approaches consistent with their own planetary protection standards to the campaign elements they each provides.

mars sample return↗

Investigating the Physical Modification of the Bennu Sample During Entry, Descent, and Landing

On September 24, 2023, the OSIRIS-REx Sample Return Capsule (SRC) entered Earth’s atmosphere and landed in the Utah Test and Training Range (UTTR). Preliminary examination of the returned Bennu sample has confirmed that OSIRIS-REx sample mass exceeds the mission requirement of 60 g of material. The sample consists of particles that range from a few centimeters to microscopic fines. During the SRC’s entry, descent, and landing (EDL) sequence, it may have experienced (i) peak decelerations of 10s of g (ii) tumbling, and (iii) touchdown at approximately 10 m/s, which could have induced physical modification of the sample. In addition, the act of sampling may have altered or biased the physical properties of the collected materials. Here, we investigate the likelihood and extent of physical modification of the sample between collection and return using observations and modeling. This work addresses the mission’s hypothesis 12, which concerns, in part, the modification of the sample during collection and Earth entry.

asteroid↗

The Scientific Value of Collecting Samples From the Jezero Crater Rim

The Mars 2020 mission has been conducting ground-based investigation of the geology, habitability, and biosignature preservation potential and collecting samples for return to Earth in Jezero crater, Mars for nearly 3 years. Analysis of these samples will address outstanding questions in Mars science including potential habitability and how and why the climate the interior of the planet evolved through time. As of December 2023, samples of 4 igneous rocks of the Jezero crater floor and 9 sedimentary rocks of the Jezero fan and inner margin, remain on the rover. 15 tubes remain to be filled to enhance the diversity of the cache and broaden the scope of the science questions that can be addressed with returned sample studies. The next step in the mission is to explore the Jezero crater rim. It will be imperative to investigate and sample the diversity of crater rim rocks because they represent materials from Mars’ most ancient crust older than those sampled in Jezero crater, a diversity of geologic processes, and potential ancient habitable environments that have not yet been investigated or sampled. Ongoing mapping efforts are using orbiter data and long-distance images from Perseverance to identify and interpret the geologic context of the crater rim. Building on this effort and the broader geologic context for the crater rim put forward by previous studies, we identify diverse targets for in situ investigation and potential sampling by Mars 2020.

Mars sample return↗

The Planetary Protection Strategy of Mars Sample Return’s Earth Return Orbiter Mission

The Mars Sample Return campaign aims to use three flight missions and one ground element to safely bring rock cores, regolith and atmospheric samples from the surface of Mars to Earth to answer key questions about the geologic and climate history of Mars, including the potential for ancient life. Since its landing in Jezero Crater in 2021, the first mission, NASA’s Mars 2020, has collected a number of samples on the crater floor and on the delta using the Perseverance rover. Subsequent missions would recover the sealed sample tubes, launch them into Mars orbit, and transport them back to Earth. The ground element would be a high-containment facility that would isolate and protect the samples during initial sample characterization, which would include sample safety assessments and time-sensitive scientific investigations. These elements are currently in the planning and design stages of development, and represent an international effort of NASA, the European Space Agency (ESA), and many industry partners. The work presented here provides an overview of the planetary protection strategy of the third flight mission, the ESA-led Earth Return Orbiter, which hosts the NASA-provided Capture, Containment, and Return System. The orbiter would detect and capture the container with up to 30 sealed tubes previously put in Martian orbit, contain them in redundant containers to ensure that no potentially hazardous Mars particles are released, and return them to Earth through an entry vehicle. Both NASA and ESA policies comply with the United Nations’ Outer Space Treaty by planning to protect Earth’s biosphere from any potential adverse effects from material returned from solar system bodies beyond the Earth-Moon system. In the conduct of Mars Sample Return, the two agencies have mutually agreed to apply approaches consistent with their own planetary protection standards to the campaign elements they each provide.

Mars Sample Return↗

Towards a Robust Sampling Approach: A Computational Review and Design

As pointed out in several other works, the estimation of the reliability of the electrical grid can not be conducted without the estimation of the stochastic phenomena of electicity demand and electricity production by variable renewable sources. Therefore, sampling procedures have become integral in the design of engineering structures and analysis. Commonly referred to as Monte Carlo uncertainty analysis or integration, the general objective of these procedures is to establish specifics about the uncertainty of an output characteristic of such an engineering system, given uncertainty about its input characteristics. These sampling procedures are applied in a context in which establishing such specifics cannot be performed through other means. Variance-reduction techniques are designed to lessen the variability among estimators to estimate statistics of those output uncertainties. Importance-sampling techniques, on the other hand, are designed to reduce the number of samples needed to estimate a particular statistic—e.g., a tail probability. The combination of these approaches can reduce the computational burden considerably for a particular estimator and statistic. Importance sampling—geared and designed as it is toward improving a particular estimator—suffers, unfortunately, from the unintended consequence of reducing the performance of other estimators in terms of their variance. The objective of this paper is to offer an alternative sampling procedure where this variance does not grow unacceptably large for a suite of estimators. Moreover, it is anticipated that, with additional knowledge of how an engineering output characteristic responds to its input characteristics, tuning parameters of the input’s sampling procedure can be set to improve the output characteristic’s estimation. This work proves the effectiveness of the suggested alternative approach. Such positive outcome will lead to a decrease of the computational burden of performing stochastic optimization of integrated energy systems (e.g., components dispatch, and portfolio composition). In particular, capturing the contribution to the overall system cost of rare and unlikely events and patterns of the electricity demand and production will become less computationally expensive. This is due to the fact that the approach demonstrated here will allow the sampling of those rare occurrences more frequently without misrepresenting their probabilistic impacts.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Evaluating Manual Sampling Locations for Regulatory and Emergency Response

Drinking water systems commonly use manual or grab sampling to monitor water quality, identify or confirm issues, and verify that corrective or emergency response actions have been effective. In this paper, the effectiveness of regulatory sampling locations for emergency response is explored. An optimization formulation based on the literature was used to identify manual sampling locations to maximize overall nodal coverage of the system. Results showed that sampling locations could be effective in confirming incidents for which they were not designed. When evaluating sampling locations optimized for emergency response against regulatory scenarios, the average performance was reduced by 3%–4%, while using optimized regulatory sampling locations for emergency response reduced performance by 7%–10%. Secondary constraints were also included in the formulation to ensure geographical and water age diversity with minimal impact on the performance. This work highlighted that regulatory sampling locations provide value in responding to an emergency for these networks.

54 ENVIRONMENTAL SCIENCES↗

Internet of Samples

Abstract Material samples are indispensable data sources in many natural science, social science, and humanity disciplines. More and more researchers recognize that samples collected in one discipline can be of great value for another. This has motivated organizations that manage a large number of samples to make their holdings accessible to the world. Currently, multiple projects are working to connect natural history and other samples managed by individual institutions or individuals into a universe of samples that follow FAIR principles. This poster reports the progress of the US NSF‐funded iSamples project, in the context of other efforts initiated by US DOE, DiSCCo, BCoN, and GBIF. By October 2021, we will also be able to present an iSamples prototype. We encourage individual organizations that hold material samples to get to know these projects and help shape these projects to realize the goal of a global linked sample cloud that connects all material samples and is accessible to all.

Richard, Stephen M.↗

Comparison of measurement techniques and sorption of radium-226 in low and high salinity aqueous samples

Human activities have the potential to redistribute radium (Ra) in the marine environment in a manner that may necessitate monitoring or management of subsequent human or environmental exposures. There is therefore a need to identify accurate and accessible techniques for Ra measurement in high salinity samples and to describe the distribution of Ra in estuarine and marine environments, but most efforts in these areas have focused on low salinity matrices. In addition, rapid and reliable measurements are crucial for time-sensitive samples such as short-lived isotopes or emergency situations. The objective of this study is to describe the limits of detection, cost, and relative ease for measurement of Ra in both low and high salinity aqueous samples via three analytical methods: liquid scintillation counting (LSC), high purity germanium (HPGe) gamma spectrometry, and inductively coupled plasma mass spectrometry (ICP-MS). To contextualize these measurements for real-world scenarios, the partitioning of 226 Ra to substrates relevant to the marine environment was also characterized. Although HPGe detection with solid phase extraction had the lowest limit of detection for low salinity samples (0.27 Bq L −1 ), poor 226 Ra recovery for high salinity samples and high materials costs make this method prohibitive for many users. Limits of detection for high salinity samples were lower for LSC (1.28 Bq L −1 ) than for ICP-MS without dilution (11.4 Bq L −1 ), but significant and unexpected degradation of the high salinity LSC standards was observed after six months. Furthermore, our preferred measurement method for high salinity Ra samples is ICP-MS with sample dilution as necessary to reduce matrix effects.

07 ISOTOPE AND RADIATION SOURCES↗

Sample environment effects on synchrotron-measured temperature profiles in an approximant of optical floating zone crystal growth

Even though the growth of crystals using optical floating zone furnaces has had an immense scientific impact, the implementation of this method remains more of an art than a science due to the difficulty of obtaining quantitative information about the sample thermal profile during crystal growth. Building on recent work demonstrating that in situ synchrotron studies can be used to map sample rod temperatures during heating, investigations were carried out to better understand how the sample environment affects the sample temperature profile. Through a combination of experimental studies and modeling efforts, it is shown that the environment in the furnace can strongly influence the sample temperature at the lamp focus, the steepness of the vertical temperature gradient, and the timescale required for sample heating and cooling - effects which can combine to produce a strong history-dependence to sample temperature profiles. It is demonstrated that the furnace effects can be effectively captured in thermal models, allowing both the steady-state and time-dependent behavior of the sample to be accurately reproduced with predictive models, and providing a launching point for improved furnace designs that can more readily deliver desired thermal profiles.

36 MATERIALS SCIENCE↗

Automated Coupling of Nanodroplet Sample Preparation with Liquid Chromatography–Mass Spectrometry for High-Throughput Single-Cell Proteomics

Single-cell proteomics can provide critical biological insight into the cellular heterogeneity that is masked by bulk-scale analysis. Here, we have developed a nanoPOTS (nanodroplet processing in one pot for trace samples) platform and demonstrated its broad applicability for single-cell proteomics. However, because of nanoliter-scale sample volumes, the nanoPOTS platform is not compatible with automated LC-MS systems, which significantly limits sample throughput and robustness. To address this challenge, we have developed a nanoPOTS autosampler allowing fully automated sample injection from nanowells to LC-MS systems. We also developed a sample drying, extraction, and loading workflow to enable reproducible and reliable sample injection. The sequential analysis of 20 samples containing 10 ng tryptic peptides demonstrated high reproducibility with correlation coefficients of >0.995 between any two samples. The nanoPOTS autosampler can provide analysis throughput of 9.6, 16, and 24 single cells per day using 120, 60, and 30 min LC gradients, respectively. As a demonstration for single-cell proteomics, the autosampler was first applied to profiling protein expression in single MCF10A cells using a label-free approach. At a throughput of 24 single cells per day, an average of 256 proteins was identified from each cell and the number was increased to 731 when the Match Between Runs algorithm of MaxQuant was used. Using a multiplexed isobaric labeling approach (TMT-11plex), ~77 single cells could be analyzed per day. We analyzed 152 cells from three acute myeloid leukemia cell lines, resulting in a total of 2558 identified proteins with 1465 proteins quantifiable (70% valid values) across the 152 cells. These data showed quantitative single-cell proteomics can cluster cells to distinct groups and reveal functionally distinct differences.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Femtosecond Laser Desorption Postionization MS vs ToF-SIMS Imaging for Uncovering Biomarkers Buried in Geological Samples

The study of lipid molecular fossils by traditional biomarker analysis requires bulk sample crushing, followed by solvent extraction, and then the analysis of the extract by gas chromatography-mass spectrometry (GC–MS). This traditional analysis mixes all organic compounds in the sample regardless of their origins, with a loss of information on the spatial distribution of organic molecules within the sample. These shortcomings can be overcome using the chemical mapping of intact samples. Spectroscopic techniques such as UV fluorescence or Raman spectroscopy, laser ablation inductively coupled plasma mass spectrometry, and time-of-flight secondary ion mass spectrometry (ToF-SIMS) are among those elemental and molecular mapping techniques. This study employed femtosecond (fs) laser ablation combined with single-photon ionization, a method called fs-laser desorption postionization mass spectrometry (fs-LDPI-MS). In this work, a pulsed ~75 fs, 800 nm laser was used to ablate the geological sample, which was then photoionized after a few microseconds by a pulsed 7.9 eV vacuum ultraviolet laser. An organic carbon-rich geological sample was used for this study to map hydrocarbon biomarkers in sediments that were previously studied by GC–MS. The petrography of this sample was examined by optical and fluorescence microscopy. It is demonstrated here that fs-LDPI-MS combined with petrography for multimodal imaging can expose buried compounds within the sample via in situ layer removal. When used in conjunction with traditional organic geochemical analysis, this method has the potential to determine the spatial distribution of organic biomarkers in geological material. Finally, fs-LDPI-MS imaging data are compared with ToF-SIMS imaging that is commonly used for such studies.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Block Design with Common Reference Samples Enables Robust Large-Scale Label-Free Quantitative Proteome Profiling

Label-free quantitative proteomics has become an increasingly popular tool for profiling global protein abundances. However, one major weakness is the potential performance drift of the LC-MS platform over time, which in turn limits its utility for analyzing large-scale sample sets. To address this, in this work we introduce an experimental and data analysis scheme based on a block-design with common controls within each block for enabling LC-MS-based large-scale label-free quantification. In this scheme, a large number of samples (e.g., >100 samples) are analyzed in smaller, and more manageable blocks, minimizing instrument drift and variability within a block. Furthermore, each designated block also contains common controls (or reference) samples for normalization within and across blocks. We demonstrated the effectiveness of this method by profiling the proteome response of human macrophage THP-1 cells to 11 engineered nanomaterials (ENMs) at two different doses. A total of 116 samples were analyzed in six blocks, yielding an average coverage of 4500 proteins per sample. The data revealed consistent quantification of proteins across all six blocks, as shown by highly stable quantification of house-keeping proteins in all samples and high levels of quantification correlation among samples from different blocks. The data also demonstrated that label-free quantification is robust and accurate enough to quantify even very subtle abundance changes as well as large fold-changes without potential ratio compression as often encountered with isobaric labeling. Our streamlined workflow is easy to implement and can be readily adapted to other large cohort studies for reproducible label-free proteome quantification.

59 BASIC BIOLOGICAL SCIENCES↗

Direct and rapid uranium isotopic analysis of environmental sample swipes via microextraction – LS-APGD/Orbitrap mass spectrometry

Accurate and precise isotopic analysis of actinides collected on environmental sample swipes from within nuclear facilities is an important safeguarding measure for detecting undeclared materials and activities. Traditional isotope ratio (IR) analysis of actinides using a bulk digestion approach can be laborious and time-consuming. Recently, direct analysis of environmental swipe samples using a microextraction approach has been explored as an alternative to conventional bulk digestion methods. The present study further evaluates this approach for accurate and precise isotopic analysis of uranium (U) in cotton swipe samples using the liquid sampling-atmospheric pressure glow discharge (LS-APGD)/Orbitrap-FTMS Booster detection system. The instrumental parameters were optimized, and quantitative capabilities were demonstrated through excellent linearity (R2 = 0.99) across a deposited mass range of 1 to 500 ng, with a limit of detection of 90 pg (in a 2 × 4 mm region) for total U; well below typical U concentrations (ng–mg range) in environmental sample swipes. The 235U/238U ratio showed excellent accuracy (−0.3% relative difference from the certificate value), with ∼20× improvement over previous microextraction-Orbitrap methods. Precision (∼1% relative standard deviation) was also greatly enhanced by ∼5×. This study also demonstrates that figures of merit achieved with neat U solution were not highly degraded when various external interfering elements/matrices were introduced. Finally, the developed method was successfully applied to analyze practical swipe samples representative of laboratory and outdoor environments. The presented figures of merit and applicability to practical samples further validate the capability of the platform for rapid analysis (<7 minutes) of environmental swipe samples.

Shrestha, Suraj [Clemson University, SC]↗

In Vitro Thrombogenicity Testing of Biomaterials in a Dynamic Flow Loop: Effects of Length and Quantity of Test Samples

Abstract The results of in vitro dynamic thrombogenicity testing of biomaterials and medical devices can be significantly impacted by test conditions. To develop and standardize a robust dynamic in vitro thrombogenicity tool, the key test parameters need to be appropriately evaluated and optimized. We used a flow loop test system previously developed in our laboratory to investigate the effects of sample length and the number of samples per test loop on the thrombogenicity results. Porcine blood heparinized to a donor-specific target concentration was recirculated at room temperature through polyvinyl chloride (PVC) tubing loops containing test materials for 1 h at 200 mL/min. Four test materials (polytetrafluoroethylene (PTFE), latex, PVC, and silicone) with various thrombotic potentials in two sample lengths (12 and 18 cm) were examined. For the 12-cm long materials, two different test configurations (one and two samples per loop) were compared. Thrombogenicity was assessed through percent thrombus surface coverage, thrombus weight, and platelet count reduction in the blood. The test system was able to effectively differentiate the thrombogenicity profile of the materials (latex > silicone > PVC ≥ PTFE) at all test configurations. Increasing test sample length by 50% did not significantly impact the test results as both 12 and 18 cm sample lengths were shown to equally differentiate thrombotic potentials between the materials. The addition of a second test sample to each loop did not increase the test sensitivity and may produce confounding results, and thus a single test sample per loop is recommended.

Engineering↗

Hierarchical Gaussian Random Field Sampling for Multilevel Markov Chain Monte Carlo: Coupling Stochastic Partial Differential Equation and the Karhunen–Loève Decomposition

This work introduces structure preserving hierarchical decompositions for sampling Gaussian random fields (GRFs) within the context of multilevel Bayesian inference in high-dimensional space. Existing scalable hierarchical sampling methods, such as those based on stochastic partial differential equations (SPDEs), often reduce the dimensionality of the sample space at the cost of accuracy of inference. Other approaches, such that those based on Karhunen-Loève (KL) expansions, offer sample space dimensionality reduction but sacrifice GRF representation accuracy and ergodicity of the Markov chain Monte Carlo (MCMC) sampler and are computationally expensive for high-dimensional problems. The proposed method integrates the dimensionality reduction capabilities of KL expansions with the scalability of SPDE-based sampling, thereby providing a robust, unified framework for high-dimensional uncertainty quantification (UQ) that is scalable and accurate, preserves ergodicity, and offers dimensionality reduction of the sample space. The hierarchy in our multilevel algorithm is derived from the geometric multigrid hierarchy. By constructing a hierarchical decomposition that maintains the covariance structure across the levels in the hierarchy, the approach enables efficient coarse-to-fine sampling while ensuring that all samples are drawn from the desired distribution. The effectiveness of the proposed method is demonstrated on a benchmark subsurface flow problem, demonstrating its effectiveness in improving computational efficiency and statistical accuracy. Furthermore, our proposed technique is more efficient and accurate and displays better convergence properties than existing methods for high-dimensional Bayesian inference problems.

Gaussian random fields↗