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 775 records · Page 43

Improving Convection Trigger Functions in Deep Convective Parameterization Schemes Using Machine Learning

Abstract Deficiencies in convection trigger functions, used in deep convection parameterizations in General Circulation Models (GCMs), have critical impacts on climate simulations. A novel convection trigger function is developed using the machine learning (ML) classification model XGBoost. The large‐scale environmental information associated with convective events is obtained from the long‐term constrained variational analysis forcing data from the Atmospheric Radiation Measurement (ARM) program at its Southern Great Plains (SGP) and Manaus (MAO) sites representing, respectively, continental mid‐latitude and tropical convection. The ML trigger is separately trained and evaluated per site, and jointly trained and evaluated at both sites as a unified trigger. The performance of the ML trigger is compared with four convective trigger functions commonly used in GCMs: dilute convective available potential energy (CAPE), undilute CAPE, dilute dynamic CAPE (dCAPE), and undilute dCAPE. The ML trigger substantially outperforms the four CAPE‐based triggers in terms of the F 1 score metric, widely used to estimate the performance of ML methods. The site‐specific ML trigger functions can achieve, respectively, 91% and 93% F 1 scores at SGP and MAO. The unified trigger also has a 91% F 1 score, with virtually no degradation from the site‐specific training, suggesting the potential of a global ML trigger function. The ML trigger alleviates a GCM deficiency regarding the overprediction of convection occurrence, offering a promising improvement to the simulation of the diurnal cycle of precipitation. Furthermore, to overcome the black box issue of the ML methods, insights derived from the ML model are discussed, which may be leveraged to improve traditional CAPE‐based triggers.

54 ENVIRONMENTAL SCIENCES↗

Constraining Black Carbon Aging in Global Models to Reflect Timescales for Internal Mixing

The radiative effects of black carbon depend critically on its atmospheric lifetime, which is controlled by the rate at which freshly emitted combustion particles become internally mixed with other aerosol components. Global aerosol models strive to represent this process, but the timescale for aerosol mixing is not easily constrained using observations. In this study, we apply a timescale parameterization derived from particle‐resolved simulations to quantify, in a global aerosol model, the timescale for internal mixing. We show that, while highly variable, the average timescale for internal mixing is approximately 3 hr, which is much shorter than the 24‐hr aging timescale traditionally applied in bulk aerosol models. We then use the mixing timescale to constrain the aging criterion in the Modal Aerosol Module. Our analysis reveals that, to best reflect timescales for internal mixing, modal models should assume that particles transition from the hydrophobic (fresh) to the hydrophilic (aged) class once they accumulate a coating thickness equal to four monolayers of sulfuric acid, as opposed to the model's current aging criterion of eight monolayers. We show that, in remote regions like the Arctic and Antarctic, predictions of black carbon loading and its seasonal variation are particularly sensitive to the model representation of aging. By constraining aging in global models to reflect mixing timescales simulated by the particle‐resolved model, we eliminate one of the free parameters governing black carbon's long‐range transport and spatiotemporal distribution.

54 ENVIRONMENTAL SCIENCES↗

Predicting the influence of particle size on the glass transition temperature and viscosity of secondary organic material

Abstract Atmospheric aerosols can assume liquid, amorphous semi-solid or glassy, and crystalline phase states. Particle phase state plays a critical role in understanding and predicting aerosol impacts on human health, visibility, cloud formation, and climate. Melting point depression increases with decreasing particle diameter and is predicted by the Gibbs–Thompson relationship. This work reviews existing data on the melting point depression to constrain a simple parameterization of the process. The parameter $$\xi $$ ξ describes the degree to which particle size lowers the melting point and is found to vary between 300 and 1800 K nm for a wide range of particle compositions. The parameterization is used together with existing frameworks for modeling the temperature and RH dependence of viscosity to predict the influence of particle size on the glass transition temperature and viscosity of secondary organic aerosol formed from the oxidation of $$\alpha $$ α -pinene. Literature data are broadly consistent with the predictions. The model predicts a sharp decrease in viscosity for particles less than 100 nm in diameter. It is computationally efficient and suitable for inclusion in models to evaluate the potential influence of the phase change on atmospheric processes. New experimental data of the size-dependence of particle viscosity for atmospheric aerosol mimics are needed to thoroughly validate the predictions.

