Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “observations model”

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 19 records
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Giant Cloud Condensation Nuclei Facilitate Drizzle Formation in Stratocumulus—Insights From a Combined Observation‐Modeling Framework

The mechanism for initiating drizzle drop remains a gap in the current understanding of warm rain formation. One prevalent hypothesis suggests that the presence of Giant Cloud Condensation Nuclei (GCCN) generates drizzle‐sized drops necessary to trigger the Collision‐Coalescence (C‐C) process. Here, in this study, this hypothesis is investigated using a novel framework that integrates in situ observations, remote sensing measurements, and idealized models. Results show that GCCN can efficiently generate drizzle drops through condensation, producing a broad Droplet Size Distribution (DSD) comparable to in situ observations. The large drizzle drop and broad DSD strongly facilitate C‐C, further accelerating drizzle initiation. To compare with observation, the model‐generated DSDs are used to generate radar Doppler spectra where radar reflectivity and Doppler skewness is estimated. The simulated radar quantities correspond well with radar observations, providing critical evidence for the GCCN‐induced drizzle initiation mechanism.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Integrating New Arctic Plant Functional Types in a Land Surface Model Using Above- and Belowground Field Observations: Modeling Archive

This Modeling Archive is in support of an NGEE Arctic publication "Integrating New Arctic Plant Functional Types in a Land Surface Model Using Above- and Belowground Field Observations". We used field observations of biomass and leaf traits (Salmon et al., 2019b) across a gradient of plant communities on the Seward Peninsula in northwest Alaska (Breen et al., 2020; Iversen et al., 2017) to replace the original ELM configuration with nine, arctic-specific PFTs. Original arctic PFTs in the model were (1) broadleaf deciduous boreal shrub and (2) C3 arctic grass. New PFTs were (1) lichen, (2) bryophyte, (3) evergreen dwarf shrub, (4) deciduous dwarf shrub, (5) deciduous low shrub, (6) deciduous low to tall shrub, (7) alder shrub, (8) forb, and (9) graminoid. New PFTs were parameterized and evaluated using site-level trait and biomass measurements (Salmon et al., 2019b), root C:N ratios and rooting depth patterns (Salmon et al., 2019c), and root biomass (Salmon et al., 2020). This archive contains forcing data and model output from ELM simulations conducted at the Kougarok Hillslope site (Kougarok Road Mile Marker 64). Simulations were conducted for the Kougarok Hillslope site using meteorological driving data from the Scenarios Network for Alaska and Arctic Planning (SNAP) downscaled climate projection dataset using NCAR-CCSM forcing (Bieniek et al., 2020 dataset; Walsh et al., 2018). The archive includes three versions of model simulations representing a progression from the original model configuration to an Arctic-specific configuration (see Output section below). Simulation 1 used grid cell data from a global E3SM configuration, including the fractional area of E3SM PFTs assigned to the grid cell containing the Kougarok Hillslope site in global simulations. This simulation did not distinguish between different plant communities on the Kougarok Hillslope, but instead used a single point simulation to represent the entire area. Simulation 2 used default E3SM PFT definitions combined with adjusted depth to bedrock for plant communities with shallow rocky layers to represent the role of abiotic soil factors in driving site differences (based on depth measurements of Iversen et al., 2019a), and adjusted the relative areas of the broadleaf deciduous boreal shrub and C3 arctic grass PFTs to reflect the observed spatial coverage of shrub and non-shrub PFTs across Kougarok Hillslope plant communities. Simulation 3 used the new arctic PFT definitions and parameterizations based on the vegetation types present at the site, along with the adjusted soil depths used in Simulation 2. Included are *.pdf, *.nc (NetCDF), *tar.gz, *.py, and *.bash files. We would like to thank Mary's Igloo Native Corporation for allowing us to perform this research on their land and giving us the opportunity to share this research with their community. Thanks to Shawn Serbin and Alistair Rogers for helpful discussions. Holly Vander Stel, Joanne Childs, and Stan Wullschleger helped with collection of field data. The Next-Generation Ecosystem Experiments: Arctic (NGEE Arctic), was a research effort to reduce uncertainty in Earth System Models by developing a predictive understanding of carbon-rich Arctic ecosystems and feedbacks to climate. NGEE Arctic was supported by the Department of Energy's Office of Biological and Environmental Research. The NGEE Arctic project had two field research sites: 1) located within the Arctic polygonal tundra coastal region on the Barrow Environmental Observatory (BEO) and the North Slope near Utqiagvik (Barrow), Alaska and 2) multiple areas on the discontinuous permafrost region of the Seward Peninsula north of Nome, Alaska. Through observations, experiments, and synthesis with existing datasets, NGEE Arctic provided an enhanced knowledge base for multi-scale modeling and contributed to improved process representation at global pan-Arctic scales within the Department of Energy's Earth system Model (the Energy Exascale Earth System Model, or E3SM), and specifically within the E3SM Land Model component (ELM).

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Learning to Count Grave Sites for Cemetery Observation Models With Satellite Imagery

