Search NASA⌕ Search

SEARCH · Search NASA

Results for “Climate Model Selection”

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 55 records · Page 3

A Systematic Framework for Investigating Climate Impact-Drivers in NASA GISS ModelE Configurations: Towards Improving ESM Development for Impact Applications

Earth System Model (ESM) development benefits from attention to the downstream applications for nature and society. The evaluation of ESMs often focuses on a limited set of climatic impact-drivers (CIDs)1 and numerous large-scale parameters that are not directly relevant for impact analysis. To address this limitation, we propose a comprehensive framework for impact-relevant model version comparison. This framework explores both the intended and unintended consequences of shifts in model components, parameters, and resolutions, emphasizing the perspectives of the climate impacts community. We have selected a core set of CIDs as a basis for comparing model versions.

earth system models↗

Icehouse Effect: A Selective Arctic Cooling Trend Current Models are Missing

The icehouse effect is a hypothesized climate feedback mechanism which could result in human-caused surface cooling trends in polar regions. Once understood in detail, it becomes apparent that these trends, which are discernable in the literature, but have been largely dismissed, do not conflict with the consensus assessment of the evidence, which infers century-scale Arctic warming. In fact, confirmation of the hypothesis would substantially strengthen the argument that there is a detectable human influence on today's climate. This apparent enigma is resolved only through careful attention to the detail of the hypothesis and the data supporting it. The posited surface cooling is entirely dependent on the existence of climate warming in layers capping the stable boundary layer. Also, the cooling is not pandemic, but is selective. It is readily revealed in properly sorted data by making use of the principles of micrometeorological similarity. Specifically, the cooling is manifest under a range of favorable turbulence conditions which can develop and disappear locally on time scales of minutes to hours because of the intrinsically intermittent nature of stable boundary layer turbulence. Because of the fine-scale nature of the processes which produce the cooling, modeling it is a difficult proposition. Vertical resolution on the order of 1 meter is required. Adequate models of intermittent surface fluxes coupled with radiation exchange do not currently exist, not as parameterizations for aggregated systems, nor in large eddy simulation (LES) models. This presentation will introduce the theory. An important testable null hypothesis emerges: the icehouse effect produces a unique signature or "fingerprint" which could not be produced by any other known process. The presence of this signature will be demonstrated using nearly all available Arctic temperature observations. Its aggregate effect is clearly found in Arctic monthly surface temperature trends when sorted by climatological stability. Using all available Arctic rawinsonde ascents - about 1.1 million profiles, "frozen moments" of the icehouse processes are captured in various states. Because turbulent time scales are so short in the stable boundary layer. each of these snapshots can be treated as independent -- their chronology is irrelevant. Micrometeorological similarity is invoked to reassemble the soundings into bins of similar stability and it is in a wide, coherent range of these stability bins where the cooling effect is revealed.

Wetzel, Peter J.↗

Unprecedented Shifts in Hydrology Are Emerging Across California's Critical Basins: An Evaluation From 0.5 to 3.5°C

With advances in climate models and downscaling techniques, stakeholders anticipate high-resolution analysis to inform regional to local changes in water management. Here, we produce hydrologic projections from an ensemble of Earth System Models (ESMs) that were selected and downscaled to support California's 5th Climate Assessment. An ensemble of 19 ESMs was downscaled to a 3-km resolution across California using a statistical-dynamical downscaling approach and subsequently run through two calibrated hydrology models. Although California has been extensively studied in the context of climate change, we provide the first evaluation of the warming thresholds at which hydroclimate metrics demonstrate statistically significant shifts. We show that present-day to near-term warming levels in Klamath and Northern Sierra Nevada basins, which serve as a critical source of water for California, show statistically significant decreases in snowfall and peak snowpack and associated decreases in summer snowmelt and runoff. More generally, shifts in these hydroclimate metrics occur for intermediate elevation basins ranging from 1,315 to 1,455 m (4,314.3–4,773.6 ft), while the warming level of emergence is delayed for lower and higher elevation basins. We also find that several basins already demonstrate significant increases in 5- to 100-year runoff intensities, primarily due to the increasing influence of precipitation falling as rain. Hydroclimate metrics with trends that demonstrate near-term warming levels of emergence are reflected in reanalysis-based observations, suggesting California is entering a fundamentally different hydroclimate regime. While this will likely stress California's water management, the research provided can support when to implement adaptation efforts.

