Search NASA⌕ Search

SEARCH · Search NASA

Results for “Model discrepancy”

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 181 records · Page 10
All sourcesNASA NTRSDOE OSTI
NASA NTRS2016

Black Carbon Absorption at the Global Scale Is Affected by Particle-Scale Diversity in Composition

Atmospheric black carbon (BC) exerts a strong, but uncertain, warming effect on the climate. BC that is coated with non-absorbing material absorbs more strongly than the same amount of BC in an uncoated particle, but the magnitude of this absorption enhancement (E(sub abs)) is not well constrained. Modelling studies and laboratory measurements have found stronger absorption enhancement than has been observed in the atmosphere. Here, using a particle-resolved aerosol model to simulate diverse BC populations, we show that absorption is overestimated by as much as a factor of two if diversity is neglected and population-averaged composition is assumed across all BC-containing particles. If, instead, composition diversity is resolved, we find E(sub abs) = 1 - 1.5 at low relative humidity, consistent with ambient observations. This study offers not only an explanation for the discrepancy between modelled and observed absorption enhancement, but also demonstrates how particle-scale simulations can be used to develop relationships for global-scale models.

Climate models↗
DOE OSTI2023

Bi-fidelity modeling of uncertain and partially unknown systems using DeepONets

Recent advances in modeling large-scale, complex physical systems have shifted research focuses towards data-driven techniques. However, generating datasets by simulating complex systems can require significant computational resources. Similarly, acquiring experimental datasets can prove difficult. For these systems, often computationally inexpensive, but in general inaccurate models, known as the low-fidelity models, are available. Here in this paper, we propose a bi-fidelity modeling approach for complex physical systems, where we model the discrepancy between the true system's response and a low-fidelity response in the presence of a small training dataset from the true system's response using a deep operator network, a neural network architecture suitable for approximating nonlinear operators. We apply the approach to systems that have parametric uncertainty and are partially unknown. Three numerical examples are used to show the efficacy of the proposed approach to model uncertain and partially unknown physical systems.

MATHEMATICS AND COMPUTING↗
DOE OSTI2025

Toward engineering lattice structures with the material point method (MPM)

This study examines the potential of two variants of the material point method—the generalized interpolation material point (GIMP) and dual domain material point (DDMP) methods—in developing a robust computational framework for engineering lattice structures under different loading conditions. The study begins with assessing the ability of the two methods in predicting elastic buckling phenomena using column geometries with and without initial geometric imperfections. The results indicate that both methods effectively capture buckling phenomena when initial geometric imperfections are introduced. After this verification step, we create several models of tetrahedral lattice structures with varying strut diameter and orientation and subject them to quasi-static loading. We then validate the numerical results using laboratory test results. The results show that, while both methods accurately predict load–displacement curves in the pre-buckling regime, their predictive capabilities diminish in the post-buckling regime. Through visual comparison between the numerical and experimental deformed shapes, it appears that the discrepancies between model and experimental results are attributed to initial geometric imperfections in the lattices that occurred during 3D printing. We then establish a second set of lattice models where different types of initial geometric imperfections are considered. The results from these models show that imperfections have a negligible influence in the pre-buckling regime but affect the behavior considerably in the post-buckling regime. As a final step in this work, we subject the lattice models to impact loading and employ hypothetical soft and stiff materials. These results show that the lattice stiffness, which depends on material stiffness, strut diameter, and orientation, significantly influences the ability of a lattice structure to resist impact. In particular, we find that a stiffer lattice (i.e., one made with a stiff material and thicker struts) is capable of absorbing more energy than a softer one during impact. Although material nonlinearities, inelasticity, and detailed contact formulations are not considered in this study, the findings obtained herein lay the groundwork for engineering lattice structures under extreme loading conditions through a simulation-driven framework based on particle-based methods.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2023

Conformer-Dependent Chemistry: Experimental Product Branching of the Vinyl Alcohol + OH + O 2 Reaction

