Search NASA⌕ Search

SEARCH · Search NASA

Results for “Earth 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 307 records · Page 17
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Evaluating long-term emission impacts of large-scale electric vehicle deployment in the US using a human-Earth systems model

While large-scale adoption of electric vehicles (EVs) globally would reduce carbon dioxide (CO 2 ) and traditional air pollutant emissions from the transportation sector, emissions from the electric sector, refineries, and potentially other sources would change in response. Here, a multi-sector human-Earth systems model is used to evaluate the net long-term emission implications of large-scale EV adoption in the US over widely differing pathways of the evolution of the electric sector. Our results indicate that high EV adoption would decrease net CO 2 emissions through 2050, even for a scenario where all electric sector capacity additions through 2050 are fossil fuel technologies. Greater net CO 2 reductions would be realized for scenarios that emphasize renewables or decarbonization of electricity production. Net air pollutant emission changes in 2050 are relatively small compared to expected overall decreases from recent levels to 2050. States participating in the Regional Greenhouse Gas Initiative experience greater CO 2 and air pollutant reductions on a percentage basis. These results suggest that coordinated, multi-sector planning can greatly enhance the climate and environmental benefits of EVs. Additional factors are identified that influence the net emission impacts of EVs, including the retirement of coal capacity, refinery operations under reduced gasoline demands, and price-induced fuel switching in residential heating and in the industrial sector.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

ORBIT: Oak Ridge Base Foundation Model for Earth System Predictability

Earth system predictability is challenged by the complexity of environmental dynamics and the multitude of variables involved. Current AI foundation models, although advanced by leveraging large and heterogeneous data, are often constrained by their size and data integration, limiting their effectiveness in addressing the full range of Earth system prediction challenges. To overcome these limitations, we introduce the Oak Ridge Base Foundation Model for Earth System Predictability (ORBIT), an advanced vision transformer model that scales up to 113 billion parameters using a novel hybrid tensor-data orthogonal parallelism technique. As the largest model of its kind, ORBIT surpasses the current climate AI foundation model size by a thousandfold. Performance scaling tests conducted on the Frontier supercomputer have demonstrated that ORBIT achieves 684 petaFLOPS to 1.6 exaFLOPS sustained throughput, with scaling efficiency maintained at 41% to 85% across 49,152 AMD GPUs. These breakthroughs establish new advances in AI-driven climate modeling and demonstrate promise to significantly improve the Earth system predictability.

Wang, Xiao↗
DOE OSTI2025

Regime-based aerosol–cloud interactions from CALIPSO-MODIS and the Energy Exascale Earth System Model version 2 (E3SMv2) over the Eastern North Atlantic

This study investigates aerosol-cloud interactions in marine boundary layer (MBL) clouds using an advanced deep-learning-driven synoptic-regime-based framework, combining satellite data (CALIPSO vertically resolved aerosol extinction and MODIS cloud properties) with 1° nudged Energy Exascale Earth System Model version 2 (E3SMv2) simulation over the Eastern North Atlantic (ENA; ∼10°×10°, 2006–2014). The E3SMv2 captures observed seasonal variations in cloud droplet number concentrations (N d ) and liquid water path (LWP), though it systematically underestimates N d . We then partition ENA meteorology into four synoptic regimes (Pre-Trough, Post-Trough, Ridge, Trough) via a deep-learning clustering of ERA5 reanalysis fields, enabling regime-dependent aerosol-cloud interactions analyses. Both satellite and E3SMv2 exhibit an inverted-V LWP-N d relationship. In Post-Trough and Ridge regimes, the satellite shows stronger negative LWP-N d sensitivities than in Pre-Trough regime. The Trough regime displays a muted satellite LWP response. In comparison, the model predicts more exaggerated LWP responses across regimes, with LWP increasing too quickly at low N d and decreasing more sharply at high N d , especially in Pre-Trough and Trough regimes. These exaggerated model LWP sensitivities may stem from uncertainties in representing drizzle processes, entrainment, and turbulent mixing. As for N d susceptibility to aerosols, N d increases with MBL aerosol extinction in both datasets, but the simulated aerosol-cloud interactions appear oversensitive to meteorological conditions. Overall, E3SMv2 better captures aerosol effects under regimes that favor stratiform clouds (Post-Trough, Ridge), but performance deteriorates for regimes with deeper, dynamically complex clouds (Trough), highlighting the need for improved representations of those cloud processes in climate models.

