Search NASA⌕ Search

SEARCH · Search NASA

Results for “Atmospheric Model Parameterizations”

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 685 records · Page 38

Leveraging Scientific Community Knowledge for Air Quality Model Chemistry Parameterizations

Air pollution contributes to adverse health outcomes. Approximately 121 million people in the United States—one third of the population—live where National Ambient Air Quality Standards (NAAQS) are violated. In most cases, the criteria pollutants exceeding standards are ozone (O 3 ) and fine particles (PM 2.5 ). In addition, 188 substances known or suspected to cause cancer or other serious health effects are designated as hazardous air pollutants (HAPs). Essentially, all O 3 and significant portions of PM 2.5 and HAPs are produced in the atmosphere through chemical and physical processes. In the case of PM 2.5 , subcomponents formed primarily from precursor gases—sulfate, nitrate, ammonium, and secondary organic aerosol (SOA)—account for 60% of the U.S. county-level annual mean concentration. In addition, 47% of the cancer risk and 25% of the noncancer risk from HAPs have been attributed to atmospheric chemistry rather than direct emissions. In this article, we introduce the role of chemical mechanisms in air quality models, a new atmospheric science community effort, and needs for further mechanism development.

Havala O.T. Pye↗

A scheme for computing surface layer turbulent fluxes from mean flow surface observations

A physical model and computational scheme are developed for generating turbulent surface stress, sensible heat flux and humidity flux from mean velocity, temperature and humidity at some fixed height in the atmospheric surface layer, where conditions at this reference level are presumed known from observations or the evolving state of a numerical atmospheric circulation model. The method is based on coupling the Monin-Obukov surface layer similarity profiles which include buoyant stability effects on mean velocity, temperature and humidity to a force-restore formulation for the evolution of surface soil temperature to yield the local values of shear stress, heat flux and surface temperature. A self-contained formulation is presented including parameterizations for solar and infrared radiant fluxes at the surface. Additional parameters needed to implement the scheme are the thermal heat capacity of the soil per unit surface area, surface aerodynamic roughness, latitude, solar declination, surface albedo, surface emissivity and atmospheric transmissivity to solar radiation.

Hoffert, M. I.↗

The contribution of Humboldt Glacier, northern Greenland, to sea-level rise through 2100 constrained by recent observations of speedup and retreat

Abstract. Humboldt Glacier, northern Greenland, has retreated and accelerated through the 21st century, raising concerns that it could be a significant contributor to future sea-level rise. We use a data-constrained ensemble of three-dimensional higher-order ice sheet model simulations to estimate the likely range of sea-level rise from the continued retreat of Humboldt Glacier. We first solve for basal traction using observed ice thickness, bed topography, and ice surface velocity from the year 2007 in a PDE-constrained (partial differential equation) optimization. Next, we impose calving rates to match mean observed retreat rates from winter 2007–2008 to winter 2017–2018 in a transient calibration of the exponent in the power-law basal friction relationship. We find that power-law exponents in the range of 1/7–1/5 – rather than the commonly used 1/3–1 – are necessary to reproduce the observed speedup over this period. We then tune an iceberg calving parameterization based on the von Mises stress yield criterion in another transient-calibration step to approximate both observed ice velocities and terminus position in 2017–2018. Finally, we use the range of basal friction relationship exponents and calving parameter values to generate the ensemble of model simulations from 2007–2100 under three climate forcing scenarios from CMIP5 (two RCP8.5 forcings, Representative Concentration Pathway) and CMIP6 (one SSP5-8.5 forcing, Shared Socioeconomic Pathway). Our simulations predict 5.2–8.7 mm of sea-level rise from Humboldt Glacier, significantly higher than a previous estimate (∼ 3.5 mm) and equivalent to a substantial fraction of the 40–140 mm predicted by ISMIP6 from the whole Greenland Ice Sheet. Our larger future sea-level rise prediction results from the transient calibration of our basal friction law to match the observed speedup, which requires a semi-plastic bed rheology. In many simulations, our model predicts the growth of a sizable ice shelf in the middle of the 21st century. Thus, atmospheric warming could lead to more retreat than predicted here if increased surface melt promotes hydrofracture of the ice shelf. Our data-constrained simulations of Humboldt Glacier underscore the sensitivity of model predictions of Greenland outlet glacier response to warming to choices of basal shear stress and iceberg calving parameterizations. Further, transient calibration of these parameterizations, which has not typically been performed, is necessary to reproduce observed behavior. Current estimates of future sea-level rise from the Greenland Ice Sheet could, therefore, contain significant biases.