Bass, B. [Univ. of California, Los Angeles, CA (Un↗

A study of model parameters associated with the urban climate using HCMM data

Progress in the study of the intensity of the urban heat island is reported. The intensity of the heat island is commonly defined as the temperature difference between the center of the city and the surrounding suburban and rural regions. The intensity is considered as a function of changes in the season and changes in meteorological conditions in order to derive various parameters which may be used in numerical models for urban climate. Twelve case studies were selected and CCT's were ordered. In situ data was obtained from sixteen stations scattered about the city of St. Louis. Upper-air meteorological data were obtained and the water vapor and the temperature data were processed. Atmospheric transmissivities were computed for each of the case studies.

Source record↗

Characteristics and Global Impact of Aerosols from Southern Africa and Eastern Asia

Supported mainly by the NASA GACP and ACMAP, we have made significant progress in the global modeling of tropospheric aerosols and their precursors in the past few years, especially in the development of the GOCART model, simulation of anthropogenic and natural aerosols, data analysis of field observations and satellite retrievals, assessment of global and regional budgets, estimate of aerosol direct radiative forcing, and aerosol forecasting and data analysis for the ACE-Asia field experiment. Our results and findings are summarized in Chin et al. Model calculated multiple-year optical thickness for individual and total aerosols are at internet. These results have been frequently used by other groups, for example, to impose initial conditions for regional models, provide dust source functions for other global models, supply aerosol fields for chemistry and climate models, help data group interpret their measurements, select monitoring sites for ground observation network, and assist satellite retrievals.

Chin, Mian↗

Habitable Moist Atmospheres on Terrestrial Planets near the Inner Edge of the Habitable Zone Around M Dwarfs

Terrestrial planets in the habitable zones (HZs) of low-mass stars and cool dwarfs have received significant scrutiny recently. Transit spectroscopy of such planets with the James Webb Space Telescope (JWST) represents our best shot at obtaining the spectrum of a habitable planet within the next decade. As these planets are likely tidally locked, improved 3D numerical simulations of such planetary atmospheres are needed to guide target selection. Here we use a 3D climate system model, updated with new water-vapor absorption coefficients derived from the HITRAN 2012 database, to study ocean-covered planets at the inner edge of the HZ around late M to mid- K stars (2600 K less than or equal to T(sub eff) less than or equal to 4500 K). Our results indicate that these updated water-vapor coefficients result in significant warming compared to previous studies, so the inner HZ around M dwarfs is not as close as suggested by earlier work. Assuming synchronously rotating Earth-sized and Earth-mass planets with background 1 bar N2 atmospheres, we find that planets at the inner HZ of stars with T(sub eff) greater than 3000 K undergo the classical ''moist greenhouse'' (H2O mixing ratio greater than 10(exp -3) in the stratosphere) at significantly lower surface temperature (approx. 280 K) in our 3D model compared with 1D climate models (approx. 340 K). This implies that some planets around low-mass stars can simultaneously undergo water loss and remain habitable. However, for stars with T(sub eff) less than or equal to 3000 K, planets at the inner HZ may directly transition to a runaway state, while bypassing the moist greenhouse water loss entirely. We analyze transmission spectra of planets in a moist greenhouse regime and find that there are several prominent H2O features, including a broad feature between 5 and 8 micron, within JWST MIRI instrument range. Thus, relying only on standard Earth-analog spectra with 24 hr rotation period around M dwarfs for habitability studies will miss the strong H2O features that one would expect to see on synchronously rotating planets around M dwarf stars, with JWST.

Kopparapu, Ravi kumar↗

21st Century Trends in Antarctic Temperature and Polar Stratospheric Cloud (PSC) Area in the GEOS Chemistry-Climate Model

This study examines trends in Antarctic temperature and APSC, a temperature proxy for the area of polar stratospheric clouds, in an ensemble of Goddard Earth Observing System (GEOS) chemistry-climate model (CCM) simulations of the 21st century. A selection of greenhouse gas, ozone-depleting substance, and sea surface temperature scenarios is used to test the trend sensitivity to these parameters. One scenario is used to compare temperature trends in two versions of the GEOS CCM. An extended austral winter season is examined in detail. In May, June, and July, the expected future increase in CO2-related radiative cooling drives temperature trends in the Antarctic lower stratosphere. At 50 hPa, a 1.3 K cooling is expected between 2000 and 2100. Ozone levels increase, despite this robust cooling signal and the consequent increase in APSC, suggesting the enhancement of stratospheric transport in future. In the lower stratosphere, the choice of climate change scenarios does not affect the magnitude of the early winter cooling. Midwinter temperature trends are generally small. In October, APSC trends have the same sign as the prescribed halogen trends. That is, there are negative APSC trends in "grealistic future" simulations, where halogen loading decreases in accordance with the Montreal Protocol and CO2 continues to increase. In these simulations, the speed of ozone recovery is not influenced by either the choice of sea surface temperature and greenhouse gas scenarios or by the model version.

Hurwitz, M. M.↗

Systematic Benchmarking of Climate Models: Methodologies, Applications, and New Directions

As climate models become increasingly complex, there is a growing need to comprehensively and systematically assess model performance with respect to observations. Given the increasing number and diversity of climate model simulations in use, the community has moved beyond simple model intercomparison and toward developing methods capable of benchmarking a large number of simulations against a suite of climate metrics. Here, we present a detailed review of evaluation and benchmarking methods and approaches developed in the last decade, focusing primarily on scientific implications for Coupled Model Intercomparison Project (CMIP) simulations and CMIP6 results that contributed to the Intergovernmental Panel on Climate Change (IPCC) Sixth Assessment Report (AR6). Based on this review, we explain the resulting contemporary philosophy of model benchmarking, and provide clear distinctions and definitions of the terms model verification, process validation, evaluation, and benchmarking. While significant progress has been made in model development based on systematic evaluation and benchmarking efforts, some climate system biases still remain. The development of open‐source community software packages has played a fundamental role in identifying areas of significant model improvement and bias reduction. We review the key features of several software packages that have been commonly used over the past decade to evaluate and benchmark global and regional climate models. Additionally, we discuss best practices for the selection of evaluation and benchmarking metrics and for interpreting the obtained results, the importance of selecting suitable sources of reference data and accurate uncertainty quantification.

Environmental sciences↗

Evaluating Diurnal and Semi-Diurnal Cycle of Precipitation in CMIP6 Models Using Satellite- and Ground-Based Observations

We report the diurnal and semi-diurnal cycle of precipitation simulated from CMIP6 models during 1996-2005 are evaluated globally between 60°S and 60°N, as well as at ten selected locations representing three categories of diurnal cycle of precipitation: (1) afternoon precipitation over land, (2) early morning precipitation over ocean, and (3) nocturnal precipitation over land. Three satellite-based and two ground-based rainfall products are used to evaluate the climate models. Globally, the ensemble mean of CMIP6 models shows a diurnal phase of 3 to 4 hours earlier over land and 1 to 2 hours earlier over ocean, when compared with the latest satellite products. These biases are in line with what were found in previous versions of climate models but reduced compared to the CMIP5 ensemble mean. Analysis at the selected locations complimented with in-situ measurements further reinforces these results. Several CMIP6 models have shown a significant improvement in the diurnal cycle of precipitation compared to their CMIP5 counterparts, notably on delaying afternoon precipitation over land. This can be attributed to the use of more sophisticated convective parameterizations. Most models are still unable to capture the nocturnal peak associated with elevated convection and propagating mesoscale convective systems, with a few exceptions that allow convection to be initiated above the boundary layer to capture nocturnal elevated convection. We also quantify an encouraging consistency between the satellite- and ground-based precipitation measurements despite differing spatiotemporal resolutions and sampling periods, which provides confidence in using them to evaluate the diurnal and semi-diurnal cycle of precipitation in climate models.

58 GEOSCIENCES↗

The Atmospheric Chemistry and Climate Model Intercomparison Project (ACCMIP): Overview and Description of Models, Simulations and Climate Diagnostics

The Atmospheric Chemistry and Climate Model Intercomparison Project (ACCMIP) consists of a series of time slice experiments targeting the long-term changes in atmospheric composition between 1850 and 2100, with the goal of documenting composition changes and the associated radiative forcing. In this overview paper, we introduce the ACCMIP activity, the various simulations performed (with a requested set of 14) and the associated model output. The 16 ACCMIP models have a wide range of horizontal and vertical resolutions, vertical extent, chemistry schemes and interaction with radiation and clouds. While anthropogenic and biomass burning emissions were specified for all time slices in the ACCMIP protocol, it is found that the natural emissions are responsible for a significant range across models, mostly in the case of ozone precursors. The analysis of selected present-day climate diagnostics (precipitation, temperature, specific humidity and zonal wind) reveals biases consistent with state-of-the-art climate models. The model-to- model comparison of changes in temperature, specific humidity and zonal wind between 1850 and 2000 and between 2000 and 2100 indicates mostly consistent results. However, models that are clear outliers are different enough from the other models to significantly affect their simulation of atmospheric chemistry.

humidity↗

The optimization of model ensemble composition and size can enhance the robustness of crop yield projections

Linked climate and crop simulation models are widely used to assess the impact of climate change on agriculture. However, it is unclear how ensemble configurations (model composition and size) influence crop yield projections and uncertainty. Here, we investigate the influences of ensemble configurations on crop yield projections and modeling uncertainty from Global Gridded Crop Models and Global Climate Models under future climate change. We performed a cluster analysis to identify distinct groups of ensemble members based on their projected outcomes, revealing unique patterns in crop yield projections and corresponding uncertainty levels, particularly for wheat and soybean. Furthermore, our findings suggest that approximately six Global Gridded Crop Models and 10 Global Climate Models are sufficient to capture modeling uncertainty, while a cluster-based selection of 3-4 Global Gridded Crop Models effectively represents the full ensemble. The contribution of individual Global Gridded Crop Models to overall uncertainty varies depending on region and crop type, emphasizing the importance of considering the impact of specific models when selecting models for local-scale applications. Our results emphasize the importance of model composition and ensemble size in identifying the primary sources of uncertainty in crop yield projections, offering valuable guidance for optimizing ensemble configurations in climate-crop modeling studies tailored to specific applications.

Agriculture↗

Field validation of isotropic analytical models for simulating fabric shades

Fabric roller shades are common shading materials used in commercial and residential buildings. Accurately characterizing and modeling shades helps practitioners select the appropriate product and its control strategy based on climate and occupants' priorities, such as visual comfort and view to outdoors. Previous studies established a generalized method for modeling complex fenestration systems using data-driven tabulated bidirectional scattering distribution functions. However, deploying such a method at scale to all fabric shading products on the market is too costly and time-consuming. Analytical models that are based on a limited set of measurements (e.g., normal-normal and normal-hemispherical visible transmittance and reflectance, and directional cut-off angles) can be used to model the wide variety of shading products on the market. This study evaluates the performance of two isotropic analytical models, Roos-Wienold and Modified-Kotey, for modeling fabric roller shades, with a focus on the model's ability to predict occupant visual comfort. The performance evaluation was conducted through laboratory and field measurements and simulations. The results showed that both models are sufficient for predicting vertical illuminance at seated eye-level. Roos-Wienold model was able to predict binary visual comfort classification (glare/no-glare) under a wide range of luminance conditions, while Modified-Kotey model did not perform as well under high-contrast low-adaptation conditions. Both models are insufficient in predicting visual comfort at a four-point scale (e.g., imperceptible, perceptible, disturbing, intolerable). The two isotropic models become less accurate when the fabric exhibits high anisotropy.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

Field-scale dynamics of planting dates in the US Corn Belt from 2000 to 2020

Crop planting dates are a dynamic feature of agricultural systems that respond to short- and long-term climate signals, crop and cultivar selection, and technology changes. Planting date records are essential for yield gap analyses, accurate crop modeling, and tracking farmer adaptations to weather and climate change. Although planting dates have high variation at local scales due to heterogeneity in farm resources and decision-making, available long-term data on planting dates is largely restricted to aggregated regional statistics or, at best, satellite-derived datasets with limited spatiotemporal extent and at resolutions unable to distinguish individual fields (> 250 m). Here, we generated retrospective annual field-scale (30 m) planting date maps for both maize and soybeans spanning 2000-2020 across a 12 state region in the United States Corn Belt based on Landsat satellite data and a large ground sample of over 28,000 maize and soybean fields. Using training data from 2015-2020 for model selection, we found that planting date predictions improved with harmonic regression of Landsat data and additional annual weather covariates. The preferred random forests model approximately doubled performance compared to a null model based on state median planting dates, capturing 47% of field-level variation for maize (mean absolute error, MAE = 7.4 days) and 44% for soybeans (MAE = 7.5 days) against held-out ground truth test data for 2008-2014. We also evaluated the full 2000-2020 dataset with state agricultural statistics, finding strong agreement with median planting dates for maize (R 2 = 0.76, MAE = 4.4 days) and slightly lower agreement for soybeans (R 2 = 0.65, MAE = 5.4 days) when aggregated to the state level. We then used this new dataset to analyze environmental determinants of planting dates at a finer-scale than previously possible, controlling for unobserved variation at the sub-state district level. We found that during 2000-2020, each standard deviation increase in rainfall delayed planting by ~ 2.5 days, and fields with higher soil productivity ratings tended to be planted earlier. We did not find meaningful trends over the last two decades in planting dates for maize or soybeans, in contrast to trends towards earlier planting dates late last century and predicted for this period in climate adaptation studies. We hypothesize increases in early season rainfall may have inhibited these shifts towards earlier planting. Remotely sensed planting dates will be a useful tool for yield gap analyses, crop simulation modeling, and ongoing assessment of climate adaptation.

54 ENVIRONMENTAL SCIENCES↗

Evaluation and projection of long period return values of extreme daily precipitation in the CMIP5 and CMIP6 models

Using a non-stationary Generalized Extreme Value statistical method, we calculate selected extreme daily precipitation indices and their 20 year return values from the CMIP5 and CMIP6 climate models over the historical and future periods. We evaluate model performance of these indices and their return values in replicating similar quantities calculated from multiple gridded observational products. Difficulties in interpreting model quality in the context of observational uncertainties are discussed. Projections are framed in terms of specified global warming target temperatures rather than at specific times and under specific emissions scenarios. The change in framing shifts projection uncertainty due to differences in model climate sensitivity from the values of the projections to the timing of the global warming target. At their standard resolutions, we find there are no meaningful differences between the two generations of models in their quality or projections of simulated extreme daily precipitation.

54 ENVIRONMENTAL SCIENCES↗

Sensitivity of Riparian Buffer Designs to Climate Change—Nutrient and Sediment Loading to Streams: A Case Study in the Albemarle-Pamlico River Basins (USA) Using HAWQS

Riparian buffer zones (RBZs) provide multiple benefits to watershed ecosystems. We aimed to conduct an extensive sensitivity analysis of the RBZ designs to climate change nutrient and sediment loadings to streams. We designed 135 simulation scenarios starting with the six baselines RBZs (grass, urban, two-zone forest, three-zone forest, wildlife, and naturalized) in three 12-digit Hydrologic Unit Code watersheds within the Albemarle-Pamlico river basin (USA). Using the hydrologic and water quality system (HAWQS), we assessed the sensitivity of the designs to five water quality indicator (WQI) parameters: dissolved oxygen (DO), total phosphorous (TP), total nitrogen (TN), sediment (SD), and biochemical oxygen demand (BD). To understand the climate mitigation potential of RBZs, we identified a subset of future climate change projection models of air temperature and precipitation using EPA’s Locating and Selecting Scenarios Online tool. Analyses revealed optimal RBZ designs for the three watersheds. In terms of watershed ecosystem services sustainability, the optimal Urban RBZ in contemporary climate (1983–2018) reduced SD from 61–96%, TN from 34–55%, TP from 9–48%, and BD from 53–99%, and raised DO from 4–10% with respect to No-RBZ in the three watersheds. The late century’s (2070–2099) extreme mean annual climate changes significantly increased the projected SD and BD; however, the addition of urban RBZs was projected to offset the climate change reducing SD from 28–94% and BD from 69–93% in the watersheds. All other types of RBZs are also projected to fully mitigate the climate change impacts on WQI parameters except three-zone RBZ.

54 ENVIRONMENTAL SCIENCES↗

Climate change drives convergent evolution of root traits on Sky Island climate relicts

Roots are essential to the strategies plants use to survive in variable environments, yet we know little of how they vary within species. Experimental conditions demonstrate that intraspecific plant root traits respond strongly to variation in the environment; however, it is unclear when these responses can be characterized as evolution in response to selective pressures of climate change over many generations. Sky Islands are model, natural climate relict ecosystems to examine climate-change driven evolution. Utilizing a common garden with replicate genotypes of Populus angustifolia (Narrowleaf cottonwood) from six Sky Island (SI) populations and nine adjacent Mountain Chain (MC) populations across three genetic provenances, we hypothesized that SI root traits have diverged due to historical isolation in warmer, drier climates. When grown in common conditions, populations originating on SI’s showed convergent evolution across three distinct genetic provenances, which was characterized by 44.16% decreased total root length, 42.64% decreased average root volume, 43.31% decreased root surface area, and significantly less root trait variation, relative to adjacent mountain chains. Convergent evolution of root traits from trees originating on SI’s is correlated with changes in mean annual precipitation and potential evapotranspiration in the field over the past ~ 125 years. These results demonstrate a consistent pattern in root trait evolution at the landscape scale and the role of climate on the evolution of root traits in a genetic and geographic context relevant to climate change.

Convergent evolution↗

The Effect of Future Ambient Air Pollution on Human Premature Mortality to 2100 Using Output from the ACCMIP Model Ensemble

Ambient air pollution from ground-level ozone and fine particulate matter (PM(sub 2.5)) is associated with premature mortality. Future concentrations of these air pollutants will be driven by natural and anthropogenic emissions and by climate change. Using anthropogenic and biomass burning emissions projected in the four Representative Concentration Pathway scenarios (RCPs), the ACCMIP ensemble of chemistry climate models simulated future concentrations of ozone and PM(sub 2.5) at selected decades between 2000 and 2100. We use output from the ACCMIP ensemble, together with projections of future population and baseline mortality rates, to quantify the human premature mortality impacts of future ambient air pollution. Future air-pollution-related premature mortality in 2030, 2050 and 2100 is estimated for each scenario and for each model using a health impact function based on changes in concentrations of ozone and PM(sub 2.5) relative to 2000 and projected future population and baseline mortality rates. Additionally, the global mortality burden of ozone and PM(sub 2.5) in 2000 and each future period is estimated relative to 1850 concentrations, using present-day and future population and baseline mortality rates. The change in future ozone concentrations relative to 2000 is associated with excess global premature mortality in some scenarios/periods, particularly in RCP8.5 in 2100 (316 thousand deaths per year), likely driven by the large increase in methane emissions and by the net effect of climate change projected in this scenario, but it leads to considerable avoided premature mortality for the three other RCPs. However, the global mortality burden of ozone markedly increases from 382000 (121000 to 728000) deaths per year in 2000 to between 1.09 and 2.36 million deaths per year in 2100, across RCPs, mostly due to the effect of increases in population and baseline mortality rates. PM(sub 2.5) concentrations decrease relative to 2000 in all scenarios, due to projected reductions in emissions, and are associated with avoided premature mortality, particularly in 2100: between 2.39 and 1.31 million deaths per year for the four RCPs. The global mortality burden of PM(sub 2.5) is estimated to decrease from 1.70 (1.30 to 2.10) million deaths per year in 2000 to between 0.95 and 1.55 million deaths per year in 2100 for the four RCPs due to the combined effect of decreases in PM(sub 2.5) concentrations and changes in population and baseline mortality rates. Trends in future air-pollution-related mortality vary regionally across scenarios, reflecting assumptions for economic growth and air pollution control specific to each RCP and region. Mortality estimates differ among chemistry climate models due to differences in simulated pollutant concentrations, which is the greatest contributor to overall mortality uncertainty for most cases assessed here, supporting the use of model ensembles to characterize uncertainty. Increases in exposed population and baseline mortality rates of respiratory diseases magnify the impact on premature mortality of changes in future air pollutant concentrations and explain why the future global mortality burden of air pollution can exceed the current burden, even where air pollutant concentrations decrease.

Silva, Raquel A.↗

Sequential Decision Making (SDM) for Mesh Refinement and Model Selection in Multiscale, Multi-Physics Applications

Intelligent automation and decision support are needed to enhance computational efficiency and robustness in multiscale and multi-physics problems, including materials science, manufacturing, and climate and weather modeling. Current scientific computing approaches for enabling decisions by scientists fail to explore the role of learning, reasoning, and probabilistic planning. Often these decisions are not performed in real-time during the computation but are made prior to the start of the computation, which must be interrupted in order to make changes to the prior choices. Such interruptions at different stages of the computation increase the total computing time and the need for a human expert to frequently monitor the results. State of art scientific computing methods consist of rule-based algorithms that cannot automatically adapt to a dynamically changing computing environment. The development of a Sequential Decision Making (SDM) framework will automate scientific computing by optimizing the policies for mesh refinement, time-stepping, model and algorithm selection, resource allocation, and pre and post-processing. Our agent SDM framework for scientific computing will consist of data-driven learning (Classifier), automated reasoning (contextual knowledge), and probabilistic planning (Reinforcement Learning). In this project, we focused on three problems to demonstrate our SDM framework on a set of ordinary and partial differential equations. Classification of Lorenz system regions using Feed-Forward Neural Networks examined learning in the SDM framework. On the other hand, reasoning and planning in the SDM framework were used in two problems: adaptive time-stepping for nonlinear ODEs using on-policy RL algorithms, and adaptive mesh refinement for 2-D PDEs using off-policy RL algorithms.

97 MATHEMATICS AND COMPUTING↗