Search NASA⌕ Search

SEARCH · Search NASA

Results for “trend 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 145 records · Page 8
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Research Trends and Applications of PMUs

This work is a survey of current trends in applications of PMUs. PMUs have the potential to solve major problems in the areas of power system estimation, protection, and stability. A variety of methods are being used for these purposes, including statistical techniques, mathematical transformations, probability, and AI. The results produced by the techniques reviewed in this work are promising, but there is work to be performed in the context of implementation and standardization. As the smart grid initiative continues to advance, the number of intelligent devices monitoring the power grid continues to increase. PMUs are at the center of this initiative, and as a result, each year more PMUs are deployed across the grid. Since their introduction, myriad solutions based on PMU-technology have been suggested. The high sampling rates and synchronized measurements provided by PMUs are expected to drive significant advancements across multiple fields, such as the protection, estimation, and control of the power grid. This work offers a review of contemporary research trends and applications of PMU technology. Most solutions presented in this work were published in the last five years, and techniques showing potential for significant impact are highlighted in greater detail. Being a relatively new technology, there are several issues that must be addressed before PMU-based solutions can be successfully implemented. This survey found that key areas where improvements are needed include the establishment of PMU-observability, data processing algorithms, the handling of heterogeneous sampling rates, and the minimization of the investment in infrastructure for PMU communication. Solutions based on Bayesian estimation, as well as those having a distributed architectures, show great promise. The material presented in this document is tailored to both new researchers entering this field and experienced researchers wishing to become acquainted with emerging trends.

42 ENGINEERING↗
DOE OSTI2021

Injection data analysis using material balance time for CO 2 storage capacity estimation in deep closed saline aquifers

Estimating the ultimate storage capacity of deep saline aquifers is important to address the formation potential to store the envisioned large volumes of CO 2 . Injection data (i.e. injection rate, bottomhole pressure, and cumulative injected volume of CO 2 ) are routinely recorded during storage operations. These data contain valuable information on the subsurface (e.g. the reservoir pore volume and the formation storage capacity) that can be extracted. In this paper, we present a two-step graphical technique to infer the pore volume and the ultimate storage capacity of closed saline aquifers by analyzing the available injection data. First, the pore volume is inferred through adapting the concept of the material balance time. Material balance time is an approximate superposition time function developed to interpret production data from oil and gas wells operating at variable pressure/rate conditions during the boundary-dominated flow period. Using material balance techniques, the ultimate storage capacity is then estimated through linear extrapolation of the average pressure trend to the maximum allowable pressure the formation can withstand. The average pressure is not available in practice, but is can be obtained from the injection data. Two approaches are presented in this study to calculate the average pressure; namely the rigorous and the approximate approaches. Unlike the rigorous approach, the approximate approach does not require a prior knowledge of some reservoir properties (e.g. relative permeability, absolute permeability, formation porosity and thickness) to calculate the average pressure. To investigate its potential and reliability in analyzing CO 2 injection data, the proposed technique is applied to four synthetic cases representing different well operating conditions. Results indicate that the approximate approach consistently overestimates the actual (simulated) storage capacity as compared to the rigorous approach. The agreement - between the inferred and the simulated reservoir pore volume, and between the analytical and numerical estimates of storage capacity - validates the potential application of the technique to CO 2 storage in closed saline aquifers. The technique is further substantiated through application to a field data set utilized from a commercial-scale geological storage (CGS) project. Finally, field data interpretation shows that the proposed technique can be utilized to identify the degree of hydraulic continuity and reservoir compartmentalization within a target formation by interpreting the corresponding pressure and rate responses.

02 PETROLEUM↗
DOE OSTI2026

Improved Spectral Classification of Local K/M Dwarfs in SDSS Surveys. I. Chemodynamic Validation and Metallicity Calibration

