Search NASA⌕ Search

SEARCH · Search NASA

Results for “general circulation model”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 973 records · Page 54
All sourcesNASA NTRSDOE OSTI
NASA NTRS1999

Global Aerosol Radiative Forcing Derived from Sea WiFS-Inferred Aerosol Optical Properties

Aerosol optical properties inferred from the Sea-viewing Wide Field-of-view Sensor (SeaWiFS) radiance measurements are used to compute the aerosol shortwave radiative forcing using a radiative transfer model. The aerosol optical thickness at the wavelength of 865-nm is taken from the SeaWIFS archive. It is found that the nominal optical thickness over oceans ranges from 0.1 to 0.2. Using a maritime aerosol model and the radiances measured at the various SeaWiFS channels, the Angstrom exponent is determined to be 0.2174, the single-scattering albedo to be 0.995, and the asymmetry factor to be 0.786. The radiative transfer model has eight bands in the visible and ultraviolet spectral regions and three bands in the near infrared. It includes the absorption due to aerosols, water vapor, carbon dioxide, and oxygen, and the scattering due to aerosols and gases (Rayleigh scattering). The radiative forcing is computed over global oceans for four months (January, April, July, and October, 1998) to represent four seasons. It is found that the aerosol radiative forcing is large and changes significantly with seasons near the continents with large-scale forest fires and desert dust. Averaged over oceans and the four months, the aerosol radiative forcing is approximately 7 W/sq m at the top of the atmosphere. This large radiative forcing is expected to have a significant cooling effect on the Earth's climate as implied from simulations of a number of general circulation models.

Chou, Ming-Dah↗
NASA NTRS1986

Initial results of the Global Thermospheric Mapping Study (GTMS)

The Global Thermospheric Mapping Study (GTMS) is a multi-technique experimental study of the thermosphere designed to map simultaneously its spatial and temporal morphology with a thoroughness and diversity of measurement techniques heretofore unachieved. The GTMS is designed around the Incoherent Scatter Radar Chain in the western hemisphere. The European incoherent scatter radars and the worldwide communities of Fabry-Perot interferometers, meteor wind radars, partial reflection drifts radars, MST radars, and satellite probes are included to extend the spatial coverage and types of measurements available. Theoretical and modeling support in the areas of thermospheric and ionospheric structure, tides, and electric fields are included to aid in program planning and data interpretation. Solar activity was low on the three observation days (F10.7 = 97, 98, 96) and magnetic conditions were unsettled to active (A = 10, 12, 20). All six incoherent scatter radar facilities collected data. Each collected F region data day and night while Saint Santin and Millstone Hill additionally collected E region data during daylight hours. Initial results from Sondrestrom and Millstone Hill are presented. Good quality Fabry Perot data were collected at Fritz Peak and San Jose dos Campos. Weather conditions produced poor results at Arequipa and Arecibo. Initial results from Fritz Peak are presented. Mesosphere/lower-thermosphere observations were conducted under the ATMAP organization. The magnetometer chains also were operational during this campaign. Initial thermospheric general circulation model predictions were made for assumed solar-geophysical conditions, and selected results are presented.

Oliver, W. L.↗
NASA NTRS2004

The Contribution of Dynamic Interannual Variability to Ozone Trends

At middle latitudes the total column ozone and the lower stratospheric ozone mixing ratio exhibit natural variability. Models and statistical analyses of observations such as SAGE ozone profiles and TOMS column measurements show that seasonal cycle, solar cycle, and interannual dynamical variability and dynamical phenomena such as the quasi- biennial oscillation all contribute to ozone variability. These must be accounted for when deriving ozone trends. Systematic or random changes in the atmospheric circulation may also contribute to ozone trends. It is presently unclear how much of the ozone trend derived from observations is due to changes in the chemical composition of the stratosphere and how much is due to changes in the atmospheric circulation. We are attempting to resolve this issue by comparing a twenty-five year simulation of ozone with fixed source gas boundary conditions with an identical simulation with time dependent source gas boundary conditions. Both simulations are driven with output from a general circulation model that produces realistic interannual variability in dynamical forcing. The model trend in ozone due to changes in composition is determined from the difference in these simulations. We compare these trends with trends determined from observations and model output using the same analysis techniques. Initial results emphasize the complications to attribution of observed ozone trends to dynamical and photochemical effects that are due to interrelationships between trends in transport, temperature, and photochemical effects. It may not be possible to describe the ozone trend as a superposition of dynamical and photochemical contributions.

