Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “ASSIMILATION”

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 1,531 records · Page 85
All sourcesNASA NTRSDOE OSTI
NASA NTRS2017

From LIMS to OMPS-LP: Limb Ozone Observations for Future Reanalyses

High vertical resolution and accuracy of ozone data from satellite-borne limb sounders has made them an invaluable tool in scientific studies of the middle and upper atmosphere. However, it was not until recently that these measurements were successfully incorporated in atmospheric reanalyses: of the major multidecadal reanalyses only ECMWF's (European Centre for Medium-Range Weather Forecasts') ERA (ECMWF Re-Analysis)-Interim/ERA5 and NASA's MERRA-2 (Modern-Era Retrospective Analysis for Research and Applications-2) use limb ozone data. Validation and comparison studies have demonstrated that the addition of observations from the Microwave Limb Sounder (MLS) on EOS (Earth Observing System) Aura greatly improved the quality of ozone fields in MERRA-2 making these assimilated data sets useful for scientific research. In this presentation, we will show the results of test experiments assimilating retrieved ozone from the Limb Infrared Monitor of the Stratosphere (LIMS, 1978/1979) and Ozone Mapping Profiler Suite Limb Profiler (OMPS-LP, 2012 to present). Our approach builds on the established assimilation methodology used for MLS in MERRA-2 and, in the case of OMPS-LP, extends the excellent record of MLS ozone assimilation into the post-EOS era in Earth observations. We will show case studies, discuss comparisons of the new experiments with MERRA-2, strategies for bias correction and the potential for combined assimilation of multiple limb ozone data types in future reanalyses for studies of multidecadal stratospheric ozone changes including trends.

Wargan, K.↗
NASA NTRS2017

Impact of Archeomagnetic Field Model Data on Modern Era Geomagnetic Forecasts

A series of geomagnetic data assimilation experiments have been carried out to demonstrate the impact of assimilating archeomagnetic data via the CALS3k.4 geomagnetic field model from the period between 10 and 1590 CE. The assimilation continues with the gufm1 model from 1590 to 1990 and CM4 model from 1990 to 2000 as observations, and comparisons between these models and the geomagnetic forecasts are used to determine an optimal maximum degree for the archeomagnetic observations, and to independently estimate errors for these observations. These are compared with an assimilation experiment that uses the uncertainties provided with CALS3k.4. Optimal 20 year forecasts in 1990 are found when the Gauss coefficients up to degree 3 are assimilated. In addition we demonstrate how a forecast and observation bias correction scheme could be used to reduce bias in modern era forecasts. Initial experiments show that this approach can reduce modern era forecast biases by as much as 50 percent.

Tangborn, Andrew↗
NASA NTRS2018

Progress Towards Integrating the Finite-Volume Cubed-Sphere (FV3) Dynamical Core Tangent Linear and Adjoint Models into JEDI

The Joint Effort for Data assimilation Integration (JEDI) -- led by the Joint Center for Satellite Data Assimilation (JCSDA) -- is an inter-organizational endeavor to develop a common framework for performing data assimilation. This extensive framework will ultimately provide solvers, observation operators, interpolation and model interfaces using object oriented modeling. Two partners involved in JEDI use or plan to use the Finite Volume Cubed-Sphere (FV3) dynamical core to produce weather forecasts; these are NASA's Global Modeling and Assimilation Office and NOAA's National Center for Environment Prediction. In this work we present an update on ongoing efforts to integrate the FV3 tangent linear and adjoint models into the prototype JEDI framework. We setup and run a simple cycled data assimilation experiment using 4DVAR on the cubed sphere grid and with the FV3 tangent linear and adjoint models. Development of the observation operators for JEDI is separately underway. Instead of using real observations a simplified set of simulated observations will be used. We discuss the steps required to bring the FV3 linearized model into the object oriented framework and consider what would be the computational requirements of running this configuration for an operational system. FV3 uses a small time-step to ensure that small scales are well resolved, however this presents design challenges when running 4DVAR with the adjoint. An approach to storing the FV3 model trajectory has been developed that maintains the flexibility of using automatic differentiation. We discuss how this approach is incorporated into the framework. Other important uses of adjoint models include computing observation impacts and singular vectors, we consider how these tools can be included in JEDI.

Holdaway, Daniel↗
NASA NTRS2018

Development of Gridded Innovations and Observations Supplement to MERRA-2