An examination of 109,276 spectra of low-mass stars in the Sloan Digital Sky Survey (SDSS) data archive, collected pre-2010, provides a broad collection of K/M (sub)dwarfs tracing the local (d ≲ 250 pc) population of the thin-disk, thick-disk, and halo, based on 3D kinematics. These populations have distinct metallicities and kinematics, which should be reflected in lower-mass members. One complication is measuring metallicities of M dwarfs from low-resolution spectroscopy or photometry alone remains a challenge, even with the availability of Gaia data. To better characterize physical parameters of low-mass stars observed by SDSS, we define a set of 536 improved, empirical K/M dwarf classification templates to more accurately determine spectral subtypes (K5.5–M8.5) and morphology classes (MCs) (0.5–12.5). We select the best-fit template to every archival SDSS spectrum in the range 5000–8000 Å based on minimum χ 2 value. We then use secondary metallicity estimators to calibrate the most likely [Fe/H] values corresponding to the assigned MC. We confirm that the most metal-rich K/M dwarfs have thin-disk kinematics, and we observe the known effect of Galactic radial metallicity migration, which we confirm is strongly imprinted in the chemodynamics of the local disk M dwarfs. Confirmation of these chemodynamic trends validates the precision of our improved templates for estimating [Fe/H] for M (sub)dwarfs, notably at the low-metallicity end (–3 < [Fe/H] < –1). Substructure in the chemodynamics of the most metal-poor K/M subdwarfs begins to distinguish between local Gaia-Enceladus and local in situ halo objects. This methodology opens the door to use the ubiquitous, local K/M (sub)dwarfs to trace local Galactic chemodynamics with the highest possible resolution.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

The mechanical behavior of metal-halide perovskites: Elasticity, plasticity, fracture, and creep

Wide ranging mechanical properties — elasticity, plasticity, fracture, and creep — most relevant to the mechanical reliability of perovskite solar cells (PSCs) are systematically investigated. High quality bulk single-crystals of the commonly studied metal-halide perovskites (MHPs) relevant to PSCs are fabricated and studied: CH 3 NH 3 PbBr 3 (MAPbBr 3 ) and CH 3 NH 3 PbI 3 (MAPbI 3 ). The first direct measurement of MHP Young's modulus (E) using uniaxial compression reveals E <100> of 13.1 ± 1.3 and 10.6 ± 1.0 GPa for MAPbBr 3 and MAPbI 3 , respectively. Further, the Vickers micro-hardness H (100) of MAPbBr 3 and MAPbI 3 is 0.54±0.02 GPa and 0.76±0.05 GPa, respectively. The Vickers micro-indentation fracture toughness K IC of MAPbBr 3 and MAPbI 3 is estimated at 0.20±0.03 MPa∙m 0.5 and 0.18±0.03 MPa∙m 0.5 , respectively. The stress-exponent, n, extracted from nanoindentation creep data is ~8 and ~10 for MAPbBr 3 and MAPbI 3 , respectively. The trends in these properties are discussed. These properties are best estimates and are recommended for use in future mechanical behavior and reliability analyses of MHPs and PSCs.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Constraining long‐term model predictions for woody growth using tropical tree rings

The strength and persistence of the tropical carbon sink hinges on the long-term responses of woody growth to climatic variations and increasing CO 2 . However, the sensitivity of tropical woody growth to these environmental changes is poorly understood, leading to large uncertainties in growth predictions. Here, we used tree ring records from a Southeast Asian tropical forest to constrain ED2.2-hydro, a terrestrial biosphere model with explicit vegetation demography. Specifically, we assessed individual-level woody growth responses to historical climate variability and increases in atmospheric CO 2 (C a ). When forced with historical C a , ED2.2-hydro reproduced the magnitude of increases in intercellular CO 2 concentration (a major determinant of photosynthesis) estimated from tree ring carbon isotope records. In contrast, simulated growth trends were considerably larger than those obtained from tree rings, suggesting that woody biomass production efficiency (WBPE = woody biomass production:gross primary productivity) was overestimated by the model. The estimated WBPE decline under increasing C a based on model-data discrepancy was comparable to or stronger than (depending on tree species and size) the observed WBPE changes from a multi-year mature-forest CO 2 fertilization experiment. In addition, we found that ED2.2-hydro generally overestimated climatic sensitivity of woody growth, especially for late-successional plant functional types. The model-data discrepancy in growth sensitivity to climate was likely caused by underestimating WBPE in hot and dry years due to commonly used model assumptions on carbon use efficiency and allocation. To our knowledge, this is the first study to constrain model predictions of individual tree-level growth sensitivity to C a and climate against tropical tree-ring data. Our results suggest that improving model processes related to WBPE is crucial to obtain better predictions of tropical forest responses to droughts and increasing C a . Finally, more accurate parameterization of WBPE will likely reduce the stimulation of woody growth by C a rise predicted by biosphere models.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Misclassification of primary liver cancer in the Life Span Study of atomic bomb survivors

