Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “yield estimation”

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 199 records · Page 11

An overview of wind-energy-production prediction bias, losses, and uncertainties

Abstract. The financing of a wind farm directly relates to the preconstruction energy yield assessments which estimate the annual energy production for the farm. The accuracy and the precision of the preconstruction energy estimates can dictate the profitability of the wind project. Historically, the wind industry tended to overpredict the annual energy production of wind farms. Experts have been dedicated to eliminating such prediction errors in the past decade, and recently the reported average energy prediction bias is declining. Herein, we present a literature review of the energy yield assessment errors across the global wind energy industry. We identify a long-term trend of reduction in the overprediction bias, whereas the uncertainty associated with the prediction error is prominent. We also summarize the recent advancements of the wind resource assessment process that justify the bias reduction, including improvements in modeling and measurement techniques. Additionally, because the energy losses and uncertainties substantially influence the prediction error, we document and examine the estimated and observed loss and uncertainty values from the literature, according to the proposed framework in the International Electrotechnical Commission 61400-15 wind resource assessment standard. From our findings, we highlight opportunities for the industry to move forward, such as the validation and reduction of prediction uncertainty and the prevention of energy losses caused by wake effect and environmental events. Overall, this study provides a summary of how the wind energy industry has been quantifying and reducing prediction errors, energy losses, and production uncertainties. Finally, for this work to be as reproducible as possible, we include all of the data used in the analysis in appendices to the article.

17 WIND ENERGY↗

Review of Neutron Yield from (α, n) Reactions: Data, Methods, and Prospects

Understanding the radiogenic neutron production rate through the (α, n) reaction is crucial in many areas of physics, including dark matter searches, neutrino studies, and nuclear astrophysics. In addition to its relevance for fundamental research, the (α, n) reaction also plays a significant role in nuclear energy technologies, for example by contributing to neutron production in subcritical systems using UO2, as well as in applications such as medical physics. This review examines the current state of (α, n) yield calculations and neutron spectra, describes the computational tools used for their estimation, and discusses the available cross-section data. We discuss the uncertainties affecting (α, n) yield estimations and propose a strategy to enhance their accuracy. Furthermore, this paper discusses and emphasizes the need for new measurements of (α, n) cross-sections for a variety of relevant materials. Such measurements are essential for improving neutron flux predictions, which are crucial for reducing uncertainties in sensitivity estimates for next-generation physics experiments operating in the keV—MeV range.

(α, n)reactions↗

Techno-economic life cycle assessment of CO 2 -EOR operations towards net negative emissions at farnsworth field unit

Optimizations of CO 2 Water Alternating Gas(WAG)- systems with multi-objectives of incremental recovery and maximization of CO 2 storage are challenging. Here, the incorporation of a total Greenhouse gas (GHG) life cycle assessment is mostly ignored leading to inaccurate estimation of overall net carbon emissions of their operations. In this study, the effect of a total GHG life cycle assessment on a multi-objective CO 2 -WAG optimization with integrated techno-economic assessment (TEA) which factors carbon tax credit is conducted. A life cycle assessment (LCA) was conducted utilizing a 20 -year optimized post history matched data from a high fidelity reservoir simulation model. Using data generated from the optimum result, a techno-economic life cycle analysis was further conducted. The first scenario classified as the base model had an estimated 81% of purchased CO 2 sequestered. The results through a comprehensive techno-economic LCA model yielded a net estimate of 73% of purchased CO 2 . The optimized forecasted model which considered key operational and reservoir factors such as WAG ratio, injection rates and periods, and well specification resulted in an improved sequestration of 92% of purchased CO 2 . However, this also dropped to 84% after taking it through LCA. These results clearly indicate a significant amount of net CO 2 is not accounted for when operations are not analyzed through LCA. From the LCA, direct flaring volumes of CO 2 , energy consumption and efficiency of unit equipment were noticed to be the major causes of these reductions. Considering ten main sources of energy as source of energy generation, a comparative techno-eco LCA was conducted. The results confirmed a lower net volume and NPV for energy sources with higher carbon footprints and vice versa. Thus, a total LCA of CO 2 -WAG greatly influences net storage factor of purchased CO 2 and hence project NPV where tax credit/incentives per ton of CO 2 sequestered is considered. Although operational conditions are optimized for best results, there are significant factors that leads to minimization of net storage factor. This study therefore provides an insightful information for optimizing CO 2 -WAG multi-objectives to achieve minimum GHG emission.