Atmospheric reanalysis have become an important source of data for weather and climate research, owing to the continuity of the data, but especially because of the multitude of observational data included (radiosondes, commercial aircraft, retrieved data products and radiances). However, the presence of assimilated observations can vary based on numerous factors, and so it is difficult or impossible for a researcher to say with any degree of certainty how many and what type of observations contributed to the reanalysis data they are using at any give point in time or space, or their contribution to the eventual analyzed fields. For example, quality control, transmission interruptions, and station outages can occasionally affect data availability. While orbital paths are known, drift in certain instruments, cloud clearing, and the large number of available instruments makes it challenging to know which satellite is observing any region at any point in the diurnal cycle. Furthermore, there is information from the statistics generated by the data assimilation that can help understand the model and the quality of the reanalysis. Typically, the assimilated observations and their innovations are saved in observation-space data formats and are not easily made available to reanalysis users.Here, we present an early version of a data set has been developed to make the MERRA-2 assimilated observations available for rapid and general use, by simplifying the data format. The observations are binned to a grid similar as MERRA-2 and saved as netCDF. This data collection includes the mean and number of observations in the bin as well as its variance. The data will also include the innovations from the data assimilation, the forecast departure and the analysis increment, as well as bias correction (for satellite radiances). In this paper, we present the data format (called MERRA-2 Gridded Innovations and Observations or GIO) and its strengths and limitations with some initial testing and validation of the methodology.

Bosilovich, Michael G.↗
NASA NTRS2019

The Ability of GeoCarb to Constrain the Interannual Variability of Carbon Gases over the Amazon

We perform a number of idealized assimilation experiments with the GEOS constituent data assimilation system to test the ability of GeoCarb retrievals of CO, CO2, and CH4 to constrain the interannual variability of these gases over the Amazon. Retrievals for instruments on other satellites which observe in similar channels (e.g. MOPITT, GOSAT, and OCO-2) are limited due to persistent cloud coverage. Given its ability to sample the same location multiple times in one day, the expectation is that GeoCarb retrievals will return more soundings than those from previous missions. The goal of the assimilation experiments is to understand which scanning strategies lead to the best sounding densities and thus have the best chance of constraining interannual variability in the carbon species. The experiments each begin by picking a given year at random from a nature run (i.e., a model simulation meant to represent the truth). The model fields are sampled according to a given strategy and then screened to account for cloud coverage. Next, we pick another year at random and assimilate the synthetic GeoCarb samples into the GEOS model for that year. The output of the assimilation, 6-hourly, 3D fields of each constituent, is then directly comparable to the nature run. This comparison allows us to evaluate the ability of GeoCarb measurements to constrain the interannual variability of each gas.

Weir, B.↗
NASA NTRS2019

Impact of Aquarius and SMAP Satellite Sea Surface Salinity Observations on Coupled El Niño/Southern Oscillation Forecasts

This study demonstrates the positive impact of including gridded Aquarius and SMAP sea surface salinity (SSS) into initialization of coupled forecasts for the tropical Indo-Pacific. An experiment that assimilates conventional ocean observations serves as the control. In a separate experiment, Aquarius and SMAP satellite SSS are additionally assimilated into the control initialization. Analysis of the initialization differences with the control indicates that SSS assimilation causes a freshening and shallowing of the mixed layer depth (MLD) near the equator and enhanced Kelvin wave amplitude. For each month from September 2011 to September 2018, 12 month coupled ENSO forecasts are initialized from both the control and satellite SSS assimilation experiments. The experiment assimilating Aquarius and SMAP SSS significantly outperforms the control relative to observed NINO3.4 sea surface temperature anomalies. This work highlights the importance of inclusion of satellite SSS for improving the initialization ENSO coupled forecasts.

Hackert, Eric C.↗
NASA NTRSDate not supplied

OCO-2 Satellite-imposed Constraints on Terrestrial Biospheric CO2 Flux over South Asia

