Search NASA⌕ Search

SEARCH · Search NASA

Results for “Vertical coordinates”

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 19 records
All sourcesNASA NTRSDOE OSTI
NASA NTRS2002

Application of the generalized vertical coordinate ocean model for better representing satellite data

It is found that two adaptive parametric functions can be introduced into the basic ocean equations for utilizing the optimal or hybrid features of commonly used z-level, terrain- following, isopycnal, and pressure coordinates in numerical ocean models. The two parametric functions are formulated by combining three techniques: the arbitrary vertical coordinate system of Kasahara (1 974), the Jacobian pressure gradient formulation of Song (1 998), and a newly developed metric factor that permits both compressible (non-Boussinesq) and incompressible (Boussinesq) approximations. Based on the new formulation, an adaptive modeling strategy is proposed and a staggered finite volume method is designed to ensure conservation of important physical properties and numerical accuracy. Implementation of the combined techniques to SCRUM (Song and Haidvogel1994) shows that the adaptive modeling strategy can be applied to any existing ocean model without incurring computational expense or altering the original numerical schemes. Such a generalized coordinate model is expected to benefit diverse ocean modelers for easily choosing optimal vertical structures and sharing modeling resources based on a common model platform. Several representing oceanographic problems with different scales and characteristics, such as coastal canyons, basin-scale circulation, and global ocean circulation, are used to demonstrate the model's capability for multiple applications. New results show that the model is capable of simultaneously resolving both Boussinesq and non-Boussinesq, and both small- and large-scale processes well. This talk will focus on its applications of multiple satellite sensing data in eddy-resolving simulations of Asian Marginal Sea and Kurosio. Attention will be given to how Topex/Poseidon SSH, TRMM SST; and GRACE ocean bottom pressure can be correctly represented in a non- Boussinesq model.

numerical ocean model generalized coordinate syste↗
NASA NTRS1991

Constant covariance in local vertical coordinates for near-circular orbits

A method is presented for devising a covariance matrix that either remains constant or grows in keeping with the presence of a period error in a rotating local-vertical coordinate system. The solution presented may prove useful in the initialization of simulation covariance matrices for near-circular-orbit problems. Use is made of the Clohessy-Wiltshire equations and the travelling-ellipse formulation.

Shepperd, Stanley W.↗
NASA NTRS2001

A Finite-Volume "Shaving" Method for Interfacing NASA/DAO''s Physical Space Statistical Analysis System to the Finite-Volume GCM with a Lagrangian Control-Volume Vertical Coordinate

Toward the development of a finite-volume Data Assimilation System (fvDAS), a consistent finite-volume methodology is developed for interfacing the NASA/DAO's Physical Space Statistical Analysis System (PSAS) to the joint NASA/NCAR finite volume CCM3 (fvCCM3). To take advantage of the Lagrangian control-volume vertical coordinate of the fvCCM3, a novel "shaving" method is applied to the lowest few model layers to reflect the surface pressure changes as implied by the final analysis. Analysis increments (from PSAS) to the upper air variables are then consistently put onto the Lagrangian layers as adjustments to the volume-mean quantities during the analysis cycle. This approach is demonstrated to be superior to the conventional method of using independently computed "tendency terms" for surface pressure and upper air prognostic variables.

Lin, Shian-Jiann↗
DOE OSTI2025

Water Mass Transformation Budgets in Finite‐Volume Generalized Vertical Coordinate Ocean Models

Water Mass Transformation (WMT) theory provides conceptual tools that in principle enable innovative analyses of numerical ocean models; in practice, however, these methods can be challenging to implement and interpret, and therefore remain under-utilized. Our aim is to demonstrate the feasibility of diagnosing all terms in the water mass budget and to exemplify their usefulness for scientific inquiry and model development by quantitatively relating water mass changes, overturning circulations, boundary fluxes, and interior mixing. We begin with a pedagogical derivation of key results of classical WMT theory. We then describe best practices for diagnosing each of the water mass budget terms from the output of Finite-Volume Generalized Vertical Coordinate (FV-GVC) ocean models, including the identification of a non-negligible remainder term as the spurious numerical mixing due to advection scheme discretization errors. We illustrate key aspects of the methodology through the analysis of a polygonal region of the Greater Baltic Sea in a regional demonstration simulation using the Modular Ocean Model v6 (MOM6). We verify the convergence of our WMT diagnostics by brute-force, comparing time-averaged (“offline”) diagnostics on various vertical grids to timestep-averaged (“online”) diagnostics on the native model grid. Finally, we briefly describe a stack of xarray-enabled Python packages for evaluating WMT budgets in FV-GVC models (culminating in the new xwmb package), which is intended to be model-agnostic and available for community use and development.