Primary liver cancer is difficult to diagnose accurately at death, due to metastases from nearby organs and to concomitant diseases, such as chronic hepatitis and cirrhosis. Trends in diagnostic accuracy could affect radiation risk estimates for incident liver cancer by altering background rates or by impacting risk modification by sex and age. We quantified the potential impact of death‐certificate inaccuracies on radiation risk estimates for liver cancer in the Life Span Study of atomic bomb survivors. True‐positive and false‐negative rates were obtained from a previous study that compared death‐certificate causes of death with those based on pathological review, from 1958 to 1987. We assumed various scenarios for misclassification rates after 1987. We obtained estimated true positives and estimated false negatives by stratified sampling from binomial distributions with probabilities given by the true‐positive and false‐negative rates, respectively. Poisson regression methods were applied to highly stratified person‐year tables of corrected case counts and accrued person years. During the study period (1958–2009), there were 1,885 cases of liver cancer, which included 383 death‐certificate‐only (DCO) cases; 1,283 cases with chronic liver disease as the underlying cause of death; and 150 DCO cases of pancreatic cancer among 105,444 study participants. Across the range of scenarios considered, radiation risk estimates based on corrected case counts were attenuated, on average, by 13–30%. Our results indicated that radiation risk estimates for liver cancer were potentially sensitive to death‐certificate inaccuracies. Additional data are needed to inform misclassification rates in recent years.

French, Benjamin↗
DOE OSTI2022

Experimental study of the partitioning of some platinum group elements (Pd and Ir) between orthopyroxene and silicate melt

Past experiments and observations on natural samples have largely focused on the roles of olivine and chromite in controlling the behaviour of the platinum-group elements (PGE) during melting and solidification, whereas other phases, such as pyroxene, have gone largely uncharacterized. Here, to address this, experiments have been done to measure the partitioning of Pd (with a subset of results for Ir), between orthopyroxene and silicate melt at 1340 °C, 0.1 MPa and log fO 2 of FMQ - 1 to FMQ + 6 (FMQ = Fayalite-Magnetite-Quartz). The X-ray Absorption Near-Edge Structure (XANES) was measured in a subset of experiment glasses. Glass concentrations of Pd (corrected to unit Pd activity) increase from ~6 to ~650 ug/g over the fO 2 range of experiments. The slope of the solubility-fO 2 relation is consistent with Pd 1+ as the dominant oxidation state, with evidence for Pd 0 and Pd 2+ at the lowest and highest experiment fO 2 , respectively. Consistent with this result, the XANES reveal spectral features similar to Pd 0 and Pd 2+ spectral reference materials (specRM) at the most reduced and oxidized synthesis conditions, respectively. Other lines of evidence require the presence of a third melt species, here interpreted to be Pd 1+ . Values of orthopyroxene/melt partition coefficients for Pd (D Pd Opx/melt ) are 0.0051 (+/-0.006) at log fO 2 < ΔFMQ + 3, increasing with fO 2 to a maximum of 0.013 at ~FMQ + 6. Sodium partition coefficients, expected to be similar to Pd, range from 0.0061 (+/-0.00061) at FMQ + 3, increase to 0.007–0.009 at higher fO 2 , but with no clear systematic trend. A value for D Ir Opx-melt of ~0.6 was measured at ~FMQ + 4, indicating significantly more compatible behaviour for Ir relative to Pd. Partitioning results are interpreted in the context of the Blundy-Wood elastic strain model in which the variation in partitioning is related to ionic radius mismatch to an optimal crystallographic site size. Based on the trend in ionic radius with oxidation state, the estimated ionic radius of Pd 1+ in octahedral coordination is similar to Na 1+ , and comparison to previous orthopyroxene-melt partitioning experiments suggests D Pd1+ opx/melt and DNaopx/melt should be nearly identical, consistent with the results of this study. The ionic radius of VI-fold Pd 2+ is close to the optimal M2 site size, so an increased proportion of this species with fO 2 accounts for the larger values of D Pd opx/melt at the highest fO 2 investigated. The much larger partition coefficient for Ir is consistent with the presence of Ir 2+ , whose estimated ionic radius is close to Fe 2+ and Mg 2+ , as well as predictions for the optimal M1 site size. With the assumption that D Pd opx/melt = D Na opx/melt , combined with a revised value for the Pd content of the primitive mantle, a melting model is presented that better reproduces the Pd concentration of high degree melts from sulfide-free mantle sources.

58 GEOSCIENCES↗
DOE OSTI2026

Growth of deuterium supersaturated surface layer with increasing ion flux and fluence in plasma-exposed tungsten