Environmental sciences↗
DOE OSTI2022

Benthic Kinetic Energy in an Earth System Model and Observations

The benthic layer of the ocean is closely correlated with every upper layer above it. Close analysis of benthic layer stress can yield an accurate depiction of the state of the ocean, as well as the atmosphere above it. Through the use of the Department of Energy’s Energy Exascale Earth System Model (E3SM), we took a look at how these features changed over time in a simulation spanning the historical period from 1958 to 2018 in the Arctic Ocean north of 60° N. Making use of E3SM’s regional refinement capability to place higher resolution in regions of interest, the E3SM-Arctic configuration has ~10 km resolution in the ocean and sea ice components in this region, coarsening up to ~60 km outside of the Arctic (Veneziani et al. 2022). In particular, we narrowly focused on the benthic kinetic energy levels, and how trends changed over time. To validate the fidelity of the E3SM model, we compared model current speed outputs over a 369-day period in the Bering Strait to actual recorded observations during this time period. Our analysis of the benthic kinetic energy mean, standard deviation, and trend throughout this time period can be used to characterize the effects of environmental impacts in the Arctic Ocean.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Characterizing the Variation and Covariation of Cloud Microphysical Properties and Implications for Simulation of Subgrid-scale Warm-Rain Processes in Earth System Models (Final DOE-ASR Report)

Warm marine boundary layer (MBL) clouds constitute an important component in the global climate system, and precipitation plays a central role in controlling the water budget, radiative effects, and lifetime of these MBL clouds. Unfortunately, because of the relatively coarse effective grid resolution of the current generation of Earth system models (ESMs), the variety of cloud microphysical processes occurring inside an ESM grid cell are often oversimplified or unconstrained by observations. For example, the warm rain processes (e.g., autoconversion and accretion) are usually parameterized as nonlinear functions of grid-mean cloud properties. Because of the nonlinear nature of these functions, neglecting variability within the ESM grid volume can lead to substantial biases in precipitation production, cloud cover, and surface radiative fluxes. In state-of-the art ESMs, the influence of subgrid-scale variability is represented as an enhancement factor (EF) coefficient to the autoconversion, and accretion rates calculated from the model variables. However, EF is typically taken to be a constant or even used as a knob to tune model cloud properties to match observations, an ad hoc approach that may yield a desired cloud outcome yet introduce compensating errors. In this project, we used the combination of in situ cloud microphysics measurements from the ACE-ENA field campaign and large-eddy simulations (LES) to characterize and understand subgrid-scale variations and co-variations of cloud microphysical properties and use the results to evaluate and improve the representation of subgrid warm-rain processes in ESMs, in particular the EF used to tune the autoconversion and accretion processes. In this final report, we summarize our research activities and main findings in Section 2, provide a list of publications (Section 3) and presentations (Section 4) resulted from our research, and briefly discuss the student activities supported by this project.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Assessing the Atmospheric Response to Subgrid Surface Heterogeneity in the Single-Column Community Earth System Model, Version 2 (CESM2)

Land-atmosphere interactions are central to the evolution of the atmospheric boundary layer and the subsequent formation of clouds and precipitation. Existing global climate models represent these connections with bulk approximations on coarse spatial scales, but observations suggest that small-scale variations in surface characteristics and co-located turbulent and momentum fluxes can significantly impact the atmosphere. Recent model development efforts have attempted to capture this phenomenon by coupling existing representations of subgrid-scale (SGS) heterogeneity between land and atmosphere models. Such approaches are in their infancy and it is not yet clear if they can produce a realistic atmospheric response to surface heterogeneity. Here, we implement a parameterization to capture the effects of SGS heterogeneity in the Community Earth System Model (CESM2), and compare single-column simulations against high-resolution Weather Research and Forecasting (WRF) large-eddy simulations (LESs), which we use as a proxy for observations. The CESM2 experiments increase the temperature and humidity variances in the lowest atmospheric levels, but the response is weaker than in WRF-LES. In part, this is attributed to an underestimate of surface heterogeneity in the land model due to a lack of SGS meteorology, a separation between deep and shallow convection schemes in the atmosphere, and a lack of explicitly represented mesoscale secondary circulations. These results highlight the complex processes involved in capturing the effects of SGS heterogeneity and suggest the need for parameterizations that communicate their influence not only at the surface but also vertically.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Climate response to off-equatorial stratospheric sulfur injections in three Earth system models – Part 1: Experimental protocols and surface changes