The spatiotemporal variability of terrestrial biospheric carbon dioxide (CO2) flux over South Asia is poorly understood. The inverse model “top-down” CO2 flux estimates which rely on atmospheric CO2 observational data are impeded by sparse in situ measurements over this region. The Orbiting Carbon Observatory 2 (OCO-2) provides much needed global retrievals of column-average CO2 dry-air mole fraction (XCO2) with the finest spatial and temporal resolution and highest sensitivity to surface fluxes available to-date. This study conducted global inverse model simulations as part of the second OCO-2 Multi-model Intercomparison Project assimilating version 9 of the OCO-2 XCO2 retrievals in land nadir (LN) and land glint (LG) observing modes and global in situ (IS) measurements. The four-dimensional variational assimilation system with the GEOS-Chem global chemical transport model was used to estimate CO2 fluxes from 2015 to 2018. We then assessed the spatiotemporal variability of optimized CO2 Net ecosystem Exchange (NEE) fluxes aggregated over the South Asia region. The most robust result found by assimilating OCO-2 observations was the constraints imposed on the seasonal cycle of terrestrial biospheric fluxes over South Asia. The seasonality of South Asian NEE estimated by assimilating OCO-2 or global IS data showed a larger seasonal cycle compared to the current understanding of NEE in this region (represented by the prior NEE used in the model). The satellite-in situ joint inversion (IS + LN + LG) led to land flux seasonal amplitude (absolute magnitude of the difference between peak and trough of monthly mean values over a year) of 4.0 PgC yr-1, compared to the prior model (1.5 PgC yr-1). Moreover, OCO-2 data imposed a phase shift in the seasonal cycle, resulting in a large CO2 source in April and a large uptake in September. Most of the features in the seasonal cycle imposed by OCO-2 data are in agreement with previous “top-down” studies assimilating regional aircraft observations. Nevertheless, we recommend a closer examination of the novel features highlighted by the OCO-2 column data over South Asia in future studies by conducting regional inverse modelling with denser in situ and vertically-resolved regional observations along with satellite data.

OCO-2↗
NASA NTRS2021

Investigating the Utility of Hyperspectral Sounders in the 9.6 μm Band to Improve Ozone Analyses

Hyperspectral sounder brightness temperatures assimilated in the Global Earth Observing System Atmospheric Data Assimilation System (GEOS-ADAS) were previously limited to assimilating temperature and moisture. The ozone-sensitive 9.6 𝜇 m region is sensed by several hyperspectral sounders including AIRS (Atmospheric InfraRed Sounder), IASI (Infrared Atmospheric Sounding Interferometer), and CrIS (Cross-track Infrared Sounder). Direct assimilation of brightness temperatures in the 9.6 𝜇 m region have been operational at ECMWF for several years. With this study, similar improvements using the GEOS-ADAS are presented. Channels were selected from available operational subsets evaluating information content and minimizing inter-channel correlation. Additionally, information such as channel selections made by other studies, and vertical sensitivities of ozone and temperature were considered. The analyses produced show improvements verified against ozonesondes taken from SHADOZ (Southern Hemisphere Additional Ozonesondes) and WOUDC (World Ozone and Ultraviolet Data Center). This work will be added to the GEOS-ADAS and will provide an improved source of ozone data in NASA Global Modeling and Assimilation Office products.

Ozone↗
NASA NTRS2022

Using OSSEs to Evaluate the Impacts of Geostationary Infrared Sounders

An observing system simulation experiment (OSSE) was performed to assess the impact of assimilating hyperspectral infrared (IR) radiances from geostationary orbit on numerical weather prediction, with a focus on the proposed sounder on board the Geostationary Extended Observations (GeoXO) program’s central satellite. Infrared sounders on a geostationary platform would fill several gaps left by IR sounders on polar-orbiting satellites, and the increased temporal resolution would allow the observation of weather phenomena evolution. The framework for this OSSE was the Global Modeling and Assimilation Office (GMAO) OSSE system, which includes a full suite of meteorological observations. The experiment additionally assimilated four identical IR sounders from geostationary orbit to create a “ring” of vertical profiling observations. Based on the experimentation, assimilation of the IR sounders provided a beneficial impact on the analyzed mass and wind fields, particularly in the tropics, and produced an error reduction in the initial 24–48 h of the subsequent forecasts. Specific attention was paid to the impact of the GeoXO Sounder (GXS) over the contiguous United States (CONUS) as this is a region that is well-observed and as such difficult to improve. The forecast sensitivity to observation impact (FSOI) metric, computed across all four synoptic times over the CONUS, reveals that the GXS had the largest impact on the 24-h forecast error of the assimilated hyperspectral infrared satellite radiances as measured using a moist energy error norm. Based on this analysis, the proposed GXS has the potential to improve numerical weather prediction globally and over the CONUS.