Deuterium supersaturated surface layer (DSSL) in tungsten, a few nm thick layer exhibiting extremely high deuterium content (> 5%), has been studied as a function of deuterium plasma ion flux and fluence. Tungsten samples were exposed at 400 K and deuterium ion energy of ∼ 60 eV. The deuterium ion flux spanned over an order of magnitude, being 7.3 × 10 20 , 4.2 × 10 21 , and 3.8 × 10 22 D/m 2 . The fluence ranged from 3.4 × 1024 to 3.4 × 1025 D/m2. The samples and their deuterium content were analyzed by nuclear reaction analysis (NRA), scanning transmission electron microscopy (STEM), and thermal desorption spectroscopy (TDS). The largest thickness of DSSL was found to be around 10.5 nm and was observed in the case of the highest exposure flux and fluence, whereas the layer was only about 3.8 nm thick in the case of the lowest value of the two parameters. The thickness of the DSSL was found to monotonically increase with both deuterium flux and fluence. Finally, similar to the thickness, the estimated D concentration seems to follow the same trend, increasing from the lowest value of around 3 at.% to the highest value of around 5 at.%.

Deuterium↗
DOE OSTI2020

The Community Earth System Model version 2 (CESM2)

An overview of the Community Earth System Model version 2 (CESM2) is provided, including a discussion of the challenges encountered during its development and how they were addressed. In addition, an evaluation of a pair of CESM2 long pre-industrial control and historical ensemble simulations is presented. These simulations were performed using the nominal 1° horizontal resolution configuration of the coupled model with both the “low-top” (40 km, with limited chemistry) and “high-top” (140 km, with comprehensive chemistry) versions of the atmospheric component. CESM2 contains many substantial science and infrastructure improvements and new capabilities since its previous major release, CESM1, resulting in improved historical simulations in comparison to CESM1 and available observations. These include major reductions in low latitude precipitation and short-wave cloud forcing biases; better representation of the Madden-Julian Oscillation; better El Niño – Southern Oscillation-related teleconnections; and a global land carbon accumulation trend that agrees well with observationally-based estimates. Most tropospheric and surface features of the low- and high-top simulations are very similar to each other, so these improvements are present in both configurations. CESM2 has an equilibrium climate sensitivity of 5.1-5.3°C, larger than in CESM1, primarily due to a combination of relatively small changes to cloud microphysics and boundary layer parameters. In contrast, CESM2’s transient climate response of 1.9-2.0°C is comparable to that of CESM1. The model outputs from these and many other simulations are available to the research community, and theyrepresent CESM2’s contributions to the Coupled Model Intercomparison Project phase 6 (CMIP6).

Danabasoglu, Gokhan↗
DOE OSTI2023

Observed increase in the peak rain rates of monsoon depressions

Most extreme precipitation in the densely populated region of central India is produced by atmospheric vortices called monsoon lows and monsoon depressions. Here we use satellite and gauge-based precipitation estimates with atmospheric reanalyses to assess 40-year trends in the rain rates of these storms, which have remained unknown. We show that rain rates increased in the rainiest quadrant of monsoon depressions, southwest of the vortex center; precipitation decreased in eastern quadrants, yielding no clear trend in precipitation averaged over the entire storm diameter. In an atmospheric reanalysis, ascent increased in the region of amplifying precipitation, but we could not detect trends in the intensity of rotational winds around the storm center. These storm changes occurred in a background environment where humidity increased rapidly over land while warming was more muted. Monsoon lows, which we show produce less precipitation than depressions, exhibit weaker trends that are less statistically robust.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Results from the Alfvén Eigenmode Active Diagnostic during the 2019-2020 JET deuterium campaign

This paper presents results of extensive analysis of mode excitation observed during the operation of the Alfvén Eigenmode Active Diagnostic (AEAD) in the JET tokamak during the 2019-2020 deuterium campaign. Here, six of eight toroidally spaced antennas, each with independent power and phasing, were successful in actively exciting stable MHD modes in 479 plasmas. In total, magnetic resonances were detected with up to fourteen fast magnetic probes. In this work, we present the calculations of resonant frequencies $f_0$, damping rates γ < 0, and toroidal mode numbers $\textit{n}$, spanning the parameter ranges $f_0$ ≈ 30–250 kHz, –γ ≈ 0–13 kHz, and $\vert n \vert \leq {30}$. In general, good agreement is seen between the resonant and the calculated toroidal Alfvén Eigenmode frequencies, and between the toroidal mode numbers applied by the AEAD and estimated of the excited resonances. We note several trends in the database: the probability of resonance detection decreases with plasma current and external heating power; the normalized damping rate increases with edge safety factor but decreases with external heating. These results provide key information to prepare future experimental campaigns and to better understand the physics of excitation and damping of Alfvén Eigenmodes in the presence of alpha particles during the upcoming DT campaign, thereby extrapolating with confidence to future tokamaks.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Potential Impacts of Accelerated Climate Change (4th Annual Report)