36 MATERIALS SCIENCE↗

Scaling Behavior of a Turbulent Kinetic Energy Closure Scheme for the Stably Stratified Atmosphere: A Steady-State Analysis

Abstract We present a turbulent kinetic energy (TKE) closure scheme for the stably stratified atmosphere in which the mixing lengths for momentum and heat are not parameterized in the same manner. The key difference is that, while the mixing length for heat tends toward the stability independent mixing length for momentum in neutrally stratified conditions, it tends toward one based on the Brunt–Väisälä time scale and square root of the TKE in the limit of large stability. This enables a unique steady-state solution for TKE to be obtained, which we demonstrate would otherwise be impossible if the mixing lengths were the same. Despite the model’s relative simplicity, it is shown to perform reasonably well with observational data from the 1999 Cooperative Atmosphere–Surface Exchange Study (CASES-99) using commonly employed model constants. Analyzing the scaling behavior of the nondimensional velocity and potential temperature gradients, or of the stability (correction) functions, reveals that for large stability the present model scales in the same manner as the first-order operational scheme of Viterbo et al. Alternatively, it appears as a blend of two cases of the TKE closure scheme of Baas et al. Critically, because a unique steady-state TKE can be obtained, the present model avoids the nonphysical behavior identified in one of the cases of Baas et al.

Meteorology & Atmospheric Sciences↗

Diabatic Eddy Forcing Increases Persistence and Opposes Propagation of the Southern Annular Mode in MERRA-2

Abstract As a dominant mode of jet variability on subseasonal time scales, the Southern Annular Mode (SAM) provides a window into how the atmosphere can produce internal oscillations on longer-than-synoptic time scales. While SAM’s existence can be explained by dry, purely barotropic theories, the time scale for its persistence and propagation is set by a lagged interaction between barotropic and baroclinic mechanisms, making the exact physical mechanisms challenging to identify and to simulate, even in latest generation models. By partitioning the eddy momentum flux convergence in MERRA-2 using an eddy–mean flow interaction framework, we demonstrate that diabatic processes (condensation and radiative heating) are the main contributors to SAM’s persistence in its stationary regime, as well as the key for preventing propagation in this regime. In SAM’s propagating regime, baroclinic and diabatic feedbacks also dominate the eddy–jet feedback. However, propagation is initiated by barotropic shifts in upper-level wave breaking and then sustained by a baroclinic response, leading to a roughly 60-day oscillation period. This barotropic propagation mechanism has been identified in dry, idealized models, but here we show evidence of this mechanism for the first time in reanalysis. The diabatic feedbacks on SAM are consistent with modulation of the storm-track latitude by SAM, altering the emission temperature and cloud cover over individual waves. Therefore, future attempts to improve the SAM time scale in models should focus on the storm-track location, as well as the roles of the cloud and moisture parameterizations. Significance Statement As they circumnavigate the planet, the tropospheric jet streams slowly drift north and south over about 30 days, longer than the normal limit of weather prediction. Understanding the source of this “memory” could improve our knowledge of how the atmosphere organizes itself and our ability to make long-term forecasts. Current theories have identified several possible internal atmospheric interactions responsible for this memory. Yet most of the theories for understanding the jets’ behavior assume that this behavior is only weakly influenced by atmospheric water vapor. We show that this assumption is not enough to understand jet persistence. Instead, clouds and precipitation are more important contributors in reanalysis data than internal “dry” mechanisms to this memory of the Southern Hemisphere jet.

54 ENVIRONMENTAL SCIENCES↗

Large-eddy simulation of an atmospheric bore and associated gravity wave effects on wind farm performance in the southern Great Plains