OSSE↗
NASA NTRSDate not supplied

Enhancing CRTM All-Sky Simulations and Implementation of A New Active Sensor Module

Microwave observations play a very important role in improving the weather forecasts. Although these observations are routinely assimilated into NWP models in clear-sky conditions, assimilation of all-sky microwave observations is very limited. Two main factors contributing to this limitation are inaccuracy in the input cloud and hydrometeor profiles used as input to the radiative transfer model and also error in scattering calculations performed by the radiative transfer model itself. The Mie theory is used by many fast RT models to estimate the optical properties of single particles. The Mie theory assumes spherical shapes for ice or snow particles with mixture of air and ice. However, hydrometeors scattering radiation at microwave frequencies have different shapes, sizes, and orientations. Therefore, using Mie theory to determine their optical properties leads to large uncertainties in all-sky radiative transfer calculations. The discrete dipole approximation (DDA) which approximates the optical properties of large objects in terms of discrete dipoles has shown promise in calculating the scattering properties of particles with different shapes in the microwave frequencies. This presentation focuses on recent advancements in the CRTM scattering calculations for frozen hydrometeors in the microwave frequencies using the DDA technique. In addition to using stand-alone CRTM calculations using collocated ATMS and reanalysis profiles, the data assimilation experiments conducted using the NOAA FV3GFS forecast system are used to evaluate the scattering improvements. Additionally, the backscattering information from the DDA database was used to implement a radar simulator into CRTM. The radar operator takes advantage of CRTM different modules to calculate clouds absorption and scattering properties. In addition to the forward model both adjoint and tangent linear of the radar simulator are implemented and evaluated as well. The radar simulator is currently being tested within the JEDI/GEOS data assimilation framework to facilitate the assimilation of radar measurements such as CloudSat CPR and GPM DPR into the NASA GEOS model.

Isaac Moradi↗
NASA NTRSDate not supplied

CRTM Microwave Cloud Scattering Lookup Tables and Radar Simulator

Microwave observations play a very important role in improving the weather forecasts. Although these observations are routinely assimilated into NWP models in clear-sky conditions, assimilation of all-sky microwave observations is very limited. Two main factors contributing to this limitation are inaccuracy in the input cloud and hydrometeor profiles used as input to the radiative transfer model and also error in scattering calculations performed by the radiative transfer model itself. The Mie theory is used by many fast RT models to estimate the optical properties of single particles. The Mie theory assumes spherical shapes for ice or snow particles with mixture of air and ice. However, hydrometeors scattering radiation at microwave frequencies have different shapes, sizes, and orientations. Therefore, using Mie theory to determine their optical properties leads to large uncertainties in all-sky radiative transfer calculations. The discrete dipole approximation (DDA) which approximates the optical properties of large objects in terms of discrete dipoles has shown promise in calculating the scattering properties of particles with different shapes in the microwave frequencies. This presentation focuses on recent advancements in the CRTM scattering calculations for frozen hydrometeors in the microwave frequencies using the DDA technique. In addition to using stand-alone CRTM calculations using collocated ATMS and reanalysis profiles, the data assimilation experiments conducted using the NOAA FV3GFS forecast system are used to evaluate the scattering improvements. Additionally, the backscattering information from the DDA database was used to implement a radar simulator into CRTM. The radar operator takes advantage of CRTM different modules to calculate clouds absorption and scattering properties. In addition to the forward model both adjoint and tangent linear of the radar simulator are implemented and evaluated as well. The radar simulator is currently being tested within the JEDI/GEOS data assimilation framework to facilitate the assimilation of radar measurements such as CloudSat CPR and GPM DPR into the NASA GEOS model.

Isaac Moradi↗
NASA NTRSDate not supplied

Examination of L-Band Brightness Temperature Forecasts in the SMAP Level-4 Soil Moisture Analysis