Understanding how people occupy open spaces is important for research in support of population modeling, policy, national security, emergency response, and sustainability. For the past decade, there has been an increase in research toward capturing and reporting population dynamics and patterns of life at the building level and in some open public spaces such as cemeteries and parks. This is done through observation models developed from local sociocultural information acquired at various spatiotemporal scales to inform night, day, and episodic population occupancy estimates (people/1000 sq ft). Sociocultural information for cemeteries and parks is scarcely available and often collected manually. The process is not only marred by inconsistencies but is laborious and time consuming. In this study, we leverage convolutional neural networks (CNNs) and satellite imagery to derive grave site counts as proxy variables to support scalable and accurate sociocultural data required in a population observation model. Through a hybrid workflow (weak localization plus regression model), we characterize a large scale automation process to counting of grave sites. We evaluate and demonstrate the efficacy of proposed workflow using out-of-data set large satellite imagery and establish its broader impact on cemetery observation models.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗
DOE OSTI2025

Hunga Tonga–Hunga Ha′apai Volcano Impact Model Observation Comparison (HTHH-MOC) project: experiment protocol and model descriptions

The 2022 Hunga volcanic eruption injected a significant amount of water vapor and a moderate amount of sulfur dioxide into the stratosphere, causing observable responses in the climate system. We have developed a model–observation comparison project to investigate the evolution of volcanic water and aerosols and their impacts on atmospheric dynamics, chemistry, and climate, using several state-of-the-art chemistry climate models. The project goals are (1) to evaluate the current chemistry–climate models to quantify their performance in comparison to observations and (2) to understand atmospheric responses in the Earth system after this exceptional event and investigate the potential impacts in the projected future. To achieve these goals, we designed specific experiments for direct comparisons to observations, for example from balloons and the Microwave Limb Sounder satellite instrument. Experiment 1 consists of two sets of free-running ensemble experiments from 2022 to 2031: one with fixed sea-surface temperatures and sea ice and one with coupled ocean. These experiments will help to understand the long-term evolution of water vapor and aerosols; quantify HTHH effects on stratospheric and mesospheric temperatures, dynamics, and transport; understand the impact of dynamic changes on ozone chemistry; quantify the net radiative forcings; and evaluate any surface climate impact. Experiment 2 is a nudged-run experiment from 2022 to 2023 using observed meteorology. To allow participation of more climate models with varying complexities of aerosol simulation, we include two sets of simulations in Experiment 2: Experiment 2a is designed for models with internally generated aerosol, while Experiment 2b is designed for models using prescribed aerosol surface area density. This experiment will help to analyze H 2 O and aerosol evolution, quantify the net radiative forcings, understand the impacts on mid-latitude and polar O 3 chemistry, and allow close comparisons with observations.