The concentration of formic acid in Earth's tropo-sphere is underestimated by detailed chemical models compared to field observations. Phototautomerization of acetaldehyde to its less stable tautomer vinyl alcohol, followed by the OH-initiated oxidation of vinyl alcohol, has been proposed as a missing source of formic acid that improves the agreement between models and field measurements. Theoretical investigations of the OH + vinyl alcohol reaction in excess O 2 conclude that OH addition to the alpha carbon of vinyl alcohol produces formaldehyde + formic acid + OH, whereas OH addition to the β site leads to glycoaldehyde + HO 2 . Furthermore, these studies predict that the conformeric structure of vinyl alcohol controls the reaction pathway, with the anti-conformer of vinyl alcohol promoting α OH addition, whereas the syn-conformer promotes β addition. However, the two theoretical studies reach different conclusions regarding which set of products dominate. We studied this reaction using time-resolved multiplexed photoionization mass spectrometry to quantify the product branching fractions. Our results, supported by a detailed kinetic model, conclude that the glycoaldehyde product channel (arising mostly from syn-vinyl alcohol) dominates over formic acid production with a 3.6:1.0 branching ratio. Additionally, this result supports the conclusion of Lei et al. that conformer-dependent hydrogen bonding at the transition state for OH-addition controls the reaction outcome. As a result, tropospheric oxidation of vinyl alcohol creates less formic acid than recently thought, increasing again the discrepancy between models and field observations of Earth's formic acid budget.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2023

Machine learning for detection of 3D features using sparse x-ray tomographic reconstruction

In many inertial confinement fusion (ICF) experiments, the neutron yield and other parameters cannot be completely accounted for with one and two dimensional models. This discrepancy suggests that there are three dimensional effects that may be significant. Sources of these effects include defects in the shells and defects in shell interfaces, the fill tube of the capsule, and the joint feature in double shell targets. Due to their ability to penetrate materials, x rays are used to capture the internal structure of objects. Methods such as computational tomography use x-ray radiographs from hundreds of projections, in order to reconstruct a three dimensional model of the object. In experimental environments, such as the National Ignition Facility and Omega-60, the availability of these views is scarce, and in many cases only consists of a single line of sight. Mathematical reconstruction of a 3D object from sparse views is an ill-posed inverse problem. These types of problems are typically solved by utilizing prior information. Neural networks have been used for the task of 3D reconstruction as they are capable of encoding and leveraging this prior information. We utilize half a dozen, different convolutional neural networks to produce different 3D representations of ICF implosions from the experimental data. Deep supervision is utilized to train a neural network to produce high-resolution reconstructions. These representations are used to track 3D features of the capsules, such as the ablator, inner shell, and the joint between shell hemispheres. Machine learning, supplemented by different priors, is a promising method for 3D reconstructions in ICF and x-ray radiography, in general.

Wolfe, Bradley T. (ORCID:0000000268301614)↗
NASA NTRS2020

Bias in CMIP6 Models as Compared to Observed Regional Dimming and Brightening

Anthropogenic aerosol emissions have increased considerably over the last century, but climate effects and quantification of the emissions are highly uncertain as one goes back in time. This uncertainty is partly due to a lack of observations in the pre-satellite era, making the observations we do have before 1990 additionally valuable. Aerosols suspended in the atmosphere scatter and absorb incoming solar radiation and thereby alter the Earth's surface energy balance. Previous studies show that Earth system models (ESMs) do not adequately represent surface energy fluxes over the historical era. We investigated global and regional aerosol effects over the time period 1961–2014 by looking at surface downwelling shortwave radiation (SDSR). We used observations from ground stations as well as multiple experiments from eight ESMs participating in the Coupled Model Intercomparison Project Version 6 (CMIP6). Our results show that this subset of models reproduces the observed transient SDSR well in Europe but poorly in China. We suggest that this may be attributed to missing emissions of sulfur dioxide in China, sulfur dioxide being a precursor to sulfate, which is a highly reflective aerosol and responsible for more reflective clouds. The emissions of sulfur dioxide used in the models do not show a temporal pattern that could explain observed SDSR evolution over China. The results from various aerosol emission perturbation experiments from DAMIP, RFMIP and AerChemMIP show that only simulations containing anthropogenic aerosol emissions show dimming, even if the dimming is underestimated. Simulated clear-sky and all-sky SDSR do not differ greatly, suggesting that cloud cover changes are not a dominant cause of the biased SDSR evolution in the simulations. Therefore we suggest that the discrepancy between modeled and observed SDSR evolution is partly caused by erroneous aerosol and aerosol precursor emission inventories. This is an important finding as it may help interpret whether ESMs reproduce the historical climate evolution for the right or wrong reason.