There is now substantial literature on climate model studies of equatorial or tropical stratospheric SO 2 injections that aim to counteract the surface warming produced by rising concentrations of greenhouse gases. Here we present the results from the first systematic intercomparison of climate responses in three Earth system models wherein the injection of SO 2 occurs at different latitudes in the lower stratosphere: CESM2-WACCM6, UKESM1.0 and GISS-E2.1-G. The first two use a modal aerosol microphysics scheme, while two versions of GISS-E2.1-G use a bulk aerosol (One-Moment Aerosol, OMA) and a two-moment (Multiconfiguration Aerosol TRacker of mIXing state, MATRIX) microphysics approach, respectively. Our aim in this work is to determine commonalities and differences between the climate model responses in terms of the distribution of the optically reflective sulfate aerosols produced from the oxidation of SO 2 and in terms of the surface response to the resulting reduction in solar radiation. A focus on understanding the contribution of characteristics of models transport alongside their microphysical and chemical schemes, and on evaluating the resulting stratospheric responses in different models, is given in the companion paper (Bednarz et al., 2023). The goal of this exercise is not to evaluate these single-point injection simulations as stand-alone proposed strategies to counteract global warming; instead we determine sources and areas of agreement and uncertainty in the simulated responses and, ultimately, the possibility of designing a comprehensive intervention strategy capable of managing multiple simultaneous climate goals through the combination of different injection locations. We find large disagreements between GISS-E2.1-G and the CESM2-WACCM6 and UKESM1.0 models regarding the magnitude of cooling per unit of aerosol optical depth (AOD) produced, which varies from 4.7 K per unit of AOD in CESM2-WACCM6 to 16.7 K in the GISS-E2.1-G version with two-moment aerosol microphysics. By normalizing the results with the global mean response in each of the models and thus assuming that the amount of SO 2 injected is a free parameter that can be managed independently, we highlight some commonalities in the overall distributions of the aerosols, in the inter-hemispheric surface temperature response and in shifts to the Intertropical Convergence Zone, as well as some areas of disagreement, such as the extent of the aerosol confinement in the equatorial region and the efficiency of the transport to polar latitudes. In conclusion, we demonstrate that it is possible to use these simulations to produce more comprehensive injection strategies in multiple climate models. However, large differences in the injection magnitudes can be expected, potentially increasing inter-model spreads in some stratospheric quantities (such as aerosol distribution) while reducing the spread in the surface response in terms of temperature and precipitation; furthermore, the selection of the injection locations may be dependent on the models' specific stratospheric transport.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Diurnal rainfall response to the physiological and radiative effects of CO 2 in tropical forests in the Energy exascale earth system model v1

Understanding how the connection between rainfall and tropical forests will respond to increasing CO$_2$ concentrations is a key element in understanding how the tropical water cycle will respond to increasing CO$_2$. The plant physiological and radiative impacts of CO$_2$ on rainfall patterns over tropical forest regions are examined in the Energy Exascale Earth System Model version 1.1 (E3SMv1.1-BGC) biogeochemistry experiments. Composite analysis reveals a dampening of the diurnal cycle of rainfall over the Amazon, Congo, and Maritime Continent in response to rising CO$_2$ levels, regardless of the sign of total rainfall change. A full factorial model experiment confirms that the CO$_2$ radiative and CO$_2$ plant physiological effects can individually or jointly reduce the magnitude of the rainfall diurnal cycle, though the physical pathway giving rise to the reduction differs between the two effects. For the physiological response, stomatal closure reduces evapotranspiration, which dries the boundary layer and raises the lifting condensation level. These effects combine to reduce deep convective rainfall during its peak occurrence in the late daytime to early nighttime period. For the radiative response, a relative reduction in daytime CAPE (consistent with a reduction in the diurnal temperature range) leads to less frequent triggering of deep convection and a reduction of rainfall diurnal amplitude. These diurnal rainfall changes are structurally similar across seasons, and show little sensitivity to representation of nutrient coupling for the land biogeochemistry. In agreement with previous findings, the physiological response has only minor impact on extreme rainfall relative to the radiative response.