02 PETROLEUM↗

Harmonising the land-use flux estimates of global models and national inventories for 2000–2020

As the focus of climate policy shifts from pledges to implementation, there is a growing need to track progress on climate change mitigation at the country level, particularly for the land-use sector. Despite new tools and models providing unprecedented monitoring opportunities, striking differences remain in estimations of anthropogenic land-use CO 2 fluxes between, on the one hand, the national greenhouse gas inventories (NGHGIs) used to assess compliance with national climate targets under the Paris Agreement and, on the other hand, the Global Carbon Budget and Intergovernmental Panel on Climate Change (IPCC) assessment reports, both based on global bookkeeping models (BMs). Recent studies have shown that these differences are mainly due to inconsistent definitions of anthropogenic CO 2 fluxes in managed forests. Countries assume larger areas of forest to be managed than BMs do, due to a broader definition of managed land in NGHGIs. Additionally, the fraction of the land sink caused by indirect effects of human-induced environmental change (e.g. fertilisation effect on vegetation growth due to increased atmospheric CO 2 concentration) on managed lands is treated as non-anthropogenic by BMs but as anthropogenic in most NGHGIs. We implement an approach that adds the CO 2 sink caused by environmental change in countries' managed forests (estimated by 16 dynamic global vegetation models, DGVMs) to the land-use fluxes from three BMs. This sum is conceptually more comparable to NGHGIs and is thus expected to be quantitatively more similar. Our analysis uses updated and more comprehensive data from NGHGIs than previous studies and provides model results at a greater level of disaggregation in terms of regions, countries and land categories (i.e. forest land, deforestation, organic soils, other land uses). Our results confirm a large difference (6.7 GtCO 2 yr —1 ) in global land-use CO 2 fluxes between the ensemble mean of the BMs, which estimate a source of 4.8 GtCO 2 yr —1 for the period 2000–2020, and NGHGIs, which estimate a sink of —1.9 GtCO 2 yr —1 in the same period. Most of the gap is found on forest land (3.5 GtCO 2 yr —1 ), with differences also for deforestation (2.4 GtCO 2 yr —1 ), for fluxes from other land uses (1.0 GtCO 2 yr —1 ) and to a lesser extent for fluxes from organic soils (0.2 GtCO 2 yr —1 ). By adding the DGVM ensemble mean sink arising from environmental change in managed forests (—6.4 GtCO 2 yr —1 ) to BM estimates, the gap between BMs and NGHGIs becomes substantially smaller both globally (residual gap: 0.3 GtCO 2 yr —1 ) and in most regions and countries. However, some discrepancies remain and deserve further investigation. For example, the BMs generally provide higher emissions from deforestation than NGHGIs and, when adjusted with the sink in managed forests estimated by DGVMs, yield a sink that is often greater than NGHGIs. In summary, this study provides a blueprint for harmonising the estimations of anthropogenic land-use fluxes, allowing for detailed comparisons between global models and national inventories at global, regional and country levels. This is crucial to increase confidence in land-use emissions estimates, support investments in land-based mitigation strategies and assess the countries' collective progress under the Global Stocktake of the Paris Agreement.

54 ENVIRONMENTAL SCIENCES↗

On the uncertainty of long-period return values of extreme daily precipitation

Methods for calculating return values of extreme precipitation and their uncertainty are compared using daily precipitation rates over the Western U.S. and Southwestern Canada from a large ensemble of climate model simulations. The roles of return-value estimation procedures and sample size in uncertainty are evaluated for various return periods. We compare two different generalized extreme value (GEV) parameter estimation techniques, namely L-moments and maximum likelihood (MLE), as well as empirical techniques. Even for very large datasets, confidence intervals calculated using GEV techniques are narrower than those calculated using empirical methods. Furthermore, the more efficient L-moments parameter estimation techniques result in narrower confidence intervals than MLE parameter estimation techniques at small sample sizes, but similar best estimates. It should be noted that we do not claim that either parameter fitting technique is better calibrated than the other to estimate long period return values. While a non-stationary MLE methodology is readily available to estimate GEV parameters, it is not for the L-moments method. Comparison of uncertainty quantification methods are found to yield significantly different estimates for small sample sizes but converge to similar results as sample size increases. Finally, practical recommendations about the length and size of climate model ensemble simulations and the choice of statistical methods to robustly estimate long period return values of extreme daily precipitation statistics and quantify their uncertainty.