surface downwelling shortwave radiation (SDSR)↗
DOE OSTI2022

Kinetic model development and Bayesian uncertainty quantification for the complete reduction of Fe-based oxygen carriers with CH 4 , CO, and H 2 for chemical looping combustion

In this work, three kinetic models are developed and calibrated for the complete multi-step reduction of an Fe-based oxygen carrier (OC) particle with CH 4 , CO, and H 2 , using data from thermogravimetric analysis. The complete reduction rate profiles exhibit complex dynamics whose trajectory is significantly different depending on the reducing gas. A Bayesian model building and parameter estimation framework is applied for simultaneous parameter and model structure uncertainty quantification. The final models show excellent agreement between model predictions and calibration data, as well as new data not used for calibration (for the reduction of the OC with CH 4 ). Parameter uncertainty is quantified by determining joint posterior distribution, and model structure uncertainty is addressed by incorporating Gaussian process stochastic functions (represented by Bayesian smoothing splines) into the kinetic models. The final kinetic models with discrepancy functions are readily employable in equation-oriented simulation and optimization platforms.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
NASA NTRS2001

Inorganic Chlorine Partitioning in the Summer Lower Stratosphere: Modeled and Measured [ClONO2/HCl] During POLARIS

We examine inorganic chlorine (Cly) partitioning in the summer lower stratosphere using in situ ER-2 aircraft observations made during the Photochemistry of Ozone Loss in the Arctic Region in Summer (POLARIS) campaign. New steady state and numerical models estimate [ClONO2]/[HCl] using currently accepted photochemistry. These models are tightly constrained by observations with OH (parameterized as a function of solar zenith angle) substituting for modeled HO2 chemistry. We find that inorganic chlorine photochemistry alone overestimates observed [ClONO2]/[HCl] by approximately 55-60% at mid and high latitudes. On the basis of POLARIS studies of the inorganic chlorine budget, [ClO]/[ClONO2], and an intercomparison with balloon observations, the most direct explanation for the model-measurement discrepancy in Cly partitioning is an error in the reactions, rate constants, and measured species concentrations linking HCl and ClO (simulated [ClO]/[HCl] too high) in combination with a possible systematic error in the ER-2 ClONO2 measurement (too low). The high precision of our simulation (+/-15% 1-sigma for [ClONO2]/[HCl], which is compared with observations) increases confidence in the observations, photolysis calculations, and laboratory rate constants. These results, along with other findings, should lead to improvements in both the accuracy and precision of stratospheric photochemical models.

Voss, P. B.↗
DOE OSTI2025

Assessing Clouds in GFDL's AM4.0 With Different Microphysical Parameterizations Using the Satellite Simulator Package COSP