Zhu, Yunqian [Univ. of Colorado, Boulder, CO (Unit↗
DOE OSTI2025

Application of PRIM for understanding patterns in carbon dioxide model-observation differences

Reducing uncertainties in regional carbon balances requires a better understanding of CO 2 transport in synoptic weather systems. Here, we apply the Patient Rule Induction Method (PRIM), a data-mining method to identify high-density regions for a target-class within an input parameter space, to airborne observations of potential temperature, wind speed, water vapor mixing ratio, and CO 2 dry mol fraction gathered during the Atmospheric Carbon and Transport (ACT)-America Summer 2016 and Winter 2017 campaigns. ACT observations were targeted at expert-designated cases of fair weather and near-frontal warm and cold sector air at atmospheric boundary-layer, lower-, and higher free tropospheric levels (ABL, LFT, and HFT, respectively). We investigate atmospheric characteristics of these pre-defined cases and associated CO 2 model-observation-differences in the mesoscale WRF-Chem model. PRIM results separate winter- and summertime observations as well as observations from ABL, LFT, and HFT with enrichment factors of 4.0–20.5 inside the PRIM box compared to the entire dataset but cannot distinguish between near-frontal warm and cold sector observations in the higher free troposphere. Analyzing of the parameter space constrained by PRIM, we find that large magnitude model observation differences preferentially associated with times when atmospheric conditions are less typical. This association suggests that PRIM could provide a useful tool for isolating atmospheric conditions with large-magnitude and non-Gaussian CO 2 -residuals for targeted transport model evaluation and to potentially improve inversion results during synoptically active periods.

Gerken, Tobias [James Madison Univ., Harrisonburg,↗
DOE OSTI2021

Understanding water and energy fluxes in the Amazonia: Lessons from an observation‐model intercomparison

Abstract Tropical forests are an important part of global water and energy cycles, but the mechanisms that drive seasonality of their land‐atmosphere exchanges have proven challenging to capture in models. Here, we (1) report the seasonality of fluxes of latent heat (LE), sensible heat ( H ), and outgoing short and longwave radiation at four diverse tropical forest sites across Amazonia—along the equator from the Caxiuanã and Tapajós National Forests in the eastern Amazon to a forest near Manaus, and from the equatorial zone to the southern forest in Reserva Jaru; (2) investigate how vegetation and climate influence these fluxes; and (3) evaluate land surface model performance by comparing simulations to observations. We found that previously identified failure of models to capture observed dry‐season increases in evapotranspiration (ET) was associated with model overestimations of (1) magnitude and seasonality of Bowen ratios (relative to aseasonal observations in which sensible was only 20%–30% of the latent heat flux) indicating model exaggerated water limitation, (2) canopy emissivity and reflectance (albedo was only 10%–15% of incoming solar radiation, compared to 0.15%–0.22% simulated), and (3) vegetation temperatures (due to underestimation of dry‐season ET and associated cooling). These partially compensating model‐observation discrepancies (e.g., higher temperatures expected from excess Bowen ratios were partially ameliorated by brighter leaves and more interception/evaporation) significantly biased seasonal model estimates of net radiation ( R n ), the key driver of water and energy fluxes (LE ~ 0.6 R n and H ~ 0.15 R n ), though these biases varied among sites and models. A better representation of energy‐related parameters associated with dynamic phenology (e.g., leaf optical properties, canopy interception, and skin temperature) could improve simulations and benchmarking of current vegetation–atmosphere exchange and reduce uncertainty of regional and global biogeochemical models.

Restrepo‐Coupe, Natalia↗
DOE OSTI2022

Space weather observations, modeling, and alerts in support of human exploration of Mars

Space weather observations and modeling at Mars have begun but they must be significantly increased to support the future of Human Exploration on the Red Planet. A comprehensive space weather understanding of a planet without a global magnetosphere and a thin atmosphere is very different from our situation at Earth so there is substantial fundamental research remaining. It is expected that the development of suitable models will lead to a comprehensive operational Mars space weather alert (MSWA) system that would provide rapid dissemination of information to Earth controllers, astronauts in transit, and those in the exploration zone (EZ) on the surface by producing alerts that are delivered rapidly and are actionable. To illustrate the importance of such a system, we use a magnetohydrodynamic code to model an extreme Carrington-type coronal mass ejection (CME) event at Mars. The results show a significant induced surface field of nearly 3,000 nT on the dayside that could radically affect unprotected electrical systems that would dramatically impact human survival on Mars. Other associated problems include coronal mass ejection (CME) shock-driven acceleration of solar energetic particles producing large doses of ionizing radiation at the Martian surface. In summary, along with working more closely with international partners, the next Heliophysics Decadal Survey must include a new initiative to meet expected demands for space weather forecasting in support of humans living and working on the surface of Mars. It will require significant effort to coordinate NASA and the international community contributions.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2020

Fingerprints of external forcings on Sahel rainfall: aerosols, greenhouse gases, and model-observation discrepancies

Over the 20th and 21st centuries, both anthropogenic greenhouse gas increases and changes in anthropogenic aerosols have affected rainfall in the Sahel. Using multiple characteristics of Sahel precipitation, we construct a multivariate fingerprint that allows us to distinguish between the model-predicted responses to greenhouse gases and anthropogenic aerosols. Models project the emergence of a detectable signal of aerosol forcing in the middle of the 20th century and a detectable signal of greenhouse gas forcing at the beginning of the 21st. However, the signals of both aerosol and greenhouse gas forcing in observations emerge earlier and are stronger than in the models, far stronger in the case of aerosols. The similarity between the response to aerosol forcing and the leading mode of internal variability makes it difficult to attribute this model-observation discrepancy to errors in the forcing, errors in the forced response, model inability to capture the amplitude of internal variability, or some combination of these. For greenhouse gases, however, the forced response is distinct from internal variability as estimated by models, and the observations are largely commensurate with the model projections.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Machine Learning for a-posteriori model-observed data fusion to enhance predictive value of ESM output

The proposed research is consistent with focal areas 2 and 3: when successful, this effort will enhance the predictive value of Earth System Model (ESM) output by identifying and quantifying its systematic deviations from available observations and by correcting and recalibrating to existing data. The product is a hybrid process & data-driven modeling tool whose estimation can also bring insights through pattern discovery recognition on model deficiencies on the one hand and data gaps on the other.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

First Direct Measurement of Mg (α,p) 25 Al and Implications for X-Ray Burst Model-Observation Comparisons

Type-I x-ray bursts can reveal the properties of an accreting neutron star system when compared with astrophysics model calculations. However, model results are sensitive to a handful of uncertain nuclear reaction rates, such as 22 Mg(α , p). We report the first direct measurement of 22 Mg(α , p), performed with the Active Target Time Projection Chamber. The corresponding astrophysical reaction rate is orders of magnitude larger than determined from a previous indirect measurement in a broad temperature range. Overall, our new measurement suggests a less-compact neutron star in the source GS1826-24.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2026

Model-observation discrepancies in Arctic moisture intrusions: causes and pathways for improved simulation

Arctic moisture intrusions (MIs), narrow filaments of strong moisture transport, are key drivers of poleward moisture flux and Arctic weather extremes, yet their representation in climate models is poorly understood. Using a new Arctic MI detection algorithm, we document persistent biases across three CMIP generations (CMIP3–CMIP6): models overestimate MI occurrence over the Pacific sector and underestimate it over the Atlantic sector. These errors stem from misrepresented midlatitude westerly jets, with an equatorward North Atlantic jet associated with too few Atlantic MIs, and a poleward, weakened North Pacific jet linked to too many Pacific MIs. Experiments that correct sea surface temperature and sea ice concentration biases and increase atmospheric resolution improve jet structure and MI statistics, while a cloud-locking simulation indicates that better high-frequency cloud–radiation–circulation interactions can yield further gains. Our results clarify pathways to reducing long-standing MI and jet biases, providing guidance for improving simulations of Arctic and midlatitude climate.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

PYOED: AN ETENSIBLE SUITE FOR DATA ASSIMILATION AND MODEL-CONSTRAINED OPTIMAL DESIGN OF EXPERIMENTS

SF-23-005 PyOED is a highly extensible scientific package that enables developing and testing model-constrained optimal experimental design (OED) for inverse problems. Specifically, PyOED aims to be a comprehensive Python toolkit for model-constrained OED. The package targets scientists and researchers interested in understanding the details of OED formulations and approaches. It is also meant to enable researchers to experiment with standard and innovative OED technologies with a wide range of test problems (e.g., simulation models). OED, inverse problems (e.g., Bayesian inversion), and data assimilation (DA) are closely related research fields, and their formulations overlap significantly. Thus, PyOED is continuously being expanded with a plethora of Bayesian inversion, DA, and OED methods as well as new scientific simulation models, observation error models, and observation operators. These pieces are added such that they can be permuted to enable testing OED methods in various settings of varying complexities. The PyOED core is completely written in Python and utilizes the inherent object-oriented capabilities; however, PyOED is meant to be extensible rather than scalable. Specifically, PyOED is developed to ``enable rapid development and benchmarking of OED methods with minimal coding effort and to maximize code reutilization.'' PyOED will be continuously expanded with a plethora of Bayesian inversion, DA, and OED methods as well as new scientific simulation models, observation error models, and observation operators.

ATTIA, AHMEDMOHAMED↗
DOE OSTI2022

Interaction between cloud–radiation, atmospheric dynamics and thermodynamics based on observational data from GoAmazon 2014/15 and a cloud-resolving model

Observational meteorological data from the field experiment GoAmazon 2014/15 and data from numerical simulations with the cloud-resolving model (CRM) called the System for Atmospheric Modeling (SAM) are used to study the interaction between the cloudiness–radiation as well as the atmospheric dynamics and thermodynamics variables for a site located in the central Amazon region (–3.2°S, –60.6°W) during the wet and dry periods. The main aims are to (a) analyze the temporal series of the integrated cloud fraction, precipitation rate and downward shortwave flux as well as (b) to determine the relationship between the integrated cloud fraction, radiative fluxes and large-scale variable anomalies as a function of the previous day's average. The temporal series of the integrated cloud fraction, precipitation rate and downward shortwave flux from SAM simulations showed physical consistency with the observations from GoAmazon 2014/15. Shallow and deep convection clouds show to have a meaningful impact on radiation fluxes in the Amazon region during wet and dry periods. Anomalies of large-scale variables (relative to the previous day's average) are physically associated with cloud formation, evolution and dissipation. SAM consistently simulated these results, where the cloud fraction vertical profile shows a pattern very close to the observed data (cloud type). Additionally, the integrated cloud fraction and large-scale variable anomalies, as a function of the previous day's average, have a good correlation. These results suggest that the memory of the large-scale dynamics from the previous day can be used to estimate the cloud fraction as well as the water content, which is a variable of the cloud itself. In general, the SAM satisfactorily simulated the interaction between cloud–radiation as well as dynamic and thermodynamic variables of the atmosphere during the periods of this study, being able to obtain atmospheric variables that are impossible to obtain in an observational way.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Revisiting AMOC Transport Estimates From Observations and Models

Abstract Reference level assumptions used to calculate the Atlantic meridional overturning circulation transports at the RAPID and Meridional Overturning Variability Experiment (MOVE) observing arrays are revisited in an eddying ocean model. Observational transport calculation methods are complemented by several alternative approaches. At RAPID, the model transports from the observational method and the model truth (based on the actual model velocities) agree well in their mean and variability. There are substantial differences among the transport estimates obtained with various methods at the MOVE site. These differences result from relatively large and time‐varying reference velocities at depth in the model, not supporting a level‐of‐no‐motion. The methods that account for these reference velocities properly at MOVE produce transports that are in good agreement with the model truth. In contrast with the observational estimates, the model transport trends at MOVE and RAPID largely agree with each other on pentadal to multi‐decadal time scales.

Danabasoglu, Gokhan↗
DOE OSTI2022

Cross-correlation of Dark Energy Survey Year 3 lensing data with ACT and P l a n c k thermal Sunyaev-Zel’dovich effect observations. II. Modeling and constraints on halo pressure profiles

Hot, ionized gas leaves an imprint on the cosmic microwave background via the thermal Sunyaev-Zel’dovich (tSZ) effect. The cross-correlation of gravitational lensing (which traces the projected mass) with the tSZ effect (which traces the projected gas pressure) is a powerful probe of the thermal state of ionized baryons throughout the Universe and is sensitive to effects such as baryonic feedback. In a companion paper (Gatti et al. Phys. Rev. D 105, 123525 (2022)), we present tomographic measurements and validation tests of the cross-correlation between Galaxy shear measurements from the first three years of observations of the Dark Energy Survey and tSZ measurements from a combination of Atacama Cosmology Telescope and Planck observations. In this work, we use the same measurements to constrain models for the pressure profiles of halos across a wide range of halo mass and redshift. We find evidence for reduced pressure in low-mass halos, consistent with predictions for the effects of feedback from active Galactic nuclei. We infer the hydrostatic mass bias (B ≡ M 500c /M SZ from our measurements, finding B = 1.8 ± 0.1 when adopting the Planck-preferred cosmological parameters. We additionally find that our measurements are consistent with a nonzero redshift evolution of B, with the correct sign and sufficient magnitude to explain the mass bias necessary to reconcile cluster count measurements with the Planck-preferred cosmology. Our analysis introduces a model for the impact of intrinsic alignments (IAs) of galaxy shapes on the shear-tSZ correlation. We show that IA can have a significant impact on these correlations at current noise levels.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2020

Observations for Model Intercomparison Project (Obs4MIPs): status for CMIP6

The Observations for Model Intercomparison Project (Obs4MIPs) was initiated in 2010 to facilitate the use of observations in climate model evaluation and research, with a particular target being the Coupled Model Intercomparison Project (CMIP), a major initiative of the World Climate Research Programme (WCRP). To this end, Obs4MIPs (1) targets observed variables that can be compared to CMIP model variables; (2) utilizes dataset formatting specifications and metadata requirements closely aligned with CMIP model output; (3) provides brief technical documentation for each dataset, designed for nonexperts and tailored towards relevance for model evaluation, including information on uncertainty, dataset merits, and limitations; and (4) disseminates the data through the Earth System Grid Federation (ESGF) platforms, making the observations searchable and accessible via the same portals as the model output. Taken together, these characteristics of the organization and structure of obs4MIPs should entice a more diverse community of researchers to engage in the comparison of model output with observations and to contribute to a more comprehensive evaluation of the climate models. At present, the number of obs4MIPs datasets has grown to about 80; many are undergoing updates, with another 20 or so in preparation, and more than 100 are proposed and under consideration. A partial list of current global satellite-based datasets includes humidity and temperature profiles; a wide range of cloud and aerosol observations; ocean surface wind, temperature, height, and sea ice fraction; surface and top-of-atmosphere longwave and shortwave radiation; and ozone (O 3 ), methane (CH 4 ), and carbon dioxide (CO 2 ) products. A partial list of proposed products expected to be useful in analyzing CMIP6 results includes the following: alternative products for the above quantities, additional products for ocean surface flux and chlorophyll products, a number of vegetation products (e.g., FAPAR, LAI, burned area fraction), ice sheet mass and height, carbon monoxide (CO), and nitrogen dioxide (NO 2 ). While most existing obs4MIPs datasets consist of monthly-mean gridded data over the global domain, products with higher time resolution (e.g., daily) and/or regional products are now receiving more attention. Along with an increasing number of datasets, obs4MIPs has implemented a number of capability upgrades including (1) an updated obs4MIPs data specifications document that provides additional search facets and generally improves congruence with CMIP6 specifications for model datasets, (2) a set of six easily understood indicators that help guide users as to a dataset's maturity and suitability for application, and (3) an option to supply supplemental information about a dataset beyond what can be found in the standard metadata. With the maturation of the obs4MIPs framework, the dataset inclusion process, and the dataset formatting guidelines and resources, the scope of the observations being considered is expected to grow to include gridded in situ datasets as well as datasets with a regional focus, and the ultimate intent is to judiciously expand this scope to any observation dataset that has applicability for evaluation of the types of Earth system models used in CMIP.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

PyOED: An Extensible Suite for Data Assimilation and Model-Constrained Optimal Design of Experiments

This article describes PyOED, a highly extensible scientific package that enables developing and testing model-constrained optimal experimental design (OED) for inverse problems. Specifically, PyOED aims to be a comprehensive Python toolkit for model-constrained OED. The package targets scientists and researchers interested in understanding the details of OED formulations and approaches. It is also meant to enable researchers to experiment with standard and innovative OED technologies with a wide range of test problems (e.g., simulation models). OED, inverse problems (e.g., Bayesian inversion), and data assimilation (DA) are closely related research fields, and their formulations overlap significantly. Thus, PyOED is continuously being expanded with a plethora of Bayesian inversion, DA, and OED methods as well as new scientific simulation models, observation error models, and observation operators. These pieces are added such that they can be permuted to enable testing OED methods in various settings of varying complexities. The PyOED core is completely written in Python and utilizes the inherent object-oriented capabilities; however, the current version of PyOED is meant to be extensible rather than scalable. Specifically, PyOED is developed to “enable rapid development and benchmarking of OED methods with minimal coding effort and to maximize code reutilization.” This article provides a brief description of the PyOED layout and philosophy and provides a set of exemplary test cases and tutorials to demonstrate the potential of the package.

97 MATHEMATICS AND COMPUTING↗