54 ENVIRONMENTAL SCIENCES↗

Assessment of Warm and Dry Bias over ARM SGP Site in E3SMv2 and E3SM-MMF

Abstract Many climate models exhibit a dry and warm bias over the central United States during the summer months, including the Energy Exascale Earth System Model (E3SM) and its Multiscale Modeling Framework (MMF) configuration. Understanding the causes of this bias is important to shine a light on this common model error and reduce the uncertainty in future projections. In this study, we use E3SMv2 and E3SM-MMF to assess how parameterized and resolved convection affect temperature and precipitation biases over the Southern Great Plains site of the Atmospheric Radiation Measurement program. Both configurations overestimate near-surface temperature and underestimate precipitation at the ARM SGP site. The bias is associated with a lack of low-level clouds during days without precipitation and too much incoming solar radiation causing the surface to warm. Low-level cloud fraction in E3SM-MMF during the nonprecipitating days is lower in comparison to E3SMv2 and observation, consistent with the larger warm bias. We also find that the underestimated precipitation can be characterized as “too frequent, too weak” in E3SMv2 and “too rare, too intense” in E3SM-MMF. These deficiencies conspire to sustain the warm and dry bias over the central United States.

54 ENVIRONMENTAL SCIENCES↗

Building an AI-enhanced modeling framework to address multiscale predictability challenges

Focal Area(s): Build an AI-enhanced modeling framework that integrates the three focus areas in the solicitation. Science Challenge: The 4M-2N complexities (Multiscale, Multiphysics, Multibody, Multidimension, Non-linearity, and Non-Gaussianality) of atmospheric aerosol-cloud-precipitation-turbulence-radiation system poses physical and computational challenges to further advance predictive models; We plan to address the challenges by developing an AI-enhanced modeling framework that facilitates automated calibration and improvement of subgrid parameterizations, enhance data assimilation of measurements to improve initial and boundary conditions used to drive the physical model, and optimally blends data-driven and physics-based forecasting models.

54 ENVIRONMENTAL SCIENCES↗

Large-scale offshore wind farm effects on weather and climate in Puerto Rico (Final Technical Report)

Puerto Rico’s current electricity generation heavily relies on imported fossil fuels. This results in an average cost of electricity higher than the U.S. mainland average in all sectors (residential, commercial, and industrial), despite abundant local offshore wind resources, which have the potential to provide secure, low-cost energy generation and consequent economic prosperity. However, effects on atmospheric and oceanic circulation resulting from large-scale deployments of offshore wind farms have not been previously studied at tropical latitude. This project addressed this knowledge gap through a computational modeling effort designed to capture the coupled dynamics of the atmosphere and the ocean in presence of offshore wind farms. Results indicate that wind farm wakes can alter wind stress, generate Ekman-driven vertical transport, and potentially affect nutrient distribution. While full model coupling remains challenging, progress in parameterization and large-eddy simulations provides a foundation for future research. The project contributes to DOE’s Earth System modeling efforts and supports STEM workforce development.

17 WIND ENERGY↗

Atmospheric Radiation Measurement (ARM)/CART site

The DOE's Atmospheric Radiation Measurement (ARM) goals are as follows: (1) to provide an experimental test bed for improving the treatment of radiative transfer in global climate models (GCM's) under all kinds of cloud cover; and (2) to improve the parameterization and modeling of cloud formation, maintenance, dissipation, and related processes in GCM's. The scientific requirements which are most critical to the objectives of ARM are discussed.

Pennell, Bill↗

Physically-Derived Dynamical Cores in Atmospheric General Circulation Models

The algorithm chosen to represent the advection in atmospheric models is often used as the primary attribute to classify the model. Meteorological models are generally classified as spectral or grid point, with the term grid point implying discretization using finite differences. These traditional approaches have a number of shortcomings that render them non-physical. That is, they provide approximate solutions to the conservation equations that do not obey the fundamental laws of physics. The most commonly discussed shortcomings are overshoots and undershoots which manifest themselves most overtly in the constituent continuity equation. For this reason many climate models have special algorithms to model water vapor advection. This talk focuses on the development of an atmospheric general circulation model which uses a consistent physically-based advection algorithm in all aspects of the model formulation. The shallow-water model of Lin and Rood (QJRMS, 1997) is generalized to three dimensions and combined with the physics parameterizations of NCAR's Community Climate Model. The scientific motivation for the development is to increase the integrity of the underlying fluid dynamics so that the physics terms can be more effectively isolated, examined, and improved. The expected benefits of the new model are discussed and results from the initial integrations will be presented.