54 ENVIRONMENTAL SCIENCES↗

Searching for Neutrino Tridents in the NOvA Near Detector

This dissertation presents a search for neutrino trident production in the NOvA near detector through the coherent ``dimuon" channel: $\nu_\mu +\hspace{1pt}\text{X} \rightarrow \nu_\mu + \mu^- + \mu^+ +\hspace{1pt}\text{X}$. Trident production is a rare, purely electroweak process with sensitivity to physics beyond the Standard Model. The theoretical background, motivation for studying the process, and previous experimental measurements are reviewed. The analysis uses data collected by the NOvA near detector (ND) from Fermilab's Neutrinos at the Main Injector (NuMI) beam between November 2014 and February 2024, corresponding to an exposure of $25.5\times 10^{20}$ protons on target. The ND is a segmented tracking calorimeter located 800~m from the beam target, receiving neutrinos with a mean energy of 2~GeV. A multi-pass background reduction strategy is implemented, including the development of a novel dimuon-specific tracking technique. Trident candidates are identified using a boost ed decision tree classifier trained on simulated signal and background events. Limited background Monte Carlo statistics necessitate the use of functional fits to sideband data, which are extrapolated to estimate backgrounds in the signal region. The unblinded data contain 9 trident-like events, with an estimated background of 5.66 $\pm$ 5.15 events. This yields a best fit estimate of 3.34 tridents compared to the Standard Model prediction of 4.66. A profiled Feldman-Cousins method is used to determine a 90\% confidence interval of [0,9.1] on the number of signal events, corresponding to an upper limit of 1.95$\times$ the Standard Model prediction. This result represents the lowest energy search for trident events to date, and the first experimental contribution to the process in 27 years.

Bowles, Reed Scott [Indiana U.]↗

Welch Method and Bootstrapping Applied to Subcritical Gamma Noise

We measured the prompt neutron decay constant 𝛼 of the CROCUS zero-power reactor at the Swiss Federal Institute of Technology Lausanne using cross-power spectral density (CPSD) analysis of gamma-gamma correlations from two trans-stilbene organic scintillators positioned near the reactor core. We measured critical and subcritical states, with water levels ranging from 960 mm (critical) to 800 mm (𝜌=−1.4 $ subcritical). Our analysis used the Welch method, dividing signal segments for fast Fourier transform (FFT) frequency analysis and applying bootstrapping uncertainty quantification that uses Welch-defined segments. Results demonstrated a clear increase in the measured 𝛼 as reactor reactivity decreased, distinguishing critical from subcritical conditions. At the 960-mm critical level, 𝛼 was estimated at 155.9 ± 0.7 s −1 , and for the 800-mm subcritical level, 𝛼 increased significantly to 367.3 ± 6.9 s –1 . A linear regression of subcritical states yielded a critical estimate of 154.0 ± 3.1 s –1 , aligning with the static 𝛼 estimate at critical. The bootstrapping method produced normally distributed 𝛼 estimates, confirming data consistency. The gamma CPSD 𝛼 estimates clearly distinguish reactor states and improve monitoring of zero-power reactors. The future deployment of modular and microreactors as potential candidates for noise analysis is demonstrated in CROCUS, particularly zero-power mock-ups of new designs. The improvement of noise analysis in the subcritical domain from this work will support experimental data for reactor deployment and procedure.

CROCUS↗

Computer Model Calibration with Time Series Data Using Deep Learning and Quantile Regression

Computer models play a key role in many scientific and engineering problems. One major source of uncertainty in computer model experiments is input parameter uncertainty. Computer model calibration is a formal statistical procedure to infer input parameters by combining information from model runs and observational data. The existing standard calibration framework suffers from inferential issues when the model output and observational data are high-dimensional dependent data, such as large time series, due to the difficulty in building an emulator and the nonidentifiability between effects from input parameters and data-model discrepancy. To overcome these challenges, we, in this work, propose a new calibration framework based on a deep neural network (DNN) with long short-term memory layers that directly emulates the inverse relationship between the model output and input parameters. Adopting the “learning with noise” idea, we train our DNN model to filter out the effects from data-model discrepancy on input parameter inference. We also formulate a new way to construct interval predictions for DNN using quantile regression to quantify the uncertainty in input parameter estimates. Through a simulation study and real data application with the Weather Research and Forecasting Model Hydrological modeling system (WRF-Hydro), we show our approach can yield accurate point estimates and well-calibrated interval estimates for input parameters.