Douglass, Anne↗
NASA NTRS2013

Wavelength Dependence of Solar Irradiance Enhancement During X-Class Flares and Its Influence on the Upper Atmosphere

The wavelength dependence of solar irradiance enhancement during flare events is one of the important factors in determining how the Thermosphere-Ionosphere (T-I) system responds to flares. To investigate the wavelength dependence of flare enhancement, the Flare Irradiance Spectral Model (FISM) was run for 61 X-class flares. The absolute and the percentage increases of solar irradiance at flare peaks, compared to pre-flare conditions, have clear wavelength dependences. The 0-14 nm irradiance increases much more (approx. 680% on average) than that in the 14-25 nm waveband (approx. 65% on average), except at 24 nm (approx. 220%). The average percentage increases for the 25-105 nm and 122-190 nm wavebands are approx. 120% and approx. 35%, respectively. The influence of 6 different wavebands (0-14 nm, 14-25 nm, 25-105 nm, 105- 120 nm, 121.56 nm, and 122-175 nm) on the thermosphere was examined for the October 28th, 2003 flare (X17-class) event by coupling FISM with the National Center for Atmospheric Research (NCAR) Thermosphere-Ionosphere-Electrodynamics General Circulation Model (TIE-GCM) under geomagnetically quiet conditions (Kp=1). While the enhancement in the 0-14 nm waveband caused the largest enhancement of the globally integrated solar heating, the impact of solar irradiance enhancement on the thermosphere at 400 km is largest for the 25-105 nm waveband (EUV), which accounts for about 33 K of the total 45 K temperature enhancement, and approx. 7.4% of the total approx. 11.5% neutral density enhancement. The effect of 122-175 nm flare radiation on the thermosphere is rather small. The study also illustrates that the high-altitude thermospheric response to the flare radiation at 0-175 nm is almost a linear combination of the responses to the individual wavebands. The upper thermospheric temperature and density enhancements peaked 3-5 h after the maximum flare radiation.

TIE-GCM↗
NASA NTRS2013

Wavelength Dependence of Solar Irradiance Enhancement During X-class Flares and Its Influence on the Upper Atmosphere

The wavelength dependence of solar irradiance enhancement during flare events is one of the important factors in determining how the Thermosphere-Ionosphere (TI) system responds to flares. To investigate the wavelength dependence of flare enhancement, the Flare Irradiance Spectral Model (FISM) was run for 61X-class flares. The absolute and the percentage increases of solar irradiance at flare peaks, compared to pre-flare conditions, have clear wavelength dependences. The 0-4 nm irradiance increases much more ((is) approximately 680 on average) than that in the 14-25 nm waveband ((is) approximately 65 on average), except at 24 nm ( (is) approximately 220). The average percentage increases for the 25-105 nm and 122-190 nm wave bands are approximately 120 and approximately 35, respectively. The influence of 6 different wavebands (0-14 nm, 14-25 nm, 25-105 nm, 105-120 nm, 121.56 nm,and122-175 nm) on the thermosphere was examined for the October 28th, 2003 flare (X17-class) event by coupling FISM with the National Center for Atmospheric Research (NCAR) Thermosphere-Ionosphere-Electrodynamics General Circulation Model(TIE-GCM) under geomagnetically quiet conditions (Kp=1). While the enhancement in the0-14nm waveband caused the largest enhancement of the globally integrated solar heating, the impact of solar irradiance enhancement on the thermosphere at 400 km is largest for the 25-105 nm waveband (EUV), which accounts for about 33 K of the total 45 K temperature enhancement, and approximately 7.4% of the total approximately 11.5% neutral density enhancement. The effect of 122-175 nm flare radiation on the thermosphere is rather small. The study also illustrates that the high-altitude thermospheric response to the flare radiation at 0-175 nm is almost a linear combination of the responses to the individual wavebands. The upper thermospheric temperature and density enhancements peaked 3-5 h after the maximum flare radiation.