Rood, Richard B.↗

The Mars Crustal Magnetic Field Control of Plasma Boundary Locations and Atmospheric Loss: MHD Prediction and Comparison with MAVEN

We present results from a global Mars time-dependent MHD simulation under constant solar wind and solar radiation impact considering inherent magnetic field variations due to continuous planetary rotation. We calculate the 3-D shapes and locations of the bow shock (BS) and the induced magnetospheric boundary (IMB) and then examine their dynamic changes with time. We develop a physics-based, empirical algorithm to effectively summarize the multidimensional crustal field distribution. It is found that by organizing the model results using this new approach, the Mars crustal field shows a clear, significant influence on both the IMB and the BS. Specifically, quantitative relationships have been established between the field distribution and the mean boundary distances and the cross-section areas in the terminator plane for both of the boundaries. The model-predicted relationships are further verified by the observations from the NASA Mars Atmosphere and Volatile EvolutioN (MAVEN) mission. Our analysis shows that the boundaries are collectively affected by the global crustal field distribution, which, however, cannot be simply parameterized by a local parameter like the widely used subsolar longitude. Our calculations show that the variability of the intrinsic crustal field distribution in Mars-centered Solar Orbital itself may account for approx.60% of the variation in total atmospheric loss, when external drivers are static. It is found that the crustal field has not only a shielding effect for atmospheric loss but also an escape-fostering effect by positively affecting the transterminator ion flow cross-section area.

Global↗

More Realistic Intermediate Depth Dry Firn Densification in the Energy Exascale Earth System Model (E3SM)

Abstract Earth system models account for seasonal snow cover, but many do not accommodate the deeper snowpack on ice sheets (aka firn) that slowly transforms to ice under accumulating snowfall. To accommodate and resolve firn depths of up to 60 m in the Energy Exascale Earth System Model's land surface model (ELM), we add 11 layers to its snowpack and evaluate three dry snow compaction equations in multi‐century simulations. After comparing results from ELM simulations (forced with atmospheric reanalysis) with empirical data, we find that implementing into ELM a two‐stage firn densification model produces more accurate dry firn densities at intermediate depths of 20–60 m. Compared to modeling firn using the equations in the (12 layer) Community Land Model (version 5), switching to the two‐stage firn densification model (with 16 layers) significantly decreases root‐mean‐square errors in upper 60 m dry firn densities by an average of 41 kg m −3 (31%). Simulations with three different firn density parameterizations show that the two‐stage firn densification model should be used for applications that prioritize accurate upper 60 m firn air content (FAC) in regions where the mean annual surface temperature is greater than roughly −31°C. Because snow metamorphism, firn density, and FAC are major components in modeling ice sheet surface albedo, melt water retention, and climatic mass balance, these developments advance broader efforts to simulate the response of land ice to atmospheric forcing in Earth system models.

58 GEOSCIENCES↗

A comparison of the impact of two time-differencing schemes on the NASA-GLAS climate model

Evidence is presented of the sensitivity of a GCM to the time-differencing scheme employed when the physical parameterizations and space discretization are not changed. Two-time marching schemes, the leapfrog and the Matsuno schemes, are analyzed and tested on the NASA-Goddard Laboratory for Atmospheric Studies fourth-order GCM in terms of the behavior and stability of two-month-averaged fields. Linear analysis indicates that Rossby waves are slightly accelerated and slightly damped when the Matsuno scheme is utilized and that these effects are scale selective, being smallest for the longest waves.

Pfeffer, Richard L.↗

Characteristics of the Martian atmosphere surface layer

