SEARCH · Engineering Papers
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.
Conference paper titled, Comparing microclimate weather conditions from the weather research and forecasting (WRF) model with typical meteorological year (TMY) data for building energy simulations accepted for Building Simulation 2025 Conference in Brisbane, Australia in Aug 2025.
To explore the integrated effects of ocean and ocean surface wave related physical processes on TC simulations, a set of three model simulations is performed. * In experiment 'A', the event is modeled using the Weather Research Forecasting (WRF) model alone with prescribed Sea Surface Temperature (SST) at 6-hour intervals. * In experiment ‘AO,’ WRF is coupled with the Finite Volume Community Ocean Model (FVCOM), enabling variable exchange between atmosphere and ocean, but without considering ocean surface wave-related physical processes. * In experiment ‘AOW’, WRF, FVCOM, and Simulating WAves Nearshore (SWAN) exchange variables with each other every hour through the OASIS3-MCT Coupler to allow direct and indirect atmosphere-ocean-wave interactions. * Observational data are also included in this dataset (Dropsonde, HRD-Radar, NDBC_wave). All simulations are initialized at 18:00 UTC on August 19, 2021, within a domain encompassing the western North Atlantic Ocean. The atmospheric domain features a horizontal resolution of 3 km. The ocean domain, which covers a substantial portion of the WRF ocean domain, employs an unstructured triangular grid with resolutions ranging from 3 km near the coast to 9 km in the open ocean, effectively resolving the complex coastline of the U.S. Northeast Coast. Initial and boundary conditions for the atmosphere model are obtained from the 6-hourly 0.25° NCEP (National Centers for Environmental Prediction) Global Forecast System (GFS; NCEP, 2015) data.
To explore the integrated effects of ocean and ocean surface wave related physical processes on TC simulations, a set of three model simulations is performed. * In experiment 'A', the event is modeled using the Weather Research Forecasting (WRF) model alone with prescribed Sea Surface Temperature (SST) at 6-hour intervals. * In experiment ‘AO,’ WRF is coupled with the Finite Volume Community Ocean Model (FVCOM), enabling variable exchange between atmosphere and ocean, but without considering ocean surface wave-related physical processes. * In experiment ‘AOW’, WRF, FVCOM, and Simulating WAves Nearshore (SWAN) exchange variables with each other every hour through the OASIS3-MCT Coupler to allow direct and indirect atmosphere-ocean-wave interactions. * Observational data are also included in this dataset (Dropsonde, HRD-Radar, NDBC_wave). All simulations are initialized at 18:00 UTC on August 19, 2021, within a domain encompassing the western North Atlantic Ocean. The atmospheric domain features a horizontal resolution of 3 km. The ocean domain, which covers a substantial portion of the WRF ocean domain, employs an unstructured triangular grid with resolutions ranging from 3 km near the coast to 9 km in the open ocean, effectively resolving the complex coastline of the U.S. Northeast Coast. Initial and boundary conditions for the atmosphere model are obtained from the 6-hourly 0.25° NCEP (National Centers for Environmental Prediction) Global Forecast System (GFS; NCEP, 2015) data.
To explore the integrated effects of ocean and ocean surface wave related physical processes on tropical cyclone simulations, a set of three model simulations is performed. * In experiment 'A', the event is modeled using the Weather Research Forecasting (WRF) model alone with prescribed Sea Surface Temperature (SST) at 6-hour intervals. * In experiment ‘AO,’ WRF is coupled with the Finite Volume Community Ocean Model (FVCOM), enabling variable exchange between atmosphere and ocean, but without considering ocean surface wave-related physical processes. * In experiment ‘AOW’, WRF, FVCOM, and Simulating WAves Nearshore (SWAN) exchange variables with each other every hour through the OASIS3-MCT Coupler to allow direct and indirect atmosphere-ocean-wave interactions. * Observational data are also included in this dataset (Dropsonde, HRD-Radar, NDBC_wave). All simulations are initialized at 18:00 UTC on August 19, 2021, within a domain encompassing the western North Atlantic Ocean. The atmospheric domain features a horizontal resolution of 3 km. The ocean domain, which covers a substantial portion of the WRF ocean domain, employs an unstructured triangular grid with resolutions ranging from 3 km near the coast to 9 km in the open ocean, effectively resolving the complex coastline of the U.S. Northeast Coast. Initial and boundary conditions for the atmosphere model are obtained from the 6-hourly 0.25° NCEP (National Centers for Environmental Prediction) Global Forecast System (GFS; NCEP, 2015) data. These CSV files are derived from the NetCDF files in the c0 dataset. Unlike the original format, where geographic coordinates were stored in a separate file, each CSV now embeds the corresponding latitude and longitude values alongside the measured or simulated variables.
Namibia is considered as one of the countries that are most vulnerable to climate change due to its generally dry climate and the percentage of its population that rely on subsistence agriculture for their livelihoods. Early warning systems are an important aspect of adapting to climate change. Weather forecasting relies on the use of numerical weather prediction models and these need to be configured properly. In this study we investigate the effects of using multi-nests and a convection scheme on the simulation of a heavy rainfall event over the Northwestern region of Kunene, Namibia. The event which was associated with a cut-off low system, was short-lived and resulted in over 45 mm of rainfall in one hour. For the multi-nest, a 9 km grid length parent domain is nested within the Global Forecast System (GFS) simulations, which in turn forces a 3 km grid spacing child domain. A different set of simulations are produced using a single nest of 3 km grid spacing, nested directly inside the GFS data. The simulations are produced with the convection scheme switched on and off. The impact of a single versus multi nest is found to be small in general, with slight differences in the location of high rainfall intensity. Switching off the convection schemes results in high rainfall intensity and increased detail in the simulations, including when a grid spacing of 9 km mis used. Using a grid spacing of 3 km with the convection scheme on, results in a loss of detail in the simulations as well as lower rainfall amounts. The study shows a need for different configurations to be tested before an optimum configuration can be selected for operational forecasting. We recommend further tests with different synoptic forcing to be conducted to identify a suitable configuration for Namibia.
Efforts to improve the prediction accuracy of high resolution (1-10 km) surface precipitation distribution and variability are of vital importance to local aspects of air pollution, wet deposition, and regional climate. However, precipitation biases and errors can occur at these spatial scales due to uncertainties in initial meteorological conditions and/or grid-scale cloud microphysics schemes. In particular, it is still unclear to what extent a subgrid-scale convection scheme could be modified to bring in scale-awareness for improving high-resolution short-term precipitation forecasts in the WRF model. To address these issues, we introduced scale-aware parameterized cloud dynamics for high-resolution forecasts by making several changes to the Kain-Fritsch (KF) convective parameterization scheme in the WRF model. These changes include subgrid-scale cloud-radiation interactions, a dynamic adjustment timescale, impacts of cloud updraft mass fluxes on grid-scale vertical velocity, and lifting condensation level-based entrainment methodology that includes scale dependency. A series of 48-hour retrospective forecasts using a combination of three treatments of convection (KF, updated KF, and the use of no cumulus parameterization), two cloud microphysics schemes and two types of initial condition datasets were performed over the U.S. southern Great Plains on 9- and 3-km grid spacings during the summers of 2002 and 2010. Results indicate that (1) the source of initial conditions play a key role in high-resolution precipitation forecasting, and (2) our updated KF scheme greatly alleviates the excessive precipitation at 9-km grid spacing and improves results at 3-km grid spacing as well. Overall, we found that the updated KF scheme incorporated into a high-resolution model does provide better forecasts for precipitation location and intensity.
This data package contains observation-based and model-simulated datasets (all provided in NetCDF format) for evaluating how wintertime lake-ice representation affects regional weather and climate over the Laurentian Great Lakes (freshwater lake ecosystem) during the high–ice-cover winter of 2009. The observational component includes: (1) Stage IV gridded precipitation at 4 km, hourly resolution for January–February 2009 over the Great Lakes region (radar–gauge multisensor precipitation analyses); (2) Great Lakes Surface Environmental Analysis (GLSEA) satellite-derived lake-ice coverage at 1.3 km, daily resolution for the 2009 winter months, providing ice coverage over Lakes Superior, Michigan, Huron, Erie, and Ontario; and (3) in situ measurements at the Standard Rock site on Lake Superior from the Great Lakes Evaporation Network (GLEN) at hourly resolution, including near-surface atmospheric variables and sensible and latent heat fluxes (air–lake exchange) at a fixed point location. The modeling component provides corresponding fields from two simulations, both archived at 4 km, hourly resolution: a standalone Weather Research Forecasting model (WRF) run driven by the European Centre for Medium-Range Weather Forecasts (ECMWF) Reanalysis 5 (ERA5), and a two-way coupled model using WRF and the Finite Volume Community Ocean Model (WRF-FVCOM, a 3-D hydrodynamic lake model). These outputs include variables relevant to air–lake interaction and lake-effect processes (e.g., near-surface temperature, humidity, wind, precipitation, and surface turbulent fluxes), enabling direct comparison with the observational datasets. Users can analyze and visualize these NetCDF files with common tools such as Python (e.g., xarray, netCDF4, numpy, pandas), NCO/CDO, Panoply, or ncview; NetCDF variables can also be converted to other formats (e.g., CSV, GeoTIFF) using these utilities.
Simulated present and future air quality is compared for the years 2006e2010 and 2048e2052 over the contiguous United States (CONUS) using the Community Multi-scale Air Quality (CMAQ) model. Regionally downscaled present and future climate results are developed using GISS and the Weather Research Forecasting (WRF) model. Present and future emissions are estimated using MARKAL 9R model. O3 and PM(sub 2.5) sensitivities to precursor emissions for the years 2010 and 2050 are calculated using CMAQDDM (Direct Decoupled Method). We find major improvements in future U.S. air quality including generally decreased MDA8 (maximum daily 8-hr average O3) mixing ratios and PM(sub 2.5) concentrations and reduced frequency of NAAQS O3 standard exceedances in most major U.S. cities. The Eastern and Pacific U.S. experience the largest reductions in summertime seasonal average MDA8 (up to 12 ppb) with localized decreases in the 4th highest MDA8 of the year, decreasing by up to 25 ppb. Results from a Climate Penalty (CP) scenario isolate the impact of climate change on air quality and show that future climate change tends to increase O3 mixing ratios in some regions of the U.S., with climate change causing increases of over 10 ppb in the annual 4th highest MDA8 in Los Angeles. Seasonal average PM(sub 2.5) decreases (2-4 microgram m(exp -3)) over the Eastern U.S. are accounted for by decreases in sulfate and nitrate concentrations resulting from reduced mobile and point source emissions of NO(sub x) and SO(sub x).
Abstract. Irrigated cultivation exerts a significant influence on the local climate and the hydrological cycle. The North China Plain (NCP) is known for its intricate agricultural system, marked by expansive cropland, high productivity, compact rotation, a semi-arid climate, and intensive irrigation practices. As a result, there has been considerable attention on the potential impact of this intensive irrigated agriculture on the local climate. However, studying the irrigation impact in this region has been challenging due to the lack of an accurate simulation of crop phenology and irrigation practices within the climate model. By incorporating double cropping with interactive irrigation, our study extends the capabilities of the Weather Research Forecast (WRF) model, which has previously demonstrated commendable performance in simulating single-cropping scenarios. This allows for two-way feedback between irrigated crops and climate, further enabling the inclusion of irrigation feedback from both ground and vegetation perspectives. The improved crop modeling system shows significant enhancement in capturing vegetation and irrigation patterns, which is evidenced by its ability to identify crop stages, estimate field biomass, predict crop yield, and project monthly leaf area index. The improved simulation of large-scale irrigated crops in the NCP can further enhance our understanding of the intricate relationship between agricultural development and climate change.
Recent GEWEX Cloud System Study (GCSS) model comparison projects have indicated that cloud-resolving models (CRMs) agree with observations better than traditional single-column models in simulating various types of clouds and cloud systems from different geographic locations. Current and future NASA satellite programs can provide cloud, precipitation, aerosol and other data at very fine spatial and temporal scales. It requires a coupled global circulation model (GCM) and cloud-scale model (termed a super-parameterization or multi-scale modeling framework, MMF) to use these satellite data to imiprove the understanding of the physical processes that are responsible for the variation in global and regional climate and hydrological systems. The use of a GCM will enable global coverage, and the use of a CRM will allow for better and more sophisticated physical parameterization. NASA satellite and field campaign cloud related datasets can provide initial conditions as well as validation for both the MMF and CRMs. A seed fund is available at NASA Goddard to build a MMF based on the 2D GCE model and the Goddard finite volume general circulation model (fvGCM). A prototype MMF will be developed by the end of 2004 and production runs will be conducted at the beginning of 2005. The purpose of this proposal is to augment the current Goddard MMF and other cloud modeling activities. I this talk, I will present: (1) A summary of the second Cloud Modeling Workshop took place at NASA Goddard, (2) A summary of the third TRMM Latent Heating Workshop took place at Nara Japan, (3) A brief discussion on the Goddard research plan of using Weather Research Forecast (WRF) model, and (4) A brief discussion on the GCE model on developing a global cloud simulator.
Weather forecasting over complex terrain with diverse land cover is challenging. Utilizing the high-resolution observations from New York State Mesonet (NYSM), we are able to evaluate the surface processes of the Weather Research Forecast (WRF) Model in a detailed, scale-dependent manner. Here, in the study, possible impacts of land–atmosphere interaction on surface meteorology and boundary layer cloud development are investigated with different model resolutions, land surface models (LSMs), and planetary boundary layer (PBL) physical parameterizations. The High-Resolution Rapid Refresh, version 3 (HRRR), forecasting model is used as a reference for the sensitivity evaluation. Results show that over complex terrain, the high-resolution simulations (1 km × 60 vertical levels) generally perform better compared to low-resolution (3 km × 50 levels) in both surface meteorology and cloud fields. LSMs play a more important role in surface meteorology compared to PBL schemes. The NoahMP land surface model exhibits daytime warmer and drier biases compared to the Rapid Update Cycle (RUC) due to better prediction of the Bowen ratio in RUC. The PBL schemes would affect the convective strength in the boundary layer. The Shin–Hong (SH) scale-aware scheme tends to produce the strongest convective strength in the PBL, while the ACM2 PBL scheme rarely resolved convection even at 1-km resolution. By considering the radiation effect of subgrid-scale (SGS) clouds, the Mellor–Yamada–Nakanishi–Niino eddy diffusivity mass flux (MYNN-EDMF) predicted the highest cloud coverage and lowest surface solar radiation bias. The configuration of SGS clouds in MYNN-EDMF would not only significantly reduce shortwave radiation bias, but also affect the convection behaviors through land surface–cloud–radiation interaction.
A modified Fast Spectral Bin Microphysics scheme (FSBM-2) embedded into the Weather Research and Forecasting Model (WRF) is used to simulate a mesoscale deep convective system observed during the Midlatitude Continental Convective Clouds Experiment (MC3E). FSBM-2 uses modified source codes as compared to the current FSBM (FSBM-1). In contrast to FSBM-1, FSBM-2 can simulate hail of several cm in diameter and includes additional processes such as spontaneous breakup of raindrops and aerosol regeneration by drop evaporation. It is shown that allowing large hail particles of diameters exceeding about 1 cm substantially increases the agreement between the simulated and observed squall-line structures both in the convective and stratiform regions. In contrast, if graupel particles are used to represent high-density hydrometeors in convective areas, the ratio of convective-to-stratiform areas diverges from the ratio seen in observations and maximum radar reflectivities are substantially underestimated. Analysis of snow size distributions in the stratiform area shows an important link between the ice crystals formed by homogeneous freezing in the convective area to ice particle number concentration in the stratiform region. Simulated raindrop size distributions in the stratiform area below the melting level from FSBM-2 show a good agreement with observations. The regeneration of CCN by droplet evaporation and detrainment of these CCN to the stratiform region increases the concentration of CCN there up to 8-9 km altitude; the additional CCN penetrate clouds and produce new droplets, leading to some convection intensification. Key words: Mesoscale convective systems (MCS), squall lines, spectral bin microphysics, large hail, in-cloud nucleation, size distributions of hydrometeors in MCS.
Abstract As risks from tropical cyclones (TCs) are fueled by climate change escalation, there is an urgent need for transformational solutions to complement traditional approaches. Seeding TCs using aerosols can be a promising method to reduce cyclone intensity, supported by theoretical understanding of the microphysical effects of aerosols on TC clouds. The ideal time to intervene effectively in TCs is likely during their initial stage, before TC wind speeds reach their peak. However, studies exploring potential aerosol effects on TC formation remain scarce. This study investigates how a TC embryo responds to the addition of aerosols of varying sizes using the Weather Research & Forecasting (WRF) model coupled with a spectral‐bin microphysics model. We found that aerosols of different sizes and concentrations distinctively affect the pre‐TC vortex's microstructure and dynamics. Fine and ultrafine aerosols enhance the latent heat of condensation, freezing, deposition, and riming, initially intensifying the vortex. However, this results in enhancement of the cold pool, thereby reducing inflow and surface fluxes, subsequently weakening the vortex. Coarse aerosols produce the opposite effect to that of fine and ultrafine aerosols. Coarse aerosols lead to a slower initial acceleration owing to enhanced warm rain. However, the resulting weaker cold pool is insufficient to effectively reduce the strength of the vortex at the later stage. This study provides critical insights into how aerosols of varying sizes and concentrations modulate the energy cascade and impact the evolution of a TC embryo, laying the groundwork for further research on TC risk management through aerosol intervention.
Numerical cloud resolving models (CRMs), which are based the non-hydrostatic equations of motion, have been extensively applied to cloud-scale and mesoscale processes during the past four decades. Recent GEWEX Cloud System Study (GCSS) model comparison projects have indicated that CRMs agree with observations in simulating various types of clouds and cloud systems from different geographic locations. Cloud resolving models now provide statistical information useful for developing more realistic physically based parameterizations for climate models and numerical weather prediction models. It is also expected that Numerical Weather Prediction (NWP) and regional scale model can be run in grid size similar to cloud resolving model through nesting technique. Current and future NASA satellite programs can provide cloud, precipitation, aerosol and other data at very fine spatial and temporal scales. It requires a coupled global circulation model (GCM) and cloud-scale model (termed a szrper-parameterization or multi-scale modeling -framework, MMF) to use these satellite data to improve the understanding of the physical processes that are responsible for the variation in global and regional climate and hydrological systems. The use of a GCM will enable global coverage, and the use of a CRM will allow for better and more sophisticated physical parameterization. NASA satellite and field campaign can provide initial conditions as well as validation through utilizing the Earth Satellite simulators. At Goddard, we have developed a multi-scale modeling system with unified physics. The modeling system consists a coupled GCM-CRM (or MMF); a state-of-the-art weather research forecast model (WRF) and a cloud-resolving model (Goddard Cumulus Ensemble model). In these models, the same microphysical schemes (2ICE, several 3ICE), radiation (including explicitly calculated cloud optical properties), and surface models are applied. In addition, a comprehensive unified Earth Satellite simulator has been developed at GSFC, which is designed to fully utilize the multi-scale modeling system. A brief review of the multi-scale modeling system with unified physics/simulator and examples is presented in this article.
One of the grand challenges of the Global Precipitation Measurement (GPM) mission is to improve precipitation measurements in mid- and high-latitudes during cold seasons through the use of high-frequency passive microwave radiometry. For this, the Weather Research Forecast (WRF) model with the Goddard microphysics scheme is coupled with the Satellite Data Simulation Unit (WRF-SDSU) that has been developed to facilitate the over-snowfall retrieval algorithm by providing virtual cloud library and microwave brightness temperature (To) measurements consistent to the GPM Microwave Imager (GMI). This study tested the Goddard cloud microphysics scheme in WRF in snowstorm events (January 20-22, 2007) over the Canadian CloudSat/CALIPSO Validation Project (C3VP) site up in Ontario, Canada. In this meeting, we will present the performance of the Goddard cloud microphysics scheme both with 2ice (ice and snow) and 3ice (ice, snow and graupel) as well as other WRF microphysics schemes. Results will be compared with the King Radar data. We will also use the WRF model outputs to drive the Goddard SDSU to calculate radiances and backscattering signals consistent to satellite direct observations. These simulated radiance are evaluated against the measurement from A-Train satellites.
The forthcoming Global Precipitation Measurement (GPM) Mission will provide next generation precipitation observations from a constellation of satellites. Since precipitation by nature has large variability and low predictability at cloud-resolving scales, the impact of precipitation data on the skills of mesoscale numerical weather prediction (NWP) is largely affected by the characterization of background and observation errors and the representation of nonlinear cloud/precipitation physics in an NWP data assimilation system. We present a data impact study on the assimilation of precipitation-affected microwave (MW) radiances from a pre-GPM satellite constellation using the Goddard WRF Ensemble Data Assimilation System (Goddard WRF-EDAS). A series of assimilation experiments are carried out in a Weather Research Forecast (WRF) model domain of 9 km resolution in western Europe. Sensitivities to observation error specifications, background error covariance estimated from ensemble forecasts with different ensemble sizes, and MW channel selections are examined through single-observation assimilation experiments. An empirical bias correction for precipitation-affected MW radiances is developed based on the statistics of radiance innovations in rainy areas. The data impact is assessed by full data assimilation cycling experiments for a storm event that occurred in France in September 2010. Results show that the assimilation of MW precipitation observations from a satellite constellation mimicking GPM has a positive impact on the accumulated rain forecasts verified with surface radar rain estimates. The case-study on a convective storm also reveals that the accuracy of ensemble-based background error covariance is limited by sampling errors and model errors such as precipitation displacement and unresolved convective scale instability.
The Weather Research and Forecasting (WRF) model is the next generation community mesoscale model designed to enhance collaboration between the research and operational sectors. The NM'S as a whole has begun a transition toward WRF as the mesoscale model of choice to use as a tool in making local forecasts. Currently, both the National Weather Service in Melbourne, FL (NWS MLB) and the Spaceflight Meteorology Group (SMG) are running the Advanced Regional Prediction System (AIRPS) Data Analysis System (ADAS) every 15 minutes over the Florida peninsula to produce high-resolution diagnostics supporting their daily operations. In addition, the NWS MLB and SMG have used ADAS to provide initial conditions for short-range forecasts from the ARPS numerical weather prediction (NWP) model. Both NM'S MLB and SMG have derived great benefit from the maturity of ADAS, and would like to use ADAS for providing initial conditions to WRF. In order to assist in this WRF transition effort, the Applied Meteorology Unit (AMU) was tasked to configure and implement an operational version of WRF that uses output from ADAS for the model initial conditions. Both agencies asked the AMU to develop a framework that allows the ADAS initial conditions to be incorporated into the WRF Environmental Modeling System (EMS) software. Developed by the NM'S Science Operations Officer (S00) Science and Training Resource Center (STRC), the EMS is a complete, full physics, NWP package that incorporates dynamical cores from both the National Center for Atmospheric Research's Advanced Research WRF (ARW) and the National Centers for Environmental Prediction's Non-Hydrostatic Mesoscale Model (NMM) into a single end-to-end forecasting system. The EMS performs nearly all pre- and postprocessing and can be run automatically to obtain external grid data for WRF boundary conditions, run the model, and convert the data into a format that can be readily viewed within the Advanced Weather Interactive Processing System. The EMS has also incorporated the WRF Standard Initialization (SI) graphical user interface (GUT), which allows the user to set up the domain, dynamical core, resolution, etc., with ease. In addition to the SI GUT, the EMS contains a number of configuration files with extensive documentation to help the user select the appropriate input parameters for model physics schemes, integration timesteps, etc. Therefore, because of its streamlined capability, it is quite advantageous to configure ADAS to provide initial condition data to the EMS software. One of the biggest potential benefits of configuring ADAS for ingest into the EMS is that the analyses could be used to initialize either the ARW or NMM. Currently, the ARPS/ADAS software has a conversion routine only for the ARW dynamical core. However, since the NIvIM runs about 2.5 times faster than the ARW, it is quite advantageous to be able to run an ADAS/NMM configuration operationally due to the increased efficiency.
The Weather Research and Forecasting (WRF) model can be used to diagnose regional land-atmosphere (L-A) coupling strength in the absence of sufficient observations but subjected to uncertainties associated with model physical parameterizations. In this study, we propose a framework to quantify and reduce model physical parameterization uncertainties associated with surface fluxes and L-A coupling. An ensemble of WRF simulations with different physical schemes is used to simulate surface fluxes and land-atmosphere coupling strength over the Amazon region. The physical parameterizations investigated include cloud microphysics (MP), land surface processes (LSM), planetary boundary layer (PBL), surface layer (SL), and cumulus (CU). We perform 120 ensemble simulations using the WRF model and different combinations of six MPs, three LSMs, six PBLs and SLs and three CUs. The measurements from the GoAMAZON field campaign and satellite data are used to evaluate model performance. A Multi-way analysis of variance (ANOVA) approach is applied to quantify the relative importance of different physics processes on L-A coupling. The Tukey’s test is used to sort schemes that have no significant differences into one group. The suite of physics that result in the best simulations of the corresponding variables are selected based on the Taylor skill score. Results show that the relative importance of processes and their interaction vary with the variables of interest. For example, CU was the most important process in modulating soil moisture, 2 m-humidity, latent heat, and net radiation. LSM showed dominant effects on 2 m-temperature and also has the largest impact on sensible heat and the lifting condensation level. The best physical parameterization ensembles show much narrower ranges of the variables of interest than the priori ensemble. Results of this study show the roles of different physical processes in modulating L-A interactions, quantify model uncertainties from physical processes, and provide insights for improving the model physics parameterizations.