The NASA Soil Moisture Active Passive (SMAP) mission [4] has been providing L-band (1.4 GHz) passive microwave brightness temperature (Tb) observations since April 2015. By assimilating the Tb observations into the NASA Catchment land surface model [5] using a spatially distributed ensemble Kalman filter, the NASA Global Modeling and Assimilation Office generates the SMAP Level-4 Soil Moisture (L4_SM) product, which provides global, 3-hourly, 9-km resolution estimates of surface (0-5 cm) and root-zone (0-100 cm) soil moisture with ~2.5-day latency for use in research and applications [6]. The L4_SM product also includes estimates of soil temperature, land surface fluxes, and assimilation diagnostics such as the model forecast and observed Tb values [7]. The output from the L4_SM system is routinely monitored by the L4_SM team. Such monitoring provides valuable information; instances of unusually large Tb observation-minus-forecast (O-F) residuals can indicate events during which soil moisture conditions are poorly described in the land modeling system [7,8]. For example, repeated occurrences of very large Tb O-F values in central Australia were traced back to deficiencies in the gauge-based precipitation product used in the land modeling system through L4_SM Version 5; this discovery prompted the use, beginning in Version 6, of satellite- and gauge-based precipitation observations outside of North America [9]. Recently, a statistical analysis of the Tb O-F residuals revealed systematic errors in the “tau-omega” L-band radiative transfer model that converts the land-model simulated soil moisture and temperature into the forecast Tb prior to the L4_SM analysis. The L4_SM analysis is built on the ensemble Kalman filter and assumes unbiased forecast errors. The presence of systematic Tb forecast errors could thus adversely impact the quality of the analyzed soil moisture. In this paper, we examine the Tb O-F residuals of the latest Version 7 L4_SM data (Science Version ID Vv7030 and Vv7032) [10,11,12]. In the L4_SM Version 7 algorithm, key parameters of the L-band radiative transfer model, including soil roughness, scattering albedo, and a (seasonally varying) climatology of vegetation opacity, are obtained from the SMAP Level-2 Radiometer retrieval product (Version 5) [1]. Additionally, we also use ground measurements of surface soil moisture and soil temperature from the SMAP core validation sites as in situ reference of soil conditions [2,3]. As part of the L4_SM system calibration, the seasonally-varying bias between the model forecast Tb and the observed values is removed prior to the assimilation of the SMAP Tb observations [7]. That is, the L4_SM system is designed to only correct errors in synoptic-scale and interannual variations from the long-term mean seasonal cycle while maintaining the model’s (potentially erroneous) climatology. The Tb O-F residuals examined here are thus computed using SMAP Tb observations after they are rescaled to the mean seasonal cycle of Tb from the modeling system without SMAP data assimilation. Consequently, the long-term average of the Tb O-F residuals has a global mean of only 0.13 K and locally small values, ranging from -1 to 3 K. Despite the small time-average values of the Tb O-F residuals, the model forecast Tb was nevertheless found to exhibit undesirable systematic errors. At some locations, the time-average Tb O-F values strongly depend on surface soil moisture (SM). At the Yanco SMAP core validation site, for example, the Tb O-F residuals typically range from 5 to 15 K under dry soil moisture conditions (SM < 0.15 m3 m-3) yet are predominantly negative under wet soil moisture conditions (SM > 0.25 m3 m-3), with values ranging from 0 to -40 K. This results in soil moisture analysis increments that persistently make the soil drier under dry SM conditions and persistently make the soil wetter under wet SM conditions, suggesting an error in the dynamic range of the simulated Tb, soil moisture or soil temperature. In this paper, we describe the higher-order systematic Tb forecast errors in more detail, examine their impact on the L4_SM product quality, and explore potential avenues to further improve the L4_SM algorithm.

Rolf Reichle↗
NASA NTRSDate not supplied

Machine-Learning-Based Adaptive Thinning of CrIS Radiances to Improve Global Tropical Cyclone Analysis and Forecasts

This work is focused on optimizing the assimilation of hyperspectral infrared (IR) radiances from the Cross-track Infrared Sounder (CrIS) with the goal of improving the representation of tropical cyclones (TCs) in global analyses and forecasts. Current operational assimilation systems rely on subsampling IR radiances on a regular thinning grid. A new and improved adaptive methodology based on machine learning (ML) recognizes TCs from geostationary satellite imagery and is implemented in the Goddard Earth Observing System (GEOS) model and data assimilation framework. The ML methodology is extensively trained on existing TC data sets and creates for each TC a dynamic mask, based on the evolving shape and life cycle of that specific event. Once a TC mask is created, a switch is then activated in the data assimilation system to alter the thinning, ingesting more CrIS radiances within the moving mask, thus increasing the TC sampling. After the TC dissipates, the assimilation of CrIS radiances reverts to normal data density. Results of TC segmentation provided by a state-of-the-art generative machine learning model known as the Denoising Diffusion Probabilistic Model (DDPM) are compared to the previously used U-Net model. The new approach surpasses the performance of the previously developed one. The methodology is applied to both clear-sky and cloud-cleared radiances. Benefits from the latter methodology, particularly in improving the structure of TCs and the intensity forecasts, are presented.