Elements of various terrestrial boundary layer models are extended to Mars in order to estimate sensible heat, latent heat, and momentum fluxes within the Martian atmospheric surface ('constant flux') layer. The atmospheric surface layer consists of an interfacial sublayer immediately adjacent to the ground and an overlying fully turbulent surface sublayer where wind-shear production of turbulence dominates buoyancy production. Within the interfacial sublayer, sensible and latent heat are transported by non-steady molecular diffusion into small-scale eddies which intermittently burst through this zone. Both the thickness of the interfacial sublayer and the characteristics of the turbulent eddies penetrating through it depend on whether airflow is aerodynamically smooth or aerodynamically rough, as determined by the Roughness Reynold's number. Within the overlying surface sublayer, similarity theory can be used to express the mean vertical windspeed, temperature, and water vapor profiles in terms of a single parameter, the Monin-Obukhov stability parameter. To estimate the molecular viscosity and thermal conductivity of a CO2-H2O gas mixture under Martian conditions, parameterizations were developed using data from the TPRC Data Series and the first-order Chapman-Cowling expressions; the required collision integrals were approximated using the Lenard-Jones potential. Parameterizations for specific heat and binary diffusivity were also determined. The Brutsart model for sensible and latent heat transport within the interfacial sublayer for both aerodynamically smooth and rough airflow was experimentally tested under similar conditions, validating its application to Martian conditions. For the surface sublayer, the definition of the Monin-Obukhov length was modified to properly account for the buoyancy forces arising from water vapor gradients in the Martian atmospheric boundary layer. It was found that under most Martian conditions, the interfacial and surface sublayers offer roughly comparable resistance to sensible heat and water vapor transport and are thus both important in determining the associated fluxes.

Clow, G. D.↗

Marine Shallow Cloud Adjustments to the Presence of Shortwave-Absorbing Aerosols: Advancing Understanding Through a Combined Analysis of Lasic Datasets and Process Modeling (Final Report)

This is a final report for this award, whose purpose was to support the analysis of data emanating from the LASIC (Layered Atlantic Smoke Interaction with Clouds) campaign. (Layered Atlantic Smoke Interactions with Clouds) is a strategy to improve our understanding of aged carbonaceous aerosol, its seasonal evolution, and the mechanisms by which clouds adjust to the presence of the aerosol. The observational strategy centers on deploying the AMF1 cloud, aerosol, and atmospheric profiling instrumentation to Ascension Island, located within the trade-wind shallow cumulus regime (14.50W, 80S) 3000 km offshore of continental Africa. The location is within the latitude zone of the maximum outflow of aerosol, with the deepening boundary layer known to entrain free-tropospheric smoke. The primary activities for LASIC are: 1) to improve current knowledge on aged biomass burning aerosol and its radiative properties as a function of the seasonal cycle; 2) to use surface-based remote sensing to sensitively interrogate the atmosphere for the relative vertical location of aerosol and clouds; 3) to improve our understanding of the cloud adjustments to the presence of shortwave-absorbing aerosol within the vertical column, both through aerosol-radiation and through aerosol-cloud interactions; 4) to aid low cloud parameterization efforts for climate models. This award serves to further these activities and includes a sub-award to Dr. Pablo Saide at UCLA for activities 3) and 4).

54 ENVIRONMENTAL SCIENCES↗

Climate and hydraulic traits interact to set thresholds for liana viability

Lianas, or woody vines, and trees dominate the canopy of tropical forests and comprise the majority of tropical aboveground carbon storage. These growth forms respond differently to contemporary variation in climate and resource availability, but their responses to future climate change are poorly understood because there are very few predictive ecosystem models representing lianas. We compile a database of liana functional traits (846 species) and use it to parameterize a mechanistic model of liana-tree competition. The substantial difference between liana and tree hydraulic conductivity represents a critical source of inter-growth form variation. Here, we show that lianas are many times more sensitive to drying atmospheric conditions than trees as a result of this trait difference. Further, we use our competition model and projections of tropical hydroclimate based on Representative Concentration Pathway 4.5 to show that lianas are more susceptible to reaching a hydraulic threshold for viability by 2100.

54 ENVIRONMENTAL SCIENCES↗

Light absorption enhancement of black carbon in a pyrocumulonimbus cloud

Abstract Pyrocumulonimbus (pyroCb) firestorm systems have been shown to inject significant amounts of black carbon (BC) to the stratosphere with a residence time of several months. Injected BC warms the local stratospheric air, consequently perturbing transport and hence spatial distributions of ozone and water vapor. A distinguishing feature of BC-containing particles residing within pyroCb smoke is their thick surface coatings made of condensed organic matter. When coated with non-refractory materials, BC’s absorption is enhanced, yet the absorption enhancement factor ( E abs ) for pyroCb BC is not well constrained. Here, we perform particle-scale measurements of BC mass, morphology, and coating thickness from inside a pyroCb cloud and quantify E abs using an established particle-resolved BC optics model. We find that the population-averaged E abs for BC asymptotes to 2.0 with increasing coating thickness. This value denotes the upper limit of E abs for thickly coated BC in the atmosphere. Our results provide observationally constrained parameterizations of BC absorption for improved radiative transfer calculations of pyroCb events.