Pacific Northwest National Laboratory (PNNL) has completed the fourth year draft Annual Letter Report that summarizes i) NRC needs for climate and hydrologic information, ii) recent research into climate and hydrologic changes that may affect nuclear power plants; iii) recent reports by federal agencies, as well as national and international organizations tasked with climate change research and developing policy with respect to climate change and mitigation; and iv) trends and advances in accommodating required climate change estimates in site specific evaluations. The Northeast U.S., a region defined by the Forth National Climate Assessment is the focus of this report. The NRC COR will review and provide comments on this draft within four weeks.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Giant kelp genetic monitoring before and after disturbance reveals stable genetic diversity in Southern California

Given the impacts of climate change and other anthropogenic stressors on marine systems, there is a need to accurately predict how species respond to changing environments and disturbance regimes. The use of genetic tools to monitor temporal trends in populations gives ecologists the ability to estimate changes in genetic diversity and effective population size that may be undetectable by traditional census methods. Although multiple studies have used temporal genetic analysis, they usually involve commercially important species, and rarely sample before and after disturbance. In this study, we run a temporal analysis of giant kelp, Macrocystis pyrifera , genetic diversity over the scope of 10 years (2008-2018) using the same microsatellite marker panel to assess the genetic consequences of disturbance in several populations of giant kelp ( Macrocystis pyrifera ) in the Southern California Bight. The study is a rare pre- and post-disturbance microsatellite analysis that included declines to giant kelp caused by the 2015/16 El Nino Southern Oscillation event. We used canopy biomass estimated by remote sensing (Landsat) to quantify the extent of disturbance to kelp beds, and sea surface temperature data to understand how kelp was pushed towards its temperature limits during this period. Despite prolonged periods with decreased canopy at several sites, no changes in genetic structure and allelic richness were observed. We argue that giant kelp in the region is best described as a “patchy population” system where true extinctions are rare. We discuss how deep refugia of subsurface sporophytes and cryptic microscopic life stages could have kept genetic diversity through disturbance. Given the increasing effects of climate change and uncertainty in modeling impacts of species with cryptic life history stages, we suggest further investigation to reveal the role such stages play in species resilience. Genetic monitoring studies of sites selected by remote census demographic and climate surveys should be continued in the future given the predicted impacts of climate change.

Klingbeil, III, William H.↗
DOE OSTI2021

Results from the Alfven Eigenmode Active Diagnostic during the 2019-2020 JET deuterium campaign

This paper presents results of extensive analysis of mode excitation observed during the operation of the Alfven Eigenmode Active Diagnostic (AEAD) in the JET tokamak during the 2019-2020 deuterium campaign. Six of eight toroidally spaced antennas, each with independent power and phasing, were successful in actively exciting stable MHD modes in 479 plasmas. In total, 4768 magnetic resonances were detected with up to fourteen fast magnetic probes. In this work, we present the calculations of resonant frequencies f0, damping rates \gamma < 0, and toroidal mode numbers n, spanning the parameter range f0 = 30 - 250 kHz, -\gamma = 0 - 13 kHz, and |n| < 30. In general, good agreement is seen between the resonant and the calculated toroidal Alfven Eigenmode frequencies, and between the toroidal mode numbers applied by the AEAD and estimated of the excited resonances. We note several trends in the database: the probability of resonance detection decreases with plasma current and external heating power; the normalized damping rate increases with edge safety factor but decreases with external heating. These results provide key information to prepare future experimental campaigns and to better understand the physics of excitation and damping of Alfven Eigenmodes in the presence of alpha particles during the upcoming DT campaign, thereby extrapolating with confidence to future tokamaks.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Changes in Climate and Its Effect on Timing of Snowmelt and Intensity-Duration-Frequency Curves