Oreste Reale↗
DOE OSTI2023

Explicitly determined sea ice emissivity and emission temperature over the Arctic for surface–sensitive microwave channels

Data assimilation of satellite microwave measurements is one of the important keys to improving weather forecasting over the Arctic region. However, the use of surface-sensitive microwave-sounding channel measurements for data assimilation or retrieval has been limited, especially during winter, due to the poorly constrained sea ice emissivity. In this study, aiming at more use of those channel measurements in the data assimilation, we propose an explicit method for specifying the surface radiative boundary conditions (namely emissivity and emitting layer temperature of snow and ice). These were explicitly determined with a radiative transfer model for snow and ice and with snow/ice physical parameters (i.e., snow/ice depths and vertical distributions of temperature, density, salinity, and grain size) simulated from the thermodynamically-driven snow/ice growth model. Here we conducted 1-D Var experiments in order to examine whether this approach can help to use the surface-sensitive microwave temperature channel measurements over the Arctic sea ice region for data assimilation. Results show that (1) the surface-sensitive microwave channels can be used in the 1-D Var retrieval, and (2) the specification of the radiative boundary condition at the surface using the snow/sea ice emission model can significantly improve the atmospheric temperature retrieval, especially in the lower troposphere (500 hPa to surface). The successful retrieval suggests that useful information can be extracted from surface-sensitive microwave-sounding channel radiances over sea ice surfaces through the explicit determination of snow/ice emissivity and emitting layer temperature.

1-D Var retrieval↗
DOE OSTI2025

Dayflow-PR: High-Resolution Streamflow Reanalysis for Puerto Rico, Version 1.0

This dataset presents a high-resolution historical streamflow reanalysis for NHDPlusV2 stream reaches across Puerto Rico (PR) spanning 1950 - 2019. The reanalysis is generated using the calibrated VIC-RAPID hydrologic modeling framework at the Hydrologic Unit Code Sub-basin (HUC08) scale, forced with sub-daily and daily meteorological forcings from Daymet. Runoff is simulated on 1- and 6-km grids, and the resulting total runoff is routed through the NHDPlusV2 river network using the RAPID routing model to produce Naturalized Streamflow Reanalysis. Where complete observational records are available over 1980 - 2019, streamflows are assimilated (substituted) and subsequently routed downstream through the river network to produce Assimilated Streamflow Reanalysis. The dataset includes streamflow outputs from eight distinct hydrologic modeling configurations along with key performanc evaluation metrics at daily and monthly scales, supporting a wide range of water resource applications. This dataset is derived to support the Non-Powered Dam Assessment, as well as 9505 Secure Water Assessment projects for the US Department of Energy (DOE) Water Power Technologies Office (WPTO). For further details, refer to Ghimire et al. (2023), Kao et al. (2024), and Ghimire et al. (2025).

13 HYDRO ENERGY↗
DOE OSTI2023

Data from: "Moisture rivals temperature in limiting photosynthesis by trees establishing beyond their cold-edge range limit under ambient and warmed conditions"