54 ENVIRONMENTAL SCIENCES↗
NASA NTRS1992

Adaptation of a matrix parameterization of the middle atmospheric radiative cooling for an arbitrary vertical coordinate grid

Matrices corresponding to a revised vertical grid are shown to be derived from the original grid by integrating the radiative cooling rate over model layers. The method is applied to temperature profiles that represent average or extreme atmospheric conditions and a range of vertical resolutions. The example makes use of the 15-micron CO2 band parameterization proposed by Akmaev and Shved (1982) for the LTE region of the middle atmosphere. Cooling rate errors provided by the present method are not more than 0.11 K/day for a vertical grid spacing of less than 5 km. The matrix transformation method is shown to be valuable for adapting effective matrix parameterizations of middle atmosphere cooling to atmospheric models with arbitrary vertical resolutions.

Akmaev, R. A.↗
NASA NTRS1990

Numerical modeling of the atmosphere with an isentropic vertical coordinate

A theta-coordinate model simulating the nonlinear evolution of a baroclinic wave is presented. In the model, vertical discretization maintains important integral constraints such as conservation of the angular momentum and total energy. A massless-layer approach is used in the treatment of the intersections of coordinate surfaces with the lower boundary. This formally eliminates the intersection problem, but raises other computational problems. Horizontal discretization of the continuity and momentum equations in the model are designed to overcome these problems. Selected results from a 10-day integration with the 25-layer, beta-plane version of the model are presented. It is concluded that the model can simulate the nonlinear evolution of a baroclinic wave and associated dynamical processes without major computational difficulties.

Hsu, Yueh-Jiuan G.↗
NASA NTRS1984

The effect of wind and moisture gradients on the arbitrary assignment of cloud motions to a vertical coordinate system in two Sesame cases

Satellite-derived cloud motion 'wind' vectors (CMV) are increasingly used in mesoscale and in global analyses, and questions have been raised regarding the uncertainty of the level assignment for the CMV. One of two major problems in selecting a level for the CMV is related to uncertainties in assigning the motion vector to either the cloud top or base. The second problem is related to the inability to transfer the 'wind' derived from the CMV at individually specified heights to a standard coordinated surface. The present investigation has the objective to determine if the arbitrary level assignment represents a serious obstacle to the use of cloud motion wind vectors in the mesoscale analysis of a severe storm environment.

Peslen, C. A.↗
NASA NTRS2005

Applications and Improvement of a Coupled, Global and Cloud-Resolving Modeling System

Recently Grabowski (2001) and Khairoutdinov and Randall (2001) have proposed the use of 2D CFWs as a "super parameterization" [or multi-scale modeling framework (MMF)] to represent cloud processes within atmospheric general circulation models (GCMs). In the MMF, a fine-resolution 2D CRM takes the place of the single-column parameterization used in conventional GCMs. A prototype Goddard MMF based on the 2D Goddard Cumulus Ensemble (GCE) model and the Goddard finite volume general circulation model (fvGCM) is now being developed. The prototype includes the fvGCM run at 2.50 x 20 horizontal resolution with 32 vertical layers from the surface to 1 mb and the 2D (x-z) GCE using 64 horizontal and 32 vertical grid points with 4 km horizontal resolution and a cyclic lateral boundary. The time step for the 2D GCE would be 15 seconds, and the fvGCM-GCE coupling frequency would be 30 minutes (i.e. the fvGCM physical time step). We have successfully developed an fvGCM-GCE coupler for this prototype. Because the vertical coordinate of the fvGCM (a terrain-following floating Lagrangian coordinate) is different from that of the GCE (a z coordinate), vertical interpolations between the two coordinates are needed in the coupler. In interpolating fields from the GCE to fvGCM, we use an existing fvGCM finite- volume piecewise parabolic mapping (PPM) algorithm, which conserves the mass, momentum, and total energy. A new finite-volume PPM algorithm, which conserves the mass, momentum and moist static energy in the z coordinate, is being developed for interpolating fields from the fvGCM to the GCE. In the meeting, we will discuss the major differences between the two MMFs (i.e., the CSU MMF and the Goddard MMF). We will also present performance and critical issues related to the MMFs. In addition, we will present multi-dimensional cloud datasets (i.e., a cloud data library) generated by the Goddard MMF that will be provided to the global modeling community to help improve the representation and performance of moist processes in climate models and to improve our understanding of cloud processes globally (the software tools needed to produce cloud statistics and to identify various types of clouds and cloud systems from both high-resolution satellite and model data will be also presented).