We evaluate cloud simulations using satellite simulators against multiple observational data sets. These simulators have been run within the Geophysical Fluid Dynamics Laboratory's Atmosphere Model version 4.0 (AM4.0), as well as an alternative configuration where a fully two‐moment Morrison‐Gettelman cloud microphysical parameterization with prognostic precipitation (MG2) is applied, denoted as AM4‐MG2. The modeled cloud spatial distributions, vertical profiles, phase partitioning, cloud‐to‐precipitation transitions, and radiative effects compare reasonably well with satellite observations. Model biases include the under‐prediction of total and low‐level clouds, especially optically thin/intermediate clouds with cloud optical depth of less than 23, but the over‐prediction of thick clouds, indicating “too few, too bright” biases. These biases counteract each other, and give rise to reasonable estimates of cloud radiative effects. The underestimate of low‐level clouds is associated with too early and too frequent drizzle/precipitation formation. The precipitation bias is improved in AM4‐MG2, where the autoconversion scheme initiates the precipitation more realistically. There also exist discrepancies between models and observations for midlevel and high‐level clouds. Additional biases include the underestimate of liquid cloud fraction and the overestimate of ice cloud fraction.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2026

Non -degenerate marginal-likelihood calibration with application to quantum characterization

Here, we propose a marginal likelihood strategy within the Kennedy-O’Hagan (KOH) Bayesian framework, where a Gaussian process (GP) models the discrepancy between a physical system and its simulator. Our approach introduces a novel marginalized likelihood by integrating out the degenerate eigenspace of the covariance matrix, rather than approximating the original likelihood. Unlike approximation methods that compromise accuracy for computational efficiency, our method defines an exact likelihood—distinct from the original but preserving all relevant information. This formulation achieves computational efficiency and stability, even for large datasets where the covariance matrix nears degeneracy. Applied to the characterization of a superconducting quantum device at Lawrence Livermore National Laboratory, the approach enhances the predictive accuracy of the Lindblad master equations for modeling Ramsey measurement data by effectively quantifying uncertainties consistent with the quantum data.

general physics↗
DOE OSTI2023

Northern Hemisphere Snow Drought in Earth System Model Simulations and ERA5‐Land Data in 1980–2014

Abstract Low snow levels over the past few decades and predictions of a low‐to‐no snow future have spurred research into snow droughts, which pose a threat to water security and management. Systematic data‐model comparisons of snow drought have been lacking, hindering our understanding of the drivers of snow drought in the past. To address this gap, we analyzed snow drought events using standardized snow water equivalent index derived from monthly results of four numerical experiments using the E3SM Land Model (ELM) and ERA5‐Land data during the period of 1980–2014. Additionally, we compared snow drought duration calculated from models with those from the ERA5‐Land data during selected El Niño‐Southern Oscillation (ENSO) years. The numerical experiments were conducted with ELM driven by two prescribed atmospheric forcings, and with the coupled land‐atmosphere configuration of E3SM with and without plant hydraulics scheme feedback. Analysis reveals that 20%–30% of snow droughts occur due to factors other than above‐normal temperature and low snowfall, such as low soil moisture, warm soil temperature, and low relative humidity, etc., especially in high latitudes (50° North). Furthermore, our study highlights the exacerbating effect of ENSO events on snow drought conditions in various regions, despite some discrepancies between model and ERA5‐Land results. We also identified limitations of the coupled land‐atmosphere models in our current configuration in capturing the spatial patterns of snow droughts. This study underscores the challenge of predicting and mitigating snow drought and the need for a comprehensive understanding of the factors contributing to snow drought.

54 ENVIRONMENTAL SCIENCES↗
NASA NTRS1996

Responses of the Tropical Pacific to Wind Forcing as Observed by Spaceborne Sensors and Simulated by an Ocean General Circulation Model

In this study, satellite observations, in situ measurements, and model simulations are combined to assess the oceanic response to surface wind forcing in the equatorial Pacific. The surface wind fields derived from observations by the spaceborne special sensor microwave imager (SSM/I) and from the operational products of the European Centre for Medium-Range Weather Forecasts (ECMWF) are compared. When SSM/I winds are used to force a primitive-equation ocean general circulation model (OGCM), they produce 3 C more surface cooling than ECMWF winds for the eastern equatorial Pacific during the cool phase of an El Nino-Southern Oscillation event. The stronger cooling by SSM/I winds is in good agreement with measurements at the moored buoys and observations by the advanced very high resolution radiometer, indicating that SSM/I winds are superior to ECMWF winds in forcing the tropical ocean. In comparison with measurements from buoys, tide gauges, and the Geosat altimeter, the OGCM simulates the temporal variations of temperature, steric, and sea level changes with reasonable realism when forced with the satellite winds. There are discrepancies between model simulations and observations that are common to both wind forcing fields, one of which is the simulation of zonal currents; they could be attributed to model deficiencies. By examining model simulations under two winds, vertical heat advection and uplifting of the thermocline are found to be the dominant factors in the anomalous cooling of the ocean mixed layer.