FISM↗
NASA NTRSDate not supplied

Secondary Volcanically-Induced Lunar Atmosphere and Lunar Volatiles: 3-D Modeling and Analysis

The fact that the Moon could have a transient secondary atmosphere due to volcanic outgassing has been known for some time, though typically such an atmosphere was believed to be extremely thin (~10-8 bar) [1]. But recent research by Needham and Kring (NK) [2] suggests that during the peak of volcanic activity ~3.5 Ga such a volcanically-outgassed atmosphere could reach ~10-2 bar of surface pressure. In similar research Wilson et al. [3] proposed a more conservative estimate, arguing that the thickness of such an atmosphere would depend on the intervals between major eruptions and may not exceed microbar densities. In either case a collisional atmosphere could be present, which would control transport of outgassed volatiles (such as H2O) and their deposition in polar regions, where they could be preserved until modern day frozen in permanently shadowed regions (PSR) or buried beneath the regolith. Here we study such a hypothetical atmosphere to investigate its stability, meteorological properties and the effect on transport of volatiles. We use the ROCKE-3D planetary 3-D General Circulation Model (GCM)[4]. The insolation and orbital parameters were set to conditions 3.5 Ga. The atmospheric composition, based on the list of outgassed species presented by NK in combination with our estimates for atmospheric escape, condensation and the results from our 1-D chemistry model, was chosen to be either CO-dominated or CO2-dominated (depending on atmospheric temperature). In this study we restricted ourselves to relatively "thick" lunar atmospheres of 1-10 mb, though we believe that our results will scale to thinner atmospheres as well. We present the results for ground and atmospheric temperature for modeled atmospheres over a wide parameter space. In particular we consider different atmospheric compositions (CO or CO2 dominated), a set of atmospheric pressures from 1 mb to 10 mb and a set of obliquities from 0o to 40o. We also present an experiment of a single major eruption [5] and show that in just 3 years ~80% of the outgassed water is deposited in polar regions. This demonstrates the efficiency of such an atmosphere in delivering volatiles. We argue that a secondary lunar atmosphere could play a significant role in forming volatile deposits currently observed in the polar regions of the Moon.

Lunar atmospheres↗
DOE OSTI2022

How May the Choice of Downscaling Techniques and Meteorological Reference Observations Affect Future Hydroclimate Projections?

Abstract We present an intercomparison of a suite of high‐resolution downscaled climate projections based on a six‐member General Circulation Model (GCM) ensemble from Coupled Models Intercomparison Project (CMIP6). The CMIP6 GCMs have been downscaled using dynamical and statistical downscaling techniques based on two meteorological reference observations over the conterminous United States. We use the regional climate model, RegCM4, for dynamical downscaling, double bias correction constructed analogs method for statistical downscaling, and Daymet and Livneh datasets as the reference observations for statistical training and bias‐correction. We evaluate the performances of downscaled data in both historical and future periods under the SSP585 scenario. While dynamical downscaling improves the simulation of some performance evaluation indices, it adds an extra bias in others, highlighting the need for statistical correction before its use in impact assessments. Downscaled datasets after bias‐correction compare exceptionally well with observations. However, the choice of downscaling techniques and the underlying reference observations influence the hydroclimate characteristics of downscaled data. For instance, the statistical downscaling generally preserves the GCMs climate change signal but overestimates the frequency of hot extremes. Similarly, simulated future changes are sensitive to the choice of reference observations, particularly for precipitation extremes that exhibit a higher projected increase in the ensembles trained and/or corrected by Daymet than Livneh. Overall, these results demonstrate that multiple factors, including downscaling techniques and reference observations, can substantially influence the outcome of downscaled climate projections and stress the need for a comprehensive understanding of such method‐based uncertainties.