Gravity waves are a common occurrence in the atmosphere, with a variety of generation mechanisms. Their impact on wind farms has only recently gained attention, with most studies focused on wind farm-induced gravity waves. In this study, the interaction between a wind farm and gravity waves generated by an atmospheric bore event is assessed using multiscale large-eddy simulations. The atmospheric bore is created by a thunderstorm downdraft from a nocturnal mesoscale convective system (MCS). The associated gravity waves impact the wind resource and power production at a nearby wind farm during the American Wake Experiment (AWAKEN) in the US southern Great Plains. A two-domain nested setup (Δx=300 and 20 m) is used in the Weather Research and Forecasting (WRF) model, forced with data from the High-Resolution Rapid Refresh model, to capture both the formation of the bore and its interaction with individual wind turbines. The MCS is resolved on the large outer domain, where the structure of the bore and the associated gravity waves are found to be especially sensitive to parameterized microphysics processes. On the finer inner domain, gravity wave interactions with individual wind turbines are resolved; wake dynamics are captured using a generalized actuator disk parameterization in WRF. The gravity waves are found to have a strong effect on the atmosphere above the wind farm; however, the effect of the waves is more nuanced closer to the surface where there is additional turbulence, both ambient and wake-generated. Notably, the gravity waves modulate the mesoscale environment by weakening and dissipating the preexisting low-level jet, which reduces hub-height wind speed and hence the simulated power output, which is confirmed by the observed supervisory control and data acquisition (SCADA) power data. Additionally, the gravity waves induce local wind direction variations correlated with fluctuations in pressure, which lead to fluctuations in the simulated power output as various turbines within the farm are subjected to waking from nearby turbines.

17 WIND ENERGY↗

Theoretical Calculation and Validation of the Water Vapor Continuum Absorption

The primary objective of this investigation is the development of an improved parameterization of the water vapor continuum absorption through the refinement and validation of our existing theoretical formalism. The chief advantage of our approach is the self-consistent, first principles, basis of the formalism which allows us to predict the frequency, temperature and pressure dependence of the continuum absorption as well as provide insights into the physical mechanisms responsible for the continuum absorption. Moreover, our approach is such that the calculated continuum absorption can be easily incorporated into satellite retrieval algorithms and climate models. Accurate determination of the water vapor continuum is essential for the next generation of retrieval algorithms which propose to use the combined constraints of multispectral measurements such as those under development for EOS data analysis (e.g., retrieval algorithms based on MODIS and AIRS measurements); current Pathfinder activities which seek to use the combined constraints of infrared and microwave (e.g., HIRS and MSU) measurements to improve temperature and water profile retrievals, and field campaigns which seek to reconcile spectrally-resolved and broad-band measurements such as those obtained as part of FIRE. Current widely used continuum treatments have been shown to produce spectrally dependent errors, with the magnitude of the error dependent on temperature and abundance which produces errors with a seasonal and latitude dependence. Translated into flux, current water vapor continuum parameterizations produce flux errors of order 10 W/sq m, which compared to the 4 W/sq m magnitude of the greenhouse gas forcing and the 1-2 W/sq m estimated aerosol forcing is certainly climatologically significant and unacceptably large. While it is possible to tune the empirical formalisms, the paucity of laboratory measurements, especially at temperatures of interest for atmospheric applications, preclude tuning, the empirical continuum models over the full spectral range of interest for remote sensing and climate applications. Thus, we propose to further develop and refine our existing, far-wing formalism to provide an improved treatment applicable from the near-infrared through the microwave. Based on the results of this investigation, we will provide to the remote sensing/climate modeling community a practical and accurate tabulation of the continuum absorption covering the near-infrared through the microwave region of the spectrum for the range of temperatures and pressures of interest for atmospheric applications.

Ma, Qiancheng↗

Toward Understanding the Simulated Phase Partitioning of Arctic Single-Layer Mixed-Phase Clouds in E3SM

21 Mixed-phase clouds are widespread in the Arctic and have substantial impacts on the 22 surface energy budget, sea ice extent, and regional climate. Significant changes are found in the 23 modeled phase partitioning of Arctic mixed-phase clouds in the newly developed U.S. 24 Department of Energy (DOE) Energy Exascale Earth System Model (E3SM) Atmosphere Model 25 version 1 (EAMv1) compared to its predecessor (CAM5/EAMv0). The mixed-phase clouds 26 simulated by EAMv1 is overly dominated by super-cool liquid which is often underestimated by 27 CAM5. Compared to CAM5, EAMv1 has made several important changes in its cloud physical 28 parameterizations, including using CLUBB to unify shallow convection, PBL turbulence, and 29 cloud macrophysics, as well as updating the MG two-moment cloud microphysics scheme and 30 using the CNT heterogeneous ice nucleation scheme to replace the Meyers scheme for mixed- 31 phase clouds. To understand which of the newly-implemented physical parameterizations is 32 responsible for the change of modeled mixed-phase cloud phase partitioning in EAMv1, 33 sensitivity tests are conducted by utilizing the short-term hindcast approach and results are 34 compared to the observations from the ARM Mixed-Phase Arctic Cloud Experiment (M-PACE) 35 field campaign. 36 It is found that although each scheme contributes to the decrease of cloud ice water 37 content (IWC) and increase of cloud liquid water content (LWC), only the combined effect of 38 CLUBB and CNT can explain the significant underestimation of IWC. By examining the budgets 39 of cloud microphysical processes, CLUBB is found to play a dominant role due to the fact that 40 CLUBB does not diagnoses cloud ice phase water in its assumed joint PDF method and thus may 41 not provide sufficient cloud ice to initialize ice phase microphysical processes such as the growth 42 of ice crystals through vapor deposition. Although a more effective heterogeneous ice nucleation 43 may alleviate the underestimation of IWC, the ice phase process in CLUBB requires urgent 44 attention in the future model development.