54 ENVIRONMENTAL SCIENCES↗

ISMIP6-based projections of ocean-forced Antarctic Ice Sheet evolution using the Community Ice Sheet Model

Abstract. The future retreat rate for marine-based regions of the Antarctic Ice Sheet is one of the largest uncertainties in sea-level projections. The Ice Sheet Model Intercomparison Project for CMIP6 (ISMIP6) aims to improve projections and quantify uncertainties by running an ensemble of ice sheet models with atmosphere and ocean forcing derived from global climate models. Here, the Community Ice Sheet Model (CISM) is used to run ISMIP6-based projections of ocean-forced Antarctic Ice Sheet evolution. Using multiple combinations of sub-ice-shelf melt parameterizations and calibrations, CISM is spun up to steady state over many millennia. During the spin-up, basal friction parameters and basin-scale thermal forcing corrections are adjusted to optimize agreement with the observed ice thickness. The model is then run forward for 550 years, from 1950–2500, applying ocean thermal forcing anomalies from six climate models. In all simulations, the ocean forcing triggers long-term retreat of the West Antarctic Ice Sheet, especially in the Filchner–Ronne and Ross sectors. Mass loss accelerates late in the 21st century and then rises steadily for several centuries without leveling off. The resulting ocean-forced sea-level rise at year 2500 varies from about 150 to 1300 mm, depending on the melt scheme and ocean forcing. Further experiments show relatively high sensitivity to the basal friction law, moderate sensitivity to grid resolution and the prescribed collapse of small ice shelves, and low sensitivity to the stress-balance approximation. The Amundsen sector exhibits threshold behavior, with modest retreat under many parameter settings but complete collapse under some combinations of low basal friction and high thermal forcing anomalies. Large uncertainties remain, as a result of parameterized sub-shelf melt rates, simplified treatments of calving and basal friction, and the lack of ice–ocean coupling.

54 ENVIRONMENTAL SCIENCES↗

The Continuous Intercomparison of Radiation Codes (CIRC): Phase I Cases

CIRC aspires to be the successor to ICRCCM (Intercomparison of Radiation Codes in Climate Models). It is envisioned as an evolving and regularly updated reference source for GCM-type radiative transfer (RT) code evaluation with the principle goal to contribute in the improvement of RT parameterizations. CIRC is jointly endorsed by DOE's Atmospheric Radiation Measurement (ARM) program and the GEWEX Radiation Panel (GRP). CIRC's goal is to provide test cases for which GCM RT algorithms should be performing at their best, i.e, well characterized clear-sky and homogeneous, overcast cloudy cases. What distinguishes CIRC from previous intercomparisons is that its pool of cases is based on observed datasets. The bulk of atmospheric and surface input as well as radiative fluxes come from ARM observations as documented in the Broadband Heating Rate Profile (BBHRP) product. BBHRP also provides reference calculations from AER's RRTM RT algorithms that can be used to select the most optimal set of cases and to provide a first-order estimate of our ability to achieve radiative flux closure given the limitations in our knowledge of the atmospheric state.

Oreopoulos, Lazaros↗

The effects of small ice crystals on cirrus infrared radiative properties

Using a model that combines single-scattering properties for spheroidal and hexagonal ice crystals, the thermal infrared radiative properties of cirrus clouds have been investigated. Infrared scattering and absorption properties for randomly oriented spheroids and hexagons are parameterized based on the anomalous diffraction theory and a geometric ray-tracing method, respectively. Using observed ice crystal size distributions, upwelling radiances at the top of cirrus cloudy atmospheres have been computed. Results show that the presence of small ice crystals can produce significant brightness temperature differences between two infrared wavelengths in the 10-micron window. Theoretical results have been compared with observed brightness temperature differences between 8.35 and 11.16 microns and between 11.16 and 12 microns. The observed values were obtained from the high-spectral resolution interferometer sounder. It is shown that the use of the present nonspherical model for ice crystals in radiative transfer calculations leads to a significantly better interpretation of the observed data than does the use of the spherical model.

Takano, Y.↗