97 MATHEMATICS AND COMPUTING↗

Notes on Regression Analysis for Radar Parameter Estimation

A fundamental task of radar, beyond merely detecting a target, is to estimate some parameters associated with it. For example, this might include range, direction, velocity, etc. In any case, multiple measurements, often noisy, need to be processed to yield a ‘best estimate’ of the parameter. A common mathematical method for doing so is called “Regression” analysis. The goal is to minimize the expected squared error in the estimate. Even when alternate algorithms are considered, the least squared-error regression analysis is the benchmark against which alternatives are compared.

47 OTHER INSTRUMENTATION↗

High deuteron and neutron yields from the interaction of a petawatt laser with a cryogenic deuterium jet

A compact high-flux, short-pulse neutron source would have applications from nuclear astrophysics to cancer therapy. Laser-driven neutron sources can achieve fluxes much higher than spallation and reactor neutron sources by reducing the volume and time in which the neutron-producing reactions occur by orders of magnitude. We report progress towards an efficient laser-driven neutron source in experiments with a cryogenic deuterium jet on the Texas Petawatt laser. Neutrons were produced both by laser-accelerated multi-MeV deuterons colliding with Be and mixed metallic catchers and by d (d,n) 3 He fusion reactions within the jet. We observed deuteron yields of 10 13 /shot in quasi-Maxwellian distributions carrying ~8-10% of the input laser energy. We obtained neutron yields greater than 10 10 /shot and found indications of a deuteron-deuteron fusion neutron source with high peak flux (>10 22 cm -2 s -1 ). The estimated fusion neutron yield in our experiment is one order of magnitude higher than any previous laser-induced dd fusion reaction. Though many technical challenges will have to be overcome to convert this proof-of-principle experiment into a consistent ultra-high flux neutron source, the neutron fluxes achieved here suggest laser-driven neutron sources can support laboratory study of the rapid neutron-capture process, which is otherwise thought to occur only in astrophysical sites such as core-collapse supernova, and binary neutron star mergers.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Blowing Snow at McMurdo Station, Antarctica During the AWARE Field Campaign: Multi-Instrument Observations of Blowing Snow

In polar regions, blowing snow (BLSN) can play a substantial role in regional thermodynamics and radiation properties. The extent of the impact depends upon the depth of the BLSN layer, a property that is difficult to measure and not well understood. The West Antarctic Radiation Experiment (AWARE) Field Campaign saw the deployment of a large suite of observational instruments to McMurdo Station, Antarctica, allowing for in-depth investigation of BLSN. During the year of the campaign, BLSN occurred ~7.4% of the time. In this study, additional remote sensing observations are used to supplement an existing ceilometer-based BLSN characterization algorithm creating a multi-instrument (MI) depth estimation algorithm to yield layer depths with more certainty. Backscatter coefficient and linear depolarization ratios from the micropulse and high spectral resolution lidar (HSRL) are incorporated along with color ratio derived from HSRL backscatter and Ka-band radar reflectivity observations. This provided estimates of BLSN depth that were less impacted by noise/artifacts in any one set of observations. When applied to the data from AWARE, the MI algorithm yielded an average depth of 168.2 m with ~60% of plumes less than 200 m in depth. Events occurring with precipitation saw the most value in the additional instruments, with the suite of observations allowing for more separation of BLSN particles from fallstreaks.

54 ENVIRONMENTAL SCIENCES↗

Developing “inverted-corona” fusion targets as high-fluence neutron sources