Zhang, Meng↗

A Model Intercomparison Study of Aerosol‐Cloud‐Turbulence Interactions in a Cloud Chamber: 1. Model Results

This study presents the first model intercomparison of aerosol‐cloud‐turbulence interactions in a controlled cloudy Rayleigh‐Bénard Convection chamber environment, utilizing the Pi Chamber at Michigan Technological University. We analyzed simulated cloud chamber‐averaged statistics of microphysics and thermodynamics in a warm‐phase, cloudy environment under steady‐state conditions at varying aerosol injection rates. Simulation results from seven distinct models (DNS, LES, and a 1D turbulence model) were compared. Our findings demonstrate that while all models qualitatively capture observed trends in droplet number concentration, mean radius, and droplet size distributions at both high and low aerosol injection rates, significant quantitative differences were observed. Notably, droplet number concentrations varied by over two orders of magnitude between models for the same injection rates, indicating sensitivities to the model treatments in droplet activation and removal and wall fluxes. Furthermore, inconsistencies in vertical relative humidity profiles and in achieving steady‐state liquid water content suggest the need for further investigation into the mechanisms driving these variations. Despite these discrepancies, the models generally reproduced consistent power‐law relationships between the microphysical variables. This model intercomparison underscores the importance of controlled cloud chamber experiments for validating and improving cloud microphysical parameterizations. Recommendations for future modeling studies are also highlighted, including constraining wall conditions and processes, investigating droplet/aerosol removal (including sidewall losses), and conducting simplified experiments to isolate specific processes contributing to model divergence and reduce model uncertainties.

54 ENVIRONMENTAL SCIENCES↗

Improving GCM Representation of Convective Cloud Microphysics by Using ARM Raman Lidar and Cloud Radar Observations

The overall objective of this Department of Energy (DOE) Atmospheric System Research (ASR) funded project is to improve the representation of convective cloud microphysics in global climate models (GCMs) and check it by comparing the model simulations with observations. We have derived the cloud ice water content by synthesizing ARM Raman lidar (RL) and cloud radar observations at the Atmospheric Radiation Measurement (ARM) sites. Noting that the simulated anvil clouds in terms of their macro- and micro-physical properties are sensitive to the parameterization of convective microphysical processes, observed ice water content in anvil clouds provides a useful constraint on these parameterizations. We have improved the convective microphysics parameterization scheme by (1) considering sedimentation for cloud ice crystals that do not fall in the original scheme, (2) applying a new terminal velocity parameterization that depends on the environmental conditions for convective snow, (3) adding a new hydrometeor category, “rimed ice,” to the original four-class (cloud liquid, cloud ice, rain, and snow) scheme, and (4) allowing convective clouds to detrain snow particles into stratiform clouds. We have examined the impact of improved convective cloud microphysics parameterization on the simulated global climate from GCM simulations.

54 ENVIRONMENTAL SCIENCES↗

Atmospheric lapse rate regimes and their parameterization

Lapse rates, moist adiabatic lapse rates and the critical lapse rate for baroclinic adjustment are calculated and compared for the mean annual, January and July states in the Northern Hemisphere. The small seasonal changes in hemispheric mean lapse rates can be attributed to counteracting seasonal changes in the strength of moist convection and baroclinic eddies. A major conclusion is that the vertical temperature structure may be well approximated by a radiative-convective equilibrium model with two critical lapse rates, viz. the moist adiabatic lapse rate and the critical lapse rate for baroclinic adjustment.