Liu, W. Timothy↗
DOE OSTI2022

Peatland Loss in Southeast Asia Contributing to U.S. Biofuel’s Greenhouse Gas Emissions

Land use change (LUC) induced by biofuel production could lead to greenhouse gas (GHG) emissions, which potentially increase biofuel’s carbon intensity. Among the sources of LUC-related emissions for soy biodiesel, the contribution from peatland loss to agricultural plantations in Southeast Asia remains uncertain. Here, in this study, we analyzed LUC in Malaysia and Indonesia and modeled its impacts on the GHG emissions of soy biodiesel produced in the United States. It shows that oil palm plantations have more than doubled over 2001–2016 and the area of palm-on-peatlands (PoP) has expanded 3.7 times. Over new palm plantations, the share of PoP is about 19% regardless of time and location and the emission factor (EF) for peatland-to-palm conversion is estimated to be 41.5 Mg CO 2 ha –1 yr –1 . With these updates on PoP and EF, the contribution of peatland loss (0.7–5.1 g CO 2 e MJ –1 ) to biodiesel emissions is only 40–65% of previous estimates, which reduces discrepancies among model simulations used by different agencies. Based on emerging evidence on LUC and related carbon changes, our analysis reexamines regional peatland loss and its impacts on LUC emissions modeling and provides new insights into the estimation of LUC impacts on biofuels’ carbon intensity.

54 ENVIRONMENTAL SCIENCES↗
NASA NTRS2012

Longwave Band-by-band Cloud Radiative Effect and its Application in GCM Evaluation

The cloud radiative effect (CRE) of each longwave (LW) absorption band of a GCM fs radiation code is uniquely valuable for GCM evaluation because (1) comparing band-by-band CRE avoids the compensating biases in the broadband CRE comparison and (2) the fractional contribution of each band to the LW broadband CRE (f(sub CRE)) is sensitive to cloud top height but largely insensitive to cloud fraction, presenting thus a diagnostic metric to separate the two macroscopic properties of clouds. Recent studies led by the first author have established methods to derive such band ]by ]band quantities from collocated AIRS and CERES observations. We present here a study that compares the observed band-by-band CRE over the tropical oceans with those simulated by three different atmospheric GCMs (GFDL AM2, NASA GEOS-5, and CCCma CanAM4) forced by observed SST. The models agree with observation on the annual ]mean LW broadband CRE over the tropical oceans within +/-1W/sq m. However, the differences among these three GCMs in some bands can be as large as or even larger than +/-1W/sq m. Observed seasonal cycles of f(sub CRE) in major bands are shown to be consistent with the seasonal cycle of cloud top pressure for both the amplitude and the phase. However, while the three simulated seasonal cycles of f(sub CRE) agree with observations on the phase, the amplitudes are underestimated. Simulated interannual anomalies from GFDL AM2 and CCCma CanAM4 are in phase with observed anomalies. The spatial distribution of f(sub CRE) highlights the discrepancies between models and observation over the low-cloud regions and the compensating biases from different bands.

Huang, Xianglei↗
NASA NTRS2021

Mid-Pliocene Atlantic Meridional Overturning Circulation simulated in PlioMIP2