We present experimental studies of inverted-corona targets as neutron sources at the OMEGA Laser Facility and the National Ignition Facility (NIF). Laser beams are directed onto the inner walls of a capsule via laser-entrance holes (LEHs), heating the target interior to fusion conditions. The fusion fuel is provided either as a wall liner, e.g., deuterated plastic (CD), or as a gas fill, e.g., D 2 gas. Such targets are robust to low-mode drive asymmetries, allowing for single-sided laser drive. On OMEGA, 1.8-mm-diameter targets with either a 10-μm CD liner or up to 2 atm of D 2 -gas fill were driven with up to 18 kJ of laser energy in a 1-ns square pulse. Neutron yields of up to 1.5 × 1010 generally followed expected trends with fill pressure or laser energy, although the data imply some mix of the CH wall into the fusion fuel for either design. Comparable performance was observed with single-sided (1x LEH) or double-sided (2x LEH) drive. NIF experiments tested the platform at scaled up dimensions and energies, combining a 15-μm CD liner and a 3-atm D2-gas fill in a 4.5-mm diameter target, laser-driven with up to 330 kJ. Neutron yields up to 2.6 × 10 12 were measured, exceeding the scaled yield expectation from the OMEGA data. The observed energy scaling on the NIF implies that the neutron production is gas dominated, suggesting a performance boost from using deuterium–tritium (DT) gas. In this work, we estimate that neutron yields exceeding 10 14 should be readily achievable using a modest laser drive of ~300 kJ with a DT fill.

47 OTHER INSTRUMENTATION↗

Constraining the masses of high-redshift clusters with weak lensing: Revised shape calibration testing for the impact of stronger shears and increased blending

Weak lensing measurements suffer from well-known shear estimation biases, which can be partially corrected for with the use of image simulations. Here we present an analysis of simulated images that mimic Hubble Space Telescope/Advance Camera for Surveys observations of high-redshift galaxy clusters, including cluster specific issues such as non-weak shear and increased blending. Our synthetic galaxies have been generated to have similar observed properties as the background-selected source samples studied in the real images. First, we used simulations with galaxies placed on a grid to determine a revised signal-to-noise-dependent (S/N KSB ) correction for multiplicative shear measurement bias, and to quantify the sensitivity of our KSB+ bias calibration to mismatches of galaxy or PSF properties between the real data and the simulations. Next, we studied the impact of increased blending and light contamination from cluster and foreground galaxies, finding it to be negligible for high-redshift (z > 0.7) clusters, whereas shear measurements can be affected at the ~1% level for lower redshift clusters given their brighter member galaxies. Finally, we studied the impact of fainter neighbours and selection bias using a set of simulated images that mimic the positions and magnitudes of galaxies in Cosmic Assembly Near-IR Deep Extragalactic Legacy Survey (CANDELS) data, thereby including realistic clustering. While the initial SExtractor object detection causes a multiplicative shear selection bias of –0.028 ± 0.002, this is reduced to –0.016 ± 0.002 by further cuts applied in our pipeline. Given the limited depth of the CANDELS data, we compared our CANDELS-based estimate for the impact of faint neighbours on the multiplicative shear measurement bias to a grid-based analysis, to which we added clustered galaxies to even fainter magnitudes based on Hubble Ultra Deep Field data, yielding a refined estimate of ~ –0.013. Our sensitivity analysis suggests that our pipeline is calibrated to an accuracy of ~0.015 once all corrections are applied, which is fully sufficient for current and near-future weak lensing studies of high-redshift clusters. As an application, we used it for a refined analysis of three highly relaxed clusters from the South Pole Telescope Sunyaev-Zeldovich survey, where we now included measurements down to the cluster core (r > 200 kpc) as enabled by our work. Compared to previously employed scales (r > 500 kpc), this tightens the cluster mass constraints by a factor 1.38 on average.

79 ASTRONOMY AND ASTROPHYSICS↗

Computational investigation of incident ion angles and material erosion at rough graphite and silicon carbide divertor surfaces

We present a computational investigation of the dependence of material erosion on the incident ion angle at rough graphite and silicon carbide divertor surfaces. Ion angle distributions (IADs) for D plasmas at NSTX-U and DIII-D divertors were calculated by an equation-of-motion model that traces the ion trajectories in the sheath. Then, the effective sputtering yields and ion shadowed area fractions were calculated by a Monte Carlo micro-patterning and roughness code that applied the calculated IADs to surface topographic data that were obtained from optical confocal microscopy of rough graphite and SiC divertor surfaces from NSTX-U and DIII-D experiments. The calculations found that the effective sputtering yields, the sputtering pattern, and the shadowed area are determined by the detailed surface topology rather than the root mean square roughness R RMS , which represents deviations from a flat surface. The suppression of the effective sputtering yields for rough surfaces compared to the yield for a smooth surface was accounted for by the change of the mean local incident ion angle (LIIA) (θ'). The mean surface inclination angle distribution (SIAD) (δ) was found to be a useful parameter to estimate the LIIA from the calculated IADs. We report global empirical formulas for the mean LIIA and fraction of the area shadowed from the main ions for D plasmas for rough surfaces with B-field incident angles α = 85°–89° as a function of the mean SIAD (δ). We propose the use of the mean LIIA (θ') to estimate the sputtering yield for rough surfaces from the angular dependence of the sputtering yield.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