54 ENVIRONMENTAL SCIENCES↗
NASA NTRSDate not supplied

Language Model For Earth Science: Exploring Potential Downstream Applications As Well As Current Challenges

The use of deep learning techniques to build transformer language models such as SciBERT and GPT3 have transformed the natural language technology (NLT) landscape. These new NLTs are being used in speech to text and vice versa, auto-mated text classification, sentiment analysis, topic modeling, text summarization, and cognitive assistants. While Earth science has no shortage of unstructured data such as journal and conference papers, little efforts have focused on harnessing NLTs for knowledge extraction and supporting the scientific process. This paper surveys the use of language models in different science. BERT-E, a new Earth science-specific language model, is presented. BERT-E is generated using a transfer learning solution. A language model that has already been trained for general Science (SciBERT) is fine-tuned using abstracts and full text extracted from various Earth science-related articles. A downstream keywords classification application is used for evaluation, and the use of BERT-E shows improved performance. The need to develop a robust set of benchmarks in evaluating the language model such as BERT-E is discussed. Finally, example applications are presented to inspire additional ideas for applications using domain-specific language models.

R Ramachandran↗
DOE OSTI2021

On the Generation of Weddell Sea Polynyas in a High-Resolution Earth System Model

Larger Weddell Sea polynyas (WSPs), differentiated in this study from the smaller Maud Rise Polynyas (MRPs), forming to the east of the prime meridian in the proximity of the Maud Rise seamount, have last been observed in the 1970s. We investigate WSPs that grow realistically out of MRPs in a high-resolution (HR) preindustrial simulation with the Energy Exascale Earth System Model version 0.1. The formation of MRPs requires HR to simulate the detailed flow around Maud Rise, while the realistic formation of WSPs requires a model to produce MRPs. Furthermore, WSPs tend to follow periods of a prolonged build-up of a heat reservoir at depth and weakly negative wind-stress curl in association with the core of the southern hemisphere westerlies at an anomalously northern position. While this scenario also leads to drier conditions over the central Weddell Sea, which some literature claims to be a necessary condition for the formation of WSPs, our model results indicate that open-ocean polynyas do not occur during periods of weakly negative wind-stress curl despite drier atmospheric conditions. Our study supports the hypothesis noted in earlier studies that a shift from a weakly negative to a strongly negative wind-stress curl over the Weddell Sea is a prerequisite for WSPs to form, together with a large heat reservoir at depth. However, the ultimate trigger is a pronounced MRP; whose associated convection creates high surface salinity anomalies that propagate westward with the flow of the Weddell Gyre. If large enough, these anomalies trigger the formation of a WSP and a pulse of newly formed Antarctic Bottom Water.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Future Changes in Active and Inactive Atlantic Hurricane Seasons in the Energy Exascale Earth System Model

Abstract North Atlantic tropical cyclones (TCs) have considerable interannual variability, with La Niña and the positive phase of the Atlantic Meridional Mode (AMM) tending to drive active hurricane seasons, and El Niño and the negative AMM often driving inactive seasons. Here, we analyze how active and inactive Atlantic hurricane seasons may change in the future using the high resolution Energy Exascale Earth System Model (E3SM). We performed atmosphere‐only simulations forced by sea‐surface temperature patterns characteristic of La Niña and the positive AMM jointly, and El Niño and the negative AMM jointly, in historical and future climates. Projected Atlantic TCs become more frequent in the future by approximately 34% during El Niño and negative AMM and by 66% during La Niña and positive AMM, with a significant increase in the portion of intense TCs. Warmer SSTs increase TC potential intensity, with reduced wind shear and increased mid‐tropospheric humidity further supporting TC activity.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Representing Global Soil Erosion and Sediment Flux in Earth System Models