54 ENVIRONMENTAL SCIENCES↗
NASA NTRS1994

Quasi-biennial oscillations of ozone and diabatic circulation in the equatorial stratosphere

The quasi-biennial oscillation (QBO) in ozone in the equatorial stratosphere is obtained by analyzing the Stratospheric Aerosol and Gas Experiment (SAGE) data from 1984 to 1989. The phase of the ozone QBO in the lower stratosphere is found to precede the zonal wind QBO by several months as opposed to the theoretically expected in-phase relationship between the two. A mechanistic model is developed to explore possible reasons for this disagreement. The model is capable of simulating the actual time evolution of the ozone QBO by introducing the observed zonal wind profile as input. The modeled results confirm the conventional view that the ozone QBO is generated by the vertical ozone advection that is driven to maintain the temperature structure against radiative damping. However, a series of experiments emphasizes the importance of the feedback of the ozone QBO to the diabatic heating through the absorption of solar radiation. Due to this effect, the phase of the ozone QBO shifts up to a quarter cycle ahead and approaches that of the temperature QBO. Because of this inphase relationship, the feedback of the ozone QBO to the diabatic heating acts to compensate for the radiative damping of the temperature structure, thus reducing the magnitude of the induced diabatic circulation. Because the reduction of the magnitude of the vertical motion facilitates downward transport of easterly momentum by the mean flow, this feedback process can help to resolve the insufficiency of the easterly momentum in driving the dynamical QBO in general circulation models (GCMs). It should be emphasized that more sophisticated models that allow for full interaction between the chemical species and radiative and dynamical processes should be developed to improve our understanding of both dynamical and ozone QBOs.

Hasebe, Fumio↗
NASA NTRS1999

Comparison of Radiative Forcing Calculations Due to Mineral Dust from a Transport Model, Satellite Measurements and an Assimilation System

This study uses information on mineral aerosol from a transport model to calculate global radiative forcing values. The transport model is driven by assimilated meteorology and outputs three-dimensional dust spatial information for various size ranges. The dust fields are input to an off-line radiative transfer calculation to obtain the direct radiative forcing due to the dust fields. During June, July and August of 1988 presence of dust 1) reduces the global net incoming radiation at the top of atmosphere (TOA) by 0.3 to 0.7 W/sq m and 2) reduces net incoming radiation at the earth's surface by 1.3 to 2.0 W/sq m. Over Africa our estimates of the reduction of radiation at the top of atmosphere compare well with TOA reductions derived from ERBE and TOMS satellite data. However, our heating rates are not consistent with analysis temperature increments produced by the assimilation system over regions of high aerosol loading. These increments are based on differences between temperature observations and temperatures from the assimilation general circulation model. One explanation is that the lower tropospheric temperatures retrieved by TOVS are being contaminated by mineral aerosol.

Weaver, Clark J.↗
NASA NTRS1995

Validation of clear-sky fluxes for tropical oceans from the Earth Radiation Budget Experiment

The existence and magnitude of a systematic bias in the clear-sky longwave fluxes from the Earth Radiation Budget Experiment (ERBE) is investigated. The bias is apparently introduced because the ERBE method for scene identification does not account for large zonal gradients in longwave absorption by water vapor. The ERBE fluxes are compared to fluxes calculated with a radiative transfer model from ship radiosonde measurements. The comparison is based upon an analysis of 5 yr of coincident satellite and radiosonde observations for equatorial ocean regions. The differences between the ERBE and model fluxes are examined as functions of sea surface temperature (SST) and relative humidity. The authors use height-mean relative humidity bar-RH as an index of atmospheric moisture. The average offset between and model ERBE fluxes ranges between +2 and +6 W/sq m for SSTs above 295 K, and the gradients with respect to SST are nearly identical. However, the difference between the model and ERBE depends significantly on the tropospheric relative humidity. ERBE fluxes exceed model fluxes for bar-RH above 70%, and the maximum offset of +9 to +12 W/sq m is consistent with previous estimates. There are also indications that the clear-sky fluxes for bar-RH below 25% may be underestimated by about 10-15 W/sq m. Since extreme values of height-mean humidity are relatively infrequent, the net bias introduced in the ERBE monthly mean clear-sky fluxes is generally less than the systematic error in estimates introduced in the ERBE monthly mean clear-sky fluxes is generally less than the systematic error in estimates of the instantaneous fluxes. These findings support earlier work on the coupling between SST and the atmospheric greenhouse effect, in particular the existence of a super greenhouse effect for oceans warmer than 300 K. Recent reports of much larger systematic differences are not supported by this analysis. The results indicate that comparison of general circulation model (GCM) and ERBE clear-sky longwave fluxes will depend explicitly on atmospheric humidity.