Snow is a critical water resource for much of the U.S. and failure to account for changes in climate could deleteriously impact military assets. In this study, we produced historical and future snow trends through modeling at three military sites (in Washington, Colorado, and North Dakota) and the Western U.S. For selected rivers, we performed seasonal trend analysis of discharge extremes. We calculated flood frequency curves and estimated the probability of occurrence of future annual maximum daily rainfall depths. Additionally, we generated intensity-duration-frequency curves (IDF) to find rainfall intensities at several return levels. Generally, our results showed a decreasing trend in historical and future snow duration, rain-on-snow events, and snowmelt runoff. This decreasing trend in snowpack could reduce water resources. A statistically significant increase in maximum streamflow for most rivers at the Washington and North Dakota sites occurred for several months of the year. In Colorado, only a few months indicated such an increase. Future IDF curves for Colorado and North Dakota indicated a slight increase in rainfall intensity whereas the Washington site had about a twofold increase. This increase in rainfall intensity could result in major flood events, demonstrating the importance of accounting for climate changes in infrastructure planning.

58 GEOSCIENCES↗
DOE OSTI2021

A Review and Outlook on Energy Consumption Estimation Models for Electric Vehicles

Electric vehicles (EVs) are critical to the transition to a low-carbon transportation system. The successful adoption of EVs heavily depends on energy consumption models that can accurately and reliably estimate electricity consumption. This paper reviews the state-of-the-art of EV energy consumption models, aiming to provide guidance for future development of EV applications. Here, we summarize influential variables of EV energy consumption into four categories: vehicle component, vehicle dynamics, traffic and environment related factors. We classify and discuss EV energy consumption models in terms of modeling scale (microscopic vs. macroscopic) and methodology (data-driven vs. rule-based). Our review shows trends of increasing macroscopic models that can be used to estimate trip-level EV energy consumption and increasing data-driven models that utilized machine learning technologies to estimate EV energy consumption based on large volume real-world data. We identify research gaps for EV energy consumption models, including the development of energy estimation models for modes other than personal vehicles (e.g., electric buses, electric trucks, and electric non-road vehicles); the development of energy estimation models that are suitable for applications related to vehicle-to-grid integration; and the development of multi-scale energy estimation models as a holistic modeling approach.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2022

Spatial and temporal prediction of radiation dose rates near Fukushima Daiichi Nuclear Power Plant

In this paper, we have developed a methodology to estimate the spatiotemporal distribution of radiation air dose rates around the Fukushima Daiichi Nuclear Power Plant (FDNPP). In our exploratory data analysis, we found that (1) the temporal evolution of dose rates is composed of a log-linear decay trend and fluctuations of air dose rates that are spatially correlated among adjacent monitoring posts; and (2) the slope of the log-linear environmental decay trend can be represented as a function of the apparent initial dose rates, coordinate position, land-use type, and soil type. From these observations, we first estimated the log-linear decay trend at each location based on these predictors, using the random forest method. We then developed a modified Kalman filter coupled with a Gaussian process model to estimate the dose-rate time series at a given location and time. We applied this method to the Fukushima evacuation zone (as of March 2017), which included 17 monitoring post locations (with monitoring datasets collected between 2014 and 2018) and generated a time series of dose-rate maps. Our results show that this approach allows us to produce accurate spatial and temporal predictions of radiation dose-rate maps using limited spatiotemporal measurements.

61 RADIATION PROTECTION AND DOSIMETRY↗
DOE OSTI2021

Optimal interpolation of global dissolved oxygen: 1965–2015

Oxygen inventory of the global ocean has declined in recent decades potentially due to the warming-induced reduction in solubility as well as the circulation and biogeochemical changes associated with ocean warming and increasing stratification. Earth System Models predict continued oxygen decline for this century with profound impacts on marine ecosystem and fisheries. Observational constraint on the rate of oxygen loss is crucial for assessing the ability of models to accurately simulate these changes. There are only a few observational assessments of the global oceanic oxygen inventory reporting a range of oxygen loss. This study develops a gridded data set of dissolved oxygen for the global oceans using optimal interpolation method. The resulting gridded product includes full-depth map of dissolved oxygen as 5-year moving average from 1965 to 2015 with uncertainty estimates. The uncertainty can come from unresolved small-scale and high-frequency variability and mapping errors. The multi-decadal trend of global dissolved oxygen is in the range of -281 to -373 Tmol/decade. This estimate is more conservative than previous works. In this study, the grid points far from the observations are essentially set equal to zero anomaly from the climatology. Calculating global inventory with this approach produces a relatively conservative estimate; thus, the results from this study likely provide a useful lower bound estimate of the global oxygen loss.

54 ENVIRONMENTAL SCIENCES↗