Stone, P. H.↗

A new parameterization of 15 micron radiative transfer for a GCM

In order to make the Wu-Kaplan longwave radiative transfer parameterization (Krishnamurthy, 1982) presently used in the 9 layer GLAS GCM more suitable for use at higher horizontal and vertical resolutions, the fixed CO2 transmittance tables and climatological O3 transmittances are replaced with appropriate models. Results of off line tests of simple models of CO2 transmittance as a function of atmospheric temperature profile and surface pressure, based on the technique used by Susskind et al. (1983) in the GLAS physical retrieval scheme are given. The models are evaluated in terms of tansmittance error, flux divergence error, and equilibrium temperature error. The transmittances of CO2 averaged over each of the spectral bands 500-660 cm-1 and 660-800 cm-1 are modeled, following Susskind et al., 1983, as products of effective layer transmittances.

Wobus, R. L.↗

The sensitivity of a general circulation model to Saharan dust heating

Evidence is accumulating that tropospheric aerosols, such as Saharan dust, can significantly influence weather and climate. The Climate Model of the Goddard Laboratory for Atmospheric Sciences is presently used to assess the possible effects of Saharan dust on the weather and climate of North Africa and the tropical Atlantic Ocean. The three-dimensional model solves the conservation equations for the horizontal wind vector, potential temperature, water vapor mixing ratio, and surface pressure. Parameterizations are included for turbulent exchange at the earth surface, cumulus convection, large scale saturation, diurnally varying solar radiation, and terrestrial radiation. Realistic simulations are thereby produced for many aspects of the observed climate.

Randall, D.↗

Comparison of cloud microphysical parameters derived from surface and satellite measurements during FIRE phase 2

Cloud microphysical properties are an important component in climate model parameterizations of water transport, cloud radiative exchange, and latent heat processes. Estimation of effective cloud particle size, liquid or ice water content, and optical depth from satellite-based instrumentation is needed to develop a climatology of cloud microphysical properties and to better understand and model cloud processes in atmospheric circulation. These parameters are estimated from two different surface data sets taken at Coffeyville, Kansas, during the First ISCCP Regional Experiment (FIRE) Phase-2 Intensive Field Observation (IFO) period (November 13 - December 7, 1991). Satellite data can also provide information about optical depth and effective particle size. This paper explores the combination of the FIRE-2 surface and satellite data to determine each of the cloud microphysical properties.

Young, David F.↗

Aerosol-Ice Formation Closure Pilot Study (AEROICESTUDY) Field Campaign Report

The aerosol-ice formation closure pilot study (AEROICESTUDY) was a field campaign conducted at the U.S. Department of Energy Atmospheric Radiation Measurement (ARM) user facility at the Southern Great Plains (SGP) observatory from October 7 to 28, 2019. The purpose of this campaign was to test a field observational approach for conducting an aerosol-ice formation closure study. In other words, the goal was to predict the number of ice nucleating particles (INPs) from the measurements of physical and chemical aerosol properties and compare those to measured INP number concentrations. This project was motivated by the fact that the aerosol community has widely conducted aerosol radiative closure, aerosol-cloud condensation nuclei (CCN), and CCN-droplet closure studies to test the physical models and parameterizations that cloud-resolving and climate models rely on to perform reliable simulations of the Earth system and energy budget. However, very few closure studies related to INPs have been conducted, and to our knowledge, none using robust, size-resolved, ambient aerosol composition measurements as model inputs

AEROICESTUDY↗

The Global Mapping of Electron Precipitation and Ionospheric Conductance From Whistler-Mode Chorus Waves