Collins, W. D.↗
NASA NTRS1980

The effect of clouds on the earth's solar and infrared radiation budgets

The effect of global cloudiness on the solar and infrared components of the earth's radiation balance is studied in general circulation model experiments. A wintertime simulation is conducted in which the cloud radiative transfer calculations use realistic cloud optical properties and are fully interactive with model-generated cloudiness. This simulation is compared to others in which the clouds are alternatively non-interactive with respect to the solar or thermal radiation calculations. Other cloud processes (formation, latent heat release, precipitation, vertical mixing) were accurately simulated in these experiments. It is concluded that on a global basis clouds increase the global radiation balance by 40 W/sq m by absorbing longwave radiation, but decrease it by 56 W/sq m by reflecting solar radiation to space. The net cloud effect is therefore a reduction of the radiation balance by 16 W/sq m, and is dominated by the cloud albedo effect. Changes in cloud frequency and distribution and in atmospheric and land temperatures are also reported for the control and for the non-interactive simulations. In general, removal of the clouds' infrared absorption cools the atmosphere and causes additional cloudiness to occur, while removal of the clouds' solar radiative properties warms the atmosphere and causes fewer clouds to form. It is suggested that layered clouds and convective clouds over water enter the climate system as positive feedback components, while convective clouds over land enter as negative components.

Herman, G. F.↗
NASA NTRSDate not supplied

Observations of the Age of Air from the Northern Hemisphere Midlatitude Surface

Understanding and modeling the interhemispheric transport (IHT) of trace gases and aerosols is important for interpreting past (and projecting future) changes in atmospheric composition. One IHT measure -- the mean age since air was last at the Northern Hemisphere (NH) midlatitude surface – represents a fundamental property of tropospheric transport that can be estimated from measurements of sulfur hexafluoride (SF6). Here we use SF6 measurements obtained from the recent NASA Atmospheric Tomography Mission (ATom) to construct observational estimates of the mean age throughout the troposphere. We compare the ATom based ages, which span all four seasons during the period 2016-2018 over both the Pacific and Atlantic Oceans, with previous in situ-based estimates obtained from the HIAPER Pole‐to‐Pole Observations (HIPPO) project, which are available only for the Pacific. In order to interpret the observed variations in mean age and their relationship with the large-scale flow we use simulations of the Goddard Earth Observing System (GEOS) general circulation model, constrained with meteorological fields derived from the Modern-Era Retrospective analysis for Research and Applications, Version 2 (MERRA-2), performed over the entire MERRA-2 period (including ATom). We confirm results from previous studies showing that the model ages are biased old compared to the observational estimates. However, we also show that this model bias is reduced in more recent years, which we interpret in terms of changes in the latitudinal distribution of SF6 emissions over the 2000s.

interhemispheric transport (IHT)↗
NASA NTRS2020

The HadGEM3-GA7.1 radiative kernel: the importance of a well-resolved stratosphere