This archive contains data files that were used to draw conclusions in “Moisture rivals temperature in limiting photosynthesis by trees establishing beyond their cold-edge range limit under ambient and warmed conditions”, by Moyes et al., 2015. All field research was completed in common garden plots set up as part of the Alpine Treeline Warming Experiment (ATWE) on Niwot Ridge, Colorado, USA.There are two main data file formats in this archive: comma-separated values (.csv), and Microsoft Excel (.xls and .xlsx). .xlsx files can be read using Microsoft Excel and Google Sheets, and .csv files can be read using any simple text editor program, such as TextEdit (Mac) and Notepad (Windows). This .pdf data user’s guide can be read using Adobe Acrobat Reader, or any other compatible software. Seedling photographs and their corresponding leaf area-processed images are available in .jpg/.JPG image format, and can be opened using Preview (Mac) and Photos (Windows). To provide additional spatial context, two types of geospatial files are also published in this data package: ESRI shapefiles (.shp) and .kml files. Shapefiles are compatible with any GIS software able to read the file type (such as QGIS or ESRI’s ArcGIS suite), and .kml files can be opened with Google Earth or Google Maps. Figures 3 and 4 in the publication contain data from Moyes et al. 2013. This publication is cited in the References section in this archive, and data files can be accessed via the Alpine Treeline Warming Experiment project portal on ESS-DIVE. ---------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------- Climate change is altering plant species distributions globally, and warming is expected to promote uphill shifts in mountain trees. However, at many cold-edge range limits, such as alpine treelines in the western United States, tree establishment may be colimited by low temperature and low moisture, making recruitment patterns with warming difficult to predict.- We measured response functions linking carbon (C) assimilation and temperature- and moisture-related microclimatic factors for limber pine (Pinus flexilis) seedlings growing in a heating × watering experiment within and above the alpine treeline. We then extrapolated these response functions using observed microclimate conditions to estimate the net effects of warming and associated soil drying on C assimilation across an entire growing season.- Moisture and temperature limitations were each estimated to reduce potential growing season C gain from a theoretical upper limit by 15–30% (c. 50% combined). Warming above current treeline conditions provided relatively little benefit to modeled net assimilation, whereas assimilation was sensitive to either wetter or drier conditions.- Summer precipitation may be at least as important as temperature in constraining C gain by establishing subalpine trees at and above current alpine treelines as seasonally dry subalpine and alpine ecosystems continue to warm.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

An entropy-based debiasing approach to quantifying experimental coverage for novel applications of interest in the nuclear community

This manuscript proposes a novel information-theoretic approach to the quantification of experimental relevance, i.e., coverage, to achieve optimal data assimilation results for nuclear engineering applications. Specifically, this work posits the need for a new metric, called coverage (q C ) of an application’s quantity of interest, i.e., eigenvalue or power peaking for an advanced reactor concept, defined herein as the theoretically maximum achievable reduction in the quantity’s uncertainty given measurements from a pool of experiments in a manner that is independent of the data assimilation procedure employed. Currently, reduction in a quantity’s uncertainty is strongly biased by the underlying assumptions of the assimilation procedure to account for the under-determined nature of such problems and the similarity criterion employed to identify relevant experiments. To address this challenge, this work has developed a coverage metric, q C , based on mutual information, which establishes a new conceptual framework for assessing coverage, one that is independent of the model parameters and responses degree of variations in both the experimental and application domains, i.e., linear vs non-linear, and their prior uncertainty distributions, i.e., Gaussian vs. non-Gaussian. The q C is an entropic measure capable of addressing coverage for general nonlinear problems with non-Gaussian uncertainties and inclusive of the measurement uncertainties from multiple experiments. Numerical experiments from manufactured analytical problems as well as a set of benchmarks from the ICSBEP handbook are employed to demonstrate its theoretical and practical performance as compared to the c k -based experiment selection methodology, commonly employed in the neutronic community. The manuscript then employs other well-known adaptations to existing data assimilation methodologies for nonlinear and non-Gaussian problems capable of achieving the coverage posited by q C .

Bayesian data assimilation↗
NASA NTRS2013

Simultaneous Estimation of Model State Variables and Observation and Forecast Biases Using a Two-Stage Hybrid Kalman Filter

In this paper, we present a two-stage hybrid Kalman filter to estimate both observation and forecast bias in hydrologic models, in addition to state variables. The biases are estimated using the discrete Kalman filter, and the state variables using the ensemble Kalman filter. A key issue in this multi-component assimilation scheme is the exact partitioning of the difference between observation and forecasts into state, forecast bias and observation bias updates. Here, the error covariances of the forecast bias and the unbiased states are calculated as constant fractions of the biased state error covariance, and the observation bias error covariance is a function of the observation prediction error covariance. In a series of synthetic experiments, focusing on the assimilation of discharge into a rainfall-runoff model, it is shown that both static and dynamic observation and forecast biases can be successfully estimated. The results indicate a strong improvement in the estimation of the state variables and resulting discharge as opposed to the use of a bias-unaware ensemble Kalman filter. Furthermore, minimal code modification in existing data assimilation software is needed to implement the method. The results suggest that a better performance of data assimilation methods should be possible if both forecast and observation biases are taken into account.

Data Assimilation↗