Tao, W.-K.↗
NASA NTRS1986

Vertical normal modes of a mesoscale model using a scaled height coordinate

Vertical modes were derived for a version of the Colorado State Regional Atmospheric Mesoscale Modeling System. The impacts of three options for dealing with the upper boundary of the model were studied. The standard model formulation holds pressure constant at a fixed altitude near the model top, and produces a fastest mode with a speed of about 90 m/sec. An alternative formulation, which allows for an external mode, could require recomputation of vertical modes for every surface elevation on the horizontal grid unless the modes are derived in a particular way. These results have bearing on the feasibility of applying vertical mode initialization to models with scaled height coordinates.

Lipton, A. E.↗
NASA NTRS1994

Measurement and interpretation of crustal deformation rates associated with postglacial rebound

Analysis of Global Positioning System (GPS) data from two sites separated by horizontal distance of only approximately 2.2 m yielded phase residuals exhibiting a systematic elevation angle dependence. One of the two GPS antennas was mounted on an approximately 1 m high concrete pillar, and the other was mounted on a standard wooden tripod. We performed elevation angle cutoff tests with these data, and established that the vertical coordinate of site position was sensitive to the minimum elevation angle (elevation cutoff) of the data analyzed. For example, the vertical coordinate of site position changed by 9.7 plus or minus 0.8 mm when the minimum elevation angle was increased from 10 to 25. We performed simulations based on a simple (ray tracing) multipath model with a single horizontal reflector, and demonstrated that the elevation angle cutoff test results and the pattern of the residual versus elevation angle could be qualitatively reproduced if the reflector were located 0.1-0.2 m beneath the antenna phase center. We therefore, hypothesized that the source of the elevation-angle-dependent error were multipath reflections and scattering and that the horizontal surface of the pillar, located a distance of approximately 0.2 m beneath the antenna phase center, was the primary reflector. We tested this hypothesis by placing microwave absorbing material between the antenna and the pillar in a number of configurations and analyzed the changes in apparent position of the antenna. The results indicate that (1) the horizontal surface of the pillar is indeed the main reflector, (2) both the concrete and the metal plate embedded in the pillar are significant reflectors, and (3) the reflection can be reduced to a great degree by the use of microwave absorbing materials. These results have significant implications for the accuracy of global GPS geodetic tracking networks which use pillar-antenna configuration identical or similar to the one used here (at the Westford WFRD GPS site).

Davis, James L.↗
NASA NTRS1985

Middle Atmosphere Reference Model Derived from Satellite Data

Zonally averaged temperature and geostrophic zonal wind for each month and for latitudes from 80 S to 80 N with pressure scale height as a vertical coordinate are given tabular form. The pressure scale height is defined as -ln (p/P sub O) where p is pressure and P sub O is surface pressure. The height interval corresponding to one pressure scale height is proportional to absolute temperature, and is 7 km at 240 K. Values in are given at intervals of 0.5 pressure scale heights, i.e., approximately 3,5 km. The geopotential height fields for the principal seasons are given. Temperature, pressure and density with geometric height as the vertical coordinate at invervals of 5 km. The temperature values were obtained from a combination of satellite data above 30 mb with values supplied by the Berlin Free University at 30 mb and the climatology derived by OORT (1983) for 50 mb and below. The geopotential height fields were obtained from these temperature fields by integrating up and down from the 30 mb geopotential height supplied by Berlin Free University. The geostrophic winds were obtained by differentiating these geopotential height fields.

Barnett, J. J.↗
NASA NTRS2014

Natural Ocean Carbon Cycle Sensitivity to Parameterizations of the Recycling in a Climate Model

Sensitivities of the oceanic biological pump within the GISS (Goddard Institute for Space Studies ) climate modeling system are explored here. Results are presented from twin control simulations of the air-sea CO2 gas exchange using two different ocean models coupled to the same atmosphere. The two ocean models (Russell ocean model and Hybrid Coordinate Ocean Model, HYCOM) use different vertical coordinate systems, and therefore different representations of column physics. Both variants of the GISS climate model are coupled to the same ocean biogeochemistry module (the NASA Ocean Biogeochemistry Model, NOBM), which computes prognostic distributions for biotic and abiotic fields that influence the air-sea flux of CO2 and the deep ocean carbon transport and storage. In particular, the model differences due to remineralization rate changes are compared to differences attributed to physical processes modeled differently in the two ocean models such as ventilation, mixing, eddy stirring and vertical advection. GISSEH(GISSER) is found to underestimate mixed layer depth compared to observations by about 55% (10 %) in the Southern Ocean and overestimate it by about 17% (underestimate by 2%) in the northern high latitudes. Everywhere else in the global ocean, the two models underestimate the surface mixing by about 12-34 %, which prevents deep nutrients from reaching the surface and promoting primary production there. Consequently, carbon export is reduced because of reduced production at the surface. Furthermore, carbon export is particularly sensitive to remineralization rate changes in the frontal regions of the subtropical gyres and at the Equator and this sensitivity in the model is much higher than the sensitivity to physical processes such as vertical mixing, vertical advection and mesoscale eddy transport. At depth, GISSER, which has a significant warm bias, remineralizes nutrients and carbon faster thereby producing more nutrients and carbon at depth, which eventually resurfaces with the global thermohaline circulation especially in the Southern Ocean. Because of the reduced primary production and carbon export in GISSEH compared to GISSER, the biological pump efficiency, i.e., the ratio of primary production and carbon export at 75 m, is half in the GISSEH of that in GISSER, The Southern Ocean emerges as a key region where the CO2 flux is as sensitive to biological parameterizations as it is to physical parameterizations. The fidelity of ocean mixing in the Southern Ocean compared to observations is shown to be a good indicator of the magnitude of the biological pump efficiency regardless of physical model choice.

Antarctic Ocean↗
NASA NTRS1985

Development and Testing of the Variable Vertical Resolution Fourth Order GCM

The vertical coordinate of the Fourth Order Model has been generalized so that the model can now run with an arbitrary number of vertical layers and so that the thicknesses of these layers can be arbitrarily specified (in the sigma coordinate). This Variable Vertical Resolution (VVR) version of the Fourth Order Model will soon replace the current production model. To assess the skill of the VVR model, it has been run with 9 equally spaced layers and compared with the current production model. In two Northern Hemispheric winter cases and one summer case, the two models were virtually identical in forecast skill for 6 to 7 days. After that the VVR model was slightly better in the winter cases and the production model was slightly better in the summer case. The only exception to this was that after 2 days the production model gave slightly more skillful 500 mb forecasts in the tropics for the summer case.

Helfand, H. M.↗
NASA NTRS1997

Wave Driven Exothermic Heating in the Mesopause Region

A full-wave propagation model was developed that describes the propagation of gravity waves from the Earth's surface to the upper boundary, which can be placed anywhere between 150 and 500 km altitude. The model includes a realistic background atmosphere, and includes the effects of mean horizontal winds and their vertical shears, mean vertical temperature gradients, the eddy and molecular diffusion of heat and momentum, and the effects of ion-drag. This model solves five coupled second-order differential equations (continuity, momentum, and energy) in the vertical coordinate to derive the perturbation variables u', v', w' (horizontal and vertical velocity components), T' (temperature) and p' (pressure). The upper boundary can be automatically selected based on tests using the radiation condition at the upper boundary, wherein the height is increased until the wave is experiencing severe dissipation at the upper boundary, ensuring that substantial absorption occurs for any waves reflected from the upper boundary. The determination of wave amplitude is a key requirement of wave energetics. Therefore, the fullwave model has been applied to airglow observations in order to determine wave amplitudes as a function of altitude. This was accomplished by using the full-wave model output to drive a chemistry perturbation module that describes minor species perturbations and the resulting airglow perturbations. The full-wave output was multiplied by an altitude-independent factor such that the modeled and observed relative airglow intensity perturbations were equal. The effects of mean winds were included in these studies, and found to be the most important model input affecting the calculations (being more important than the choice of eddy diffusion profiles and chemical kinetic coefficients). In one study (Hickey et al., 1997a) these winds could not be well estimated from the measurements, whereas in the second study (Hickey et al.,1997b) the mean were well defined with a sodium wind-temperature lidar.

Hickey, Michael P.↗