We present top-of-atmosphere and surface radiative kernels based on the atmospheric component (GA7.1) of the HadGEM3 general circulation model developed by the UK Met Office. We show that the utility of radiative kernels for forcing adjustments in idealised CO2 perturbation experiments is greatest where there is sufficiently high resolution in the stratosphere in both the target climate model and the radiative kernel. This is because stratospheric cooling to a CO2 perturbation continues to increase with height, and low-resolution or low-top kernels or climate model output are unable to fully resolve the full stratospheric temperature adjustment. In the sixth phase of the Coupled Model Intercomparison Project (CMIP6), standard atmospheric model data are available up to 1 hPa on 19 pressure levels, which is a substantial advantage compared to CMIP5. We show in the IPSL-CM6A-LR model where a full set of climate diagnostics are available that the HadGEM3-GA7.1 kernel exhibits linear behaviour and the residual error term is small, as well as from a survey of kernels available in the literature that in general low-top radiative kernels underestimate the stratospheric temperature response. The HadGEM3-GA7.1 radiative kernels are available at https://doi.org/10.5281/zenodo.3594673 (Smith, 2019).

radiative kernels↗
NASA NTRS2004

Stratospheric Impact of Varying Sea Surface Temperatures

The Finite-Volume General Circulation Model (FVGCM) has been run in 50 year simulations with the: 1) 1949-1999 Hadley Centre sea surface temperatures (SST), and 2) a fixed annual cycle of SSTs. In this presentation we first show that the 1949-1999 FVGCM simulation produces a very credible stratosphere in comparison to an NCEP/NCAR reanalysis climatology. In particular, the northern hemisphere has numerous major and minor stratospheric warming, while the southern hemisphere has only a few over the 50-year simulation. During the northern hemisphere winter, temperatures are both warmer in the lower stratosphere and the polar vortex is weaker than is found in the mid-winter southern hemisphere. Mean temperature differences in the lower stratosphere are shown to be small (less than 2 K), and planetary wave forcing is found to be very consistent with the climatology. We then will show the differences between our varying SST simulation and the fixed SST simulation in both the dynamics and in two parameterized trace gases (ozone and methane). In general, differences are found to be small, with subtle changes in planetary wave forcing that lead to reduced temperatures in the SH and increased temperatures in the NH.

Newman, Paul A.↗
NASA NTRS1994

A stochastic precipitation disaggregation scheme for GCM applications

In the surface hydrologic parameterization of general circulation models (GCMs), it is commonly assumed that the precipitation processes are homogeneous over a GCM grid square and that the precipitation intensity is uniformly distributed. Based on evidence that the spatial distribution of precipitation within a GCM grid square is crucial for the land surface hydrology parameterization, a few researchers have explored the impacts of assuming that the precipitation is exponentially distributed. This paper explores the suitability of the aforementioned assumptions. First, a statistical analysis is conducted of historical precipitation data for three GCM grids in different regions of the United States. The analysis suggests that neither the uniform nor the exponential distribution assumption may be suitable at the GCM grid scale and, that instead, the spatial variability in precipitation is characterized by statistical patterns that are inhomogeneous. These patterns vary from grid to grid and are induced by the interaction between atmospheric conditions and various land surface characteristics, such as topographical features, surface properties, etc. Within the same grid square, however, the statistical patterns are generally constant from year to year. Based on this analysis, a computationally viable (i.e., usable with GCMs) stochastic precipitation disaggregation scheme that utilizes these stable statistical patterns is proposed. The method was used to generate spatially distributed hourly rainfall for a summer season in the southwestern region of the continental United States. Analysis of the results shows that the methodology preserves the seasonal characteristics of spatial variability in precipitation that is observed in the long-term historical data.

Gao, Xiaogang↗
NASA NTRS2021

A Phenomenon-Based Decomposition of Model-Based Estimates of Boreal Winter ENSO Variability