An Overview of Wind Energy Production Prediction Bias, Losses, and Uncertainties

The financing of a wind farm directly relates to the preconstruction energy yield assessments which estimate the annual energy production for the farm. The accuracy and the precision of the preconstruction energy estimates can dictate the profitability of the wind project. Historically, the wind industry tended to overpredict the annual energy production of wind farms. Experts have been dedicated to eliminating such prediction errors in the past decade, and recently the industry is recording near-zero average energy prediction bias. Herein, we present an overview of the energy yield assessment errors across the global wind energy industry. We identify a long-term trend of reduction in the overprediction bias, whereas the uncertainty associated with the prediction error is prominent. We also summarize the recent advancements of the wind resource assessment process that justify the bias reduction, including the improvements in modeling and measurement techniques. Additionally, because the energy losses and uncertainties substantially influence the prediction error, we document and examine the estimated and observed loss and uncertainty values from the literature, according to the proposed framework in the International Electrotechnical Commission 61400-15 wind resource assessment standard. From our findings, we highlight the opportunities for the industry to move forward, such as the validation and reduction of prediction uncertainty, and the prevention of energy losses caused by wake effect and environmental events. Overall, this study provides a summary on how the wind energy industry has been quantifying and reducing prediction errors, energy losses, and production uncertainties. Finally, for this work to be as reproducible as possible, we include all of the data used in the analysis in appendices to the manuscript.

literature review↗

A Bayesian nonparametric analysis for zero-inflated multivariate count data with application to microbiome study

High-throughput sequencing technology has enabled researchers to profile microbial communities from a variety of environments, but analysis of multivariate taxon count data remains challenging. Here, we develop a Bayesian nonparametric (BNP) regression model with zero inflation to analyse multivariate count data from microbiome studies. A BNP approach flexibly models microbial associations with covariates, such as environmental factors and clinical characteristics. The model produces estimates for probability distributions which relate microbial diversity and differential abundance to covariates, and facilitates community comparisons beyond those provided by simple statistical tests. We compare the model to simpler models and popular alternatives in simulation studies, showing, in addition to these additional community-level insights, it yields superior parameter estimates and model fit in various settings. The model's utility is demonstrated by applying it to a chronic wound microbiome data set and a Human Microbiome Project data set, where it is used to compare microbial communities present in different environments.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Microbial communities in the liver and brain are informative for postmortem submersion interval estimation in the late phase of decomposition: A study in mouse cadavers recovered from freshwater

Introduction Bodies recovered from water, especially in the late phase of decomposition, pose difficulties to the investigating authorities. Various methods have been proposed for postmortem submersion interval (PMSI) estimation and drowning identification, but some limitations remain. Many recent studies have proved the value of microbiota succession in viscera for postmortem interval estimation. Nevertheless, the visceral microbiota succession and its application for PMSI estimation and drowning identification require further investigation. Methods In the current study, mouse drowning and CO 2 asphyxia models were developed, and cadavers were immersed in freshwater for 0 to 14 days. Microbial communities in the liver and brain were characterized via 16S rDNA high-throughput sequencing. Results Only livers and brains collected from 5 to 14 days postmortem were qualified for sequencing. There was significant variation between microbiota from liver and brain. Differences in microbiota between the cadavers of mice that had drowned and those only subjected to postmortem submersion decreased over the PMSI. Significant successions in microbial communities were observed among the different subgroups within the late phase of the PMSI in livers and brains. Eighteen taxa in the liver which were mainly related to Clostridium_sensu_stricto and Aeromonas , and 26 taxa in the brain which were mainly belonged to Clostridium_sensu_stricto , Acetobacteroides , and Limnochorda , were selected as potential biomarkers for PMSI estimation based on a random forest algorithm. The PMSI estimation models established yielded accurate prediction results with mean absolute errors ± the standard error of 1.282 ± 0.189 d for the liver and 0.989 ± 0.237 d for the brain. Conclusions The present study provides novel information on visceral postmortem microbiota succession in corpses submerged in freshwater which sheds new light on PMSI estimation based on the liver and brain in forensic practice.

Wang, Linlin↗