In the Pliocene Model Intercomparison Project Phase 2 (PlioMIP2), coupled climate models have been used to simulate an interglacial climate during the mid-Piacenzian warm period (mPWP; 3.264 to 3.025 Ma). Here, we compare the Atlantic Meridional Overturning Circulation (AMOC), poleward ocean heat transport and sea surface warming in the Atlantic simulated with these models. In PlioMIP2, all models simulate an intensified mid-Pliocene AMOC. However, there is no consistent response in the simulated Atlantic ocean heat transport nor in the depth of the Atlantic overturning cell. The models show a large spread in the simulated AMOC maximum, the Atlantic ocean heat transport and the surface warming in the North Atlantic. Although a few models simulate a surface warming of ~8–12°C in the North Atlantic, similar to the reconstruction from Pliocene Research, Interpretation and Synoptic Mapping (PRISM) version 4, most models appear to underestimate this warming. The large model spread and model–data discrepancies in the PlioMIP2 ensemble do not support the hypothesis that an intensification of the AMOC, together with an increase in northward ocean heat transport, is the dominant mechanism for the mid-Pliocene warm climate over the North Atlantic.

Pliocene Model Intercomparison Project Phase 2 (Pl↗
NASA NTRS2012

Carbon Consequences of Forest Disturbance and Recovery Across the Conterminous United States

Forests of North America are thought to constitute a significant long term sink for atmospheric carbon. The United States Forest Service Forest Inventory and Analysis (FIA) program has developed a large data base of stock changes derived from consecutive estimates of growing stock volume in the US. These data reveal a large and relatively stable increase in forest carbon stocks over the last two decades or more. The mechanisms underlying this national increase in forest stocks may include recovery of forests from past disturbances, net increases in forest area, and growth enhancement driven by climate or fertilization by CO2 and Nitrogen. Here we estimate the forest recovery component of the observed stock changes using FIA data on the age structure of US forests and carbon stocks as a function of age. The latter are used to parameterize forest disturbance and recovery processes in a carbon cycle model. We then apply resulting disturbance/recovery dynamics to landscapes and regions based on the forest age distributions. The analysis centers on 28 representative climate settings spread about forested regions of the conterminous US. We estimate carbon fluxes for each region and propagate uncertainties in calibration data through to the predicted fluxes. The largest recovery-driven carbon sinks are found in the South central, Pacific Northwest, and Pacific Southwest regions, with spatially averaged net ecosystem productivity (NEP) of about 100 g C / square m / a driven by forest age structure. Carbon sinks from recovery in the Northeast and Northern Lake States remain moderate to large owing to the legacy of historical clearing and relatively low modern disturbance rates from harvest and fire. At the continental scale, we find a conterminous U.S. forest NEP of only 0.16 Pg C/a from age structure in 2005, or only 0.047 Pg C/a of forest stock change after accounting for fire emissions and harvest transfers. Recent estimates of NEP derived from inventory stock change, harvest, and fire data show twice the NEP sink we derive from forest age distributions. We discuss possible reasons for the discrepancies including modeling errors and the possibility of climate and/or fertilization (CO2 or N) growth enhancements.

Williams, Christopher A.↗
DOE OSTI2023

Evaluation of liquid cloud albedo susceptibility in E3SM using coupled eastern North Atlantic surface and satellite retrievals

The impact of aerosol number concentration on cloud albedo is a persistent source of spread in global climate predictions due to multi-scale, interactive atmospheric processes that remain difficult to quantify. We use 5 years of geostationary satellite and surface retrievals at the US Department of Energy (DOE) Atmospheric Radiation Measurement (ARM) eastern North Atlantic (ENA) site in the Azores to evaluate the representation of liquid cloud albedo susceptibility for overcast cloud scenes in the DOE Energy Exascale Earth System Model version 1 (E3SMv1) and provide possible reasons for model–observation discrepancies. The overall distribution of surface 0.2 % CCN concentration values is reasonably simulated, but simulated liquid water path (LWP) is lower than observed and layer mean droplet concentration (N d ) comparisons are highly variable depending on the N d retrieval technique. E3SMv1's cloud albedo is greater than observed for given LWP and N d values due to a lower cloud effective radius than observed. However, the simulated albedo response to N d is suppressed due to a correlation between the solar zenith angle (SZA) and N d created by the seasonal cycle that is not observed. Controlling for this effect by examining the cloud optical depth (COD) shows that E3SMv1's COD response to CCN concentration is greater than observed. For surface-based retrievals, this is only true after controlling for cloud adiabaticity because E3SMv1's adiabaticities are much lower than observed. Assuming a constant adiabaticity in surface retrievals as done in top-of-atmosphere (TOA) retrievals narrows the retrieved ln N d distribution, which increases the cloud albedo sensitivity to ln N d to match the TOA sensitivity. The greater sensitivity of COD to CCN is caused by a greater Twomey effect in which the sensitivity of N d to CCN is greater than observed for TOA-retrieved N d , and once model–observation cloud adiabaticity differences are removed, this is also true for surface-retrieved N d . The LWP response to N d in E3SMv1 is overall negative as observed. Despite reproducing the observed LWP–N d relationship, observed clouds become much more adiabatic as N d increases, while E3SMv1 clouds do not, associated with more heavily precipitating clouds that are partially but not completely caused by deeper clouds and weaker inversions in E3SMv1. These cloud property differences indicate that the negative LWP–N d relationship is likely not caused by the same mechanisms in E3SMv1 and observations. The negative simulated LWP response also fails to mute the excessively strong Twomey effect, highlighting potentially important confounding factor effects that likely render the LWP–N d relationship non-causal. N d retrieval scales and assumptions, particularly related to cloud adiabaticity, contribute to substantial spreads in the model–observation comparisons, though enough consistency exists to suggest that aerosol activation, drizzle, and entrainment processes are critical areas to focus E3SMv1 development for improving the fidelity of aerosol–cloud interactions in E3SM.

54 ENVIRONMENTAL SCIENCES↗
NASA NTRS2006

Radiative Effect of Clouds on Tropospheric Chemistry in a Global Three-Dimensional Chemical Transport Model

Clouds exert an important influence on tropospheric photochemistry through modification of solar radiation that determines photolysis frequencies (J-values). We assess the radiative effect of clouds on photolysis frequencies and key oxidants in the troposphere with a global three-dimensional (3-D) chemical transport model (GEOS-CHEM) driven by assimilated meteorological observations from the Goddard Earth Observing System data assimilation system (GEOS DAS) at the NASA Global Modeling and Assimilation Office (GMAO). We focus on the year of 2001 with the GEOS-3 meteorological observations. Photolysis frequencies are calculated using the Fast-J radiative transfer algorithm. The GEOS-3 global cloud optical depth and cloud fraction are evaluated and generally consistent with the satellite retrieval products from the Moderate Resolution Imaging Spectroradiometer (MODIS) and the International Satellite Cloud Climatology Project (ISCCP). Results using the linear assumption, which assumes linear scaling of cloud optical depth with cloud fraction in a grid box, show global mean OH concentrations generally increase by less than 6% because of the radiative effect of clouds. The OH distribution shows much larger changes (with maximum decrease of approx.20% near the surface), reflecting the opposite effects of enhanced (weakened) photochemistry above (below) clouds. The global mean photolysis frequencies for J[O1D] and J[NO2] in the troposphere change by less than 5% because of clouds; global mean O3 concentrations in the troposphere increase by less than 5%. This study shows tropical upper tropospheric O3 to be less sensitive to the radiative effect of clouds than previously reported (approx.5% versus approx.20-30%). These results emphasize that the dominant effect of clouds is to influence the vertical redistribution of the intensity of photochemical activity while global average effects remain modest, again contrasting with previous studies. Differing vertical distributions of clouds may explain part, but not the majority, of these discrepancies between models. Using an approximate random overlap or a maximum-random overlap scheme to take account of the effect of cloud overlap in the vertical reduces the impact of clouds on photochemistry but does not significantly change our results with respect to the modest global average effect.

Liu, Hongyu↗