Climate models are now routinely being used to simulate and predict climate variability on time scales ranging from sub-seasonal to seasonal and longer. As such, there are now long histories of such simulations and predictions spanning multiple decades and multiple ensemble members, both of which are crucial for separating climate signal from climate noise. A key focus of such runs has been the El Niño-Southern Oscillation (ENSO), spurred by recent improvements in our ability to predict such events, though questions remain as to how well climate models do beyond simply always predicting the “canonical” atmospheric response to an ENSO event—something simple statistical models already do reasonably well. This is a critical issue that needs addressing, given the importance of event-to-event differences for predicting regional impacts of ENSO teleconnections, and the need to justify the expense of running sophisticated climate models. Unfortunately, current diagnostic tools are not well suited for quantifying the different sources of variability associated with specific phenomena such as ENSO. More generally, while much effort has focused on addressing model bias, less has been done to address errors in second moment statistics—an issue whose importance is gaining increased attention particularly as we build climate prediction systems capable of taking advantage of forecasts of opportunity—a capability that requires reliable estimates of forecast uncertainty. In this report, we outline a phenomenon-based statistical decomposition of climate variance(in essence a detailed variance budget)that is specifically tailored to address the above questions by separating the variability (both the signal and noise) into that tied to the long-term average impact of a particular phenomenon(the composite mean) and the event-to-event(E2E) variability about the composite mean. In addition, we provide related decompositions of the correlations that allow us to quantify how much of the agreement with observations (the skill) comes from the composite mean as opposed to from the E2Evariability. As an example, we present the results of such a decomposition for ENSO based on simulations with the GEOS atmospheric general circulation model (AGCM), with a focus on the monthly mean impacts over North America during boreal winter(December –March). Here we take advantage of existing GEOS AGCM simulations that were produced as companion simulations to MERRA-2for the period 1980-2016. Comparisons are made throughout with MERRA-2.

ENSO↗
NASA NTRS2017

High-order Accurate Finite-Volume Formulations for the Pressure Gradient Force in Layered Ocean Models

Discretisation of the horizontal pressure gradient force in layered ocean models is a challenging task, with non-trivial interactions between the thermodynamics of the fluid and the geometry of the layers often leading to numerical difficulties. We present two new finite-volume schemes for the pressure gradient operator designed to address these issues. In each case, the horizontal acceleration is computed as an integration of the contact pressure force that acts along the perimeter of an associated momentum control-volume. A pair of new schemes are developed by exploring different control-volume geometries. Non-linearities in the underlying equation-of-state definitions and thermodynamic profiles are treated using a high-order accurate numerical integration framework, designed to preserve hydrostatic balance in a non-linear manner. Numerical experiments show that the new methods achieve high levels of consistency, maintaining hydrostatic and thermobaric equilibrium in the presence of strongly-sloping layer geometries, non-linear equations-of-state and non-uniform vertical stratification profiles. These results suggest that the new pressure gradient formulations may be appropriate for general circulation models that employ hybrid vertical coordinates andor terrain-following representations.

pressure gradient force↗
DOE OSTI2023

Improving the representation of shallow cumulus convection with the simplified-higher-order-closure–mass-flux (SHOC+MF v1.0) approach

Abstract. Parameterized boundary layer turbulence and moist convection remain some of the largest sources of uncertainty in general circulation models. High-resolution climate modeling aims to reduce that uncertainty by explicitly attempting to resolve deep moist convective motions. An example of such a model is the Simple Cloud-Resolving E3SM Atmosphere Model (SCREAM) with a target global resolution of 3.25 km, allowing for a more accurate representation of complex mesoscale deep convective dynamics. Yet, small-scale planetary boundary layer turbulence and shallow convection still need to be parameterized, which in SCREAM is accomplished through the turbulent-kinetic-energy-based (TKE-based) simplified higher-order closure (SHOC) – a simplified version of the assumed-double-Gaussian-PDF (probability density function) higher-order-closure method. In this paper, we implement a stochastic-multiplume mass-flux (MF) parameterization of dry and shallow convection in SCREAM to go beyond the limitations of double-Gaussian-PDF closures and couple it to SHOC (SHOC+MF). The new parameterization implemented in a single-column model type version of SCREAM produces results for two shallow cumulus convection cases (marine and continental shallow convection) that agree well with the reference data from large-eddy simulations, thus improving the general representation of the thermodynamic quantities and their turbulent fluxes as well as cloud macrophysics in the model. Furthermore, SHOC+MF parameterization shows weak sensitivity to the vertical grid resolution and model time step.

54 ENVIRONMENTAL SCIENCES↗