Auroral precipitation is the second major energy source after solar irradiation that ionizes the Earth’s upper atmosphere. Diffuse electron aurora caused by wave-particle interaction in the inner magnetosphere (L < 8) takes over 60% of total auroral energy flux, strongly contributing to the ionospheric conductance and thus to the ionosphere-thermosphere dynamics. This paper quantifies the impact of chorus waves on the diffuse aurora and the ionospheric conductance during quiet, medium, and strong geomagnetic activities, parameterized by AE <100, 100 < AE < 300, and AE > 300, respectively. Using chorus wave statistics and inner-magnetosphere plasma conditions from Timed History Events and Macroscale Interactions during Substorms (THEMIS) observations, we directly derive the energy spectrum of diffuse electron precipitation under quasi-linear theory. We then calculate the height-integrated conductance from the wave-driven aurora spectrum using the electron impact ionization model of Fang et al. (Geophys. Res. Lett., 2010, 37) and the MSIS atmosphere model. By utilizing Fang’s ionization model, the US Naval Research Laboratory Mass Spectrometer and Incoherent Scattar Radar (NRLMSISE-00) model from 2000s for the neutral atmosphere components, and the University of California, Los Angeles (UCLA) Full Diffusion Code, we improve upon the standard generalization of Maxwellian diffuse electron precipitation patterns and their resulting ionosphere conductance. Our study of global auroral precipitation and ionospheric conductance from chorus wave statistics is the first statistical model of its kind. We show that the total electron flux and conductance pattern from our results agree with those of Ovation Prime model over the pre-midnight to post-dawn sector as geomagnetic activity increases. Our study examines the relative contributions of upper band chorus (UBC) and lower band chorus wave (LBC) driven conductance in the ionosphere. We found LBC waves drove diffuse electron precipitation significantly more than UBC waves, however it is possible that THEMIS data may have underestimated the upper chorus band wave observations for magnetic latitudes below 65 degrees. ° .

diffuse auroral precipitation↗

Improved Characterization of PSC Processes Derived from a Third-Generation CALIOP and MLS Detection and Composition Classification Algorithm

The new 3-year CloudSat and CALIPSO Science Team project described in this poster will use a unique combination of data from the Cloud-Aerosol LIdar with Orthogonal Polarization (CALIOP) instrument on CALIPSO and the Microwave Limb Sounder (MLS) on Aura, in conjunction with supporting meteorological information and detailed modeling studies, to advance our understanding of polar stratospheric cloud (PSC) processes and their role in ozone depletion. We will develop a third-generation (Gen3) PSC detection and composition algorithm that incorporates a new, more robust two-dimensional, multi-channel CALIOP feature detection scheme (2D-McDA). We will also devise and implement an improved two-dimensional PSC composition classification scheme that utilizes multiple parameters (e.g., CALIOP 532-nm parallel and perpendicular scattering ratios, CALIOP 1064-nm total scattering ratio, MLS HNO3 and H2O, ambient temperature, and temperature histories) in a Bayesian approach to determine the most likely PSC composition and help constrain solid PSC particle number density and size/shape. The combined CALIOP/MLS analyses will allow us to study in detail the full life cycle of PSCs and their resulting impact on gas-phase HNO3 and H2O, which should lead to improved parameterizations of PSC microphysics in global CCMs where detailed particle information is not available. The Gen3 CALIOP PSC algorithm will be a natural stepping-stone toward the analysis of data collected during future spaceborne lidar missions, such as NASA’s Atmosphere Observation System (AtmOS) mission currently scheduled for launch late in this decade. We will also investigate possible trends in PSC occurrence and composition over the entire CALIOP data record and through further comparisons with the Stratospheric Aerosol Measurement (SAM) II solar occultation PSC record from 1979-1989. Finally, we will validate the mountain-wave parameterization and PSC schemes used in the UM-UKCA (Unified Model coupled to the United Kingdom Chemistry and Aerosol module) chemistry-climate model through detailed comparisons with earlier CALIOP PSC data products and those developed under this proposal.

CALIPSO↗

Investigate wave-mean flow interaction and transport in the extratropical winter stratosphere

The grant supported studies using several models along with observations in order to investigate some questions of wave-mean flow interaction and transport in the extratropical winter stratosphere. A quasi-geostrophic wave model was used to investigate the possibility that resonant growth of planetary wave 2 may have played a role in the sudden stratospheric warming of February 1979. The results of the time-dependent integration support the interpretation of resonance during February, 1979. Because of the possibility that the model treatment of critical line interactions exerted a controlling influence on the atmospheric dynamics, a more accurate model was needed for wave-mean flow interaction studies. A new model was adapted from the 3-dimensional primitive equation model developed by K. Rose and G. Brasseur. In its present form the model is global, rather than hemispheric; it contains an infrared cooling algorithm and a parameterized solar heating; it has parameterized gravity wave drag; and the chemistry has been entirely revised.

Smith, Anne K.↗