Abstract Soil erosion produces enormous amounts of sediment, carbon, and nutrient fluxes from land to rivers, thus playing crucial roles in global biogeochemical cycles and food security. To predict soil erosion in the context of climate and land use changes, we explicitly parameterize cropland management actions (i.e., irrigation, conserved agriculture, and crop residue management) and geological factors (i.e., lithology and glacier) in the Energy Exascale Earth System Model (E3SM) soil erosion module. The erosion model is calibrated using a global‐scale regionalized parameter calibration method. The spatial variabilities of the modeled and the Revised Universal Soil Loss Equation (RUSLE) based benchmark soil erosion are consistent across vegetation, climate and soil properties. Compared with independent data, our model shows a bias reduction in 59% of the observations relative to the RUSLE‐based soil erosion, with 53% of the bias reduction exceeding 50%. This improvement is mainly due to a better representation of the topographic effect on soil erosion. Our results indicate that conserved agriculture practices have effectively reduced soil erosion in cropland by over 25% in the United States and Argentina. In contrast, irrigation has increased soil erosion in many Asian countries. For upland sediment flux, our model is consistent with the WBMsed benchmark data in inter‐basin variability but could be more skillful in simulating intra‐basin variability because it couples soil erosion and sediment flux explicitly. The developed model provides useful skills for more realistic predictions of soil erosion and river sediment dynamics under environmental changes.

58 GEOSCIENCES↗
DOE OSTI2021

Dataset and scripts for manuscript "Using Neural Network Ensembles to Separate Ocean Biogeochemical and Physical Drivers of Phytoplankton Biogeography in Earth System Models"

Please note: The title of this version contains an updated title for the manuscript compared to the previous version of this dataset. This is only due to title updates during the peer review process for the manuscript. The zip file contains the scripts, functions, and source files for the manuscript titled "Using Neural Network Ensembles to Separate Ocean Biogeochemical and Physical Drivers of Phytoplankton Biogeography in Earth System Models." The manuscript has been submitted for peer review. Please consult the README file for information on the specifications of the files. These files may occasionally be updated to add annotations to the scripts to make them more user friendly and to correct any errors.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Clusters of Regional Precipitation Seasonality Change in the Community Earth System Model, Version 2

Abstract The likely changes to precipitation seasonality with warming are both impactful and not well understood. This work aims to describe areas that experience similar changes to seasonal precipitation irrespective of the original underlying precipitation seasonality. We train a self-organizing map on the difference between the seasonal cycle of precipitation in the past and in a high-warming future climate as represented by the Community Earth System Model, version 2, to create regions with similar changes in precipitation seasonality. This method is applied separately over land and ocean surfaces because of the differing processes leading to precipitation over each. This method indicates that future changes in seasonal precipitation are most varied in the tropics because of a southward shift in the intertropical convergence zone. The seasonal shifts found over midlatitude oceans indicate a poleward shift in atmospheric river activity. We find a correspondence between certain land-based precipitation changes and Köppen climate classification. The seasonality of large-scale and convective precipitation is examined for each region. The relationship between the seasonal changes to precipitation and associated atmospheric processes is discussed. These processes include atmospheric rivers, the intertropical convergence zone, tropical cyclones, and monsoons.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Investigating Future Arctic Sea Ice Loss and Near–Surface Wind Speed Changes Related to Surface Roughness Using the Community Earth System Model

The Arctic is undergoing a pronounced and rapid transformation in response to changing greenhouse gasses, including reduction in sea ice extent and thickness. There are also projected increases in near-surface Arctic wind. This study addresses how the winds trends may be driven by changing surface roughness and/or stability in different Arctic regions and seasons, something that has not yet been thoroughly investigated. We analyze 50 experiments from the Community Earth System Model Version 2 (CESM2) Large Ensemble and five experiments using CESM2 with an artificially decreased sea ice roughness to match that of the open ocean. We find that with a smoother surface there are higher mean wind speeds and slower mean ice speeds in the autumn, winter, and spring. The artificially reduced surface roughness also strongly impacts the wind speed trends in autumn and winter, and we find that atmospheric stability changes are also important contributors to driving wind trends in both experiments. In contrast to the clear impacts on winds, the sea ice mean state and trends are statistically indistinguishable, suggesting that near-surface winds are not major drivers of Arctic sea ice loss. Two major results of this work are: (a) the near-surface wind trends are driven by changes in both surface roughness and near-surface atmospheric stability that are themselves changing from sea ice loss, and (b) the sea ice mean state and trends are driven by the overall warming trend due to increasing greenhouse gas emissions and not significantly impacted by coupled feedbacks with the surface winds.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

An overview of cloud–radiation denial experiments for the Energy Exascale Earth System Model version 1

Abstract. The interaction between clouds and radiation is a key process within the climate system, and assessing the impacts of that interaction provides valuable insights into both the present-day climate and future projections. Many modeling experiments have been designed over the years to probe the impact of the cloud radiative effect (CRE) on the climate, including those that seek to disrupt the mean CRE effect and those that only disrupt the covariance of the CRE with the circulation. Seven such experimental designs have been added to the Energy Exascale Earth System Model version 1 (E3SMv1) of the US Department of Energy. These experiments include both the first and second iterations of the Clouds On/Off Klimate Intercomparison Experiment (COOKIE) experimental design, as well as the cloud-locking method. This paper documents the code changes necessary to implement such experiments and also provides detailed instructions for how to run them. Analyses across experiment types provide valuable insights and confirm the findings of prior studies, including the role of cloud radiative heating toward intensifying the monsoon, intensifying rain rates, and poleward expansion of the general circulation owing to cloud feedbacks.

Harrop, Bryce E. (ORCID:0000000339524525)↗
NASA NTRSDate not supplied

Coupled Aerosol-Chemistry Simulations of the January 2022 Eruptions of Hunga Tonga-Hunga Ha’apai in the NASA GEOS Earth System Model

The January 2022 eruptions of the underwater Hunga Tonga-Hunga Ha’apai volcano injected more than 100 Tg of water vapor and approximately 0.5 Tg of sulfur dioxide into the stratosphere. Injected materials reached as high as ~55 km in altitude, while the main plume from the eruption travelled west over Australia and the Indian Ocean between about 20 – 30 km altitude. We investigate the transport and the impact of the erupted materials in coupled aerosol-chemistry simulations performed with the NASA Goddard Earth Observing System (GEOS) Earth system model. A sulfur mechanism introduced into the Global Modeling Initiative (GMI) stratospheric-tropospheric chemistry package allows for an interactive simulation of the water vapor-chemistry impacts, and the large water vapor perturbation results in rapid conversion of sulfur dioxide to sulfate aerosol in our aerosol mechanism. We report on the results of a multi-year ensemble of simulations performed with this system that include a control ensemble (no eruption), a water vapor-only injection ensemble, and a water vapor and sulfur dioxide injection ensemble. Simulations of the near-field aerosol and chemistry transport are compared to available measurements from OMPS-LP, OMPS-NM, CALIOP, MLS, and other sensors. The extended impact of the volcanic materials on the stratospheric composition and chemistry over the next several years is further investigated in forecast simulations.

Peter Colarco↗
DOE OSTI2021

An efficient method to estimate the probability density of seismic Green's functions

We present a computationally efficient method to approximate the probability distribution of seismic Green's functions given the uncertainty of an Earth model. The method is based on the Karhunen-Loève (KL) theorem and an approximation of the Green's function (or seismogram) covariance. Using Monte Carlo (MC) simulations as a control case, we demonstrate that our KL-based method can accurately reproduce a probability distribution of seismograms that results from an uncertain Earth model for a MC-derived seismogram covariance. We then describe a method to estimate the covariance of the seismograms resulting from those Earth models that is not based on MC simulations. We use the estimated Green's function covariance in conjunction with our KL-based method to produce a Green's function probability distribution, and compare that distribution to a Green's function probability distribution produced using a MC finite difference method. We find that the Green's function probability distribution approximated using our KL-based method generally mimics that produced using the MC simulations, especially for direct-arriving body waves. However the accuracy of the KL-based method generally decreases for later times in the simulated Green's function distribution.

58 GEOSCIENCES↗