SEARCH · Engineering Papers

Results for “Primitive equations”

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.

SourceResearch connectionFrom yearThrough year

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

At least 199 records · Page 11

All sources NASA NTRS DOE OSTI

NASA NTRS1992

Simulations of the seasonal evolution of planetary waves in the stratosphere

Results are presented of experiments performed with the 3D primitive equation model of the middle atmosphere to elucidate the connections between the behavior of planetary waves in the stratosphere and in the upper troposphere during the seasonal cycle. In these experiments the geopotential height fields at the model's lower boundary consisted of a time-varying zonally symmetric component plus a time-independent wave 1. Experiments were also performed with the same time-varying zonally symmetric component plus a time-independent combination of wave 1 and 2. In each of the wave 1-only experiments, the wave 1 amplitude in the middle stratosphere has two maxima, one in fall and the other in spring, separated by a midwinter minimum. It is suggested that some important transient features in the seasonal evolution of stratospheric planetary waves do not depend on time variations in the amplitude or phase of the corresponding tropospheric waves.

Farrara, John D.↗

NASA NTRS1992

The coupling of the troposphere and stratosphere

A linear, time-dependent, primitive-equation, 3D numerical model is used to investigate the dynamical coupling between the troposphere and stratosphere through the vertical propagation of planetry waves. Attention is given to the roles of the tropopause in affecting the propagation of planetary waves and to the extent to which the stratosphere can determine its own wave driving, i.e., the sensitivity of the generation of the wave activity in the troposphere to the change of zonal wind in the stratosphere. Two maxima in Eliassen-Palm (EP) flux convergence are found for wave 1 in high latitudes, one in the upper stratosphere and the other in the upper troposphere. For wave 2 there is a third maximum north of the zero wind line in the tropical stratosphere. The formation of the maximum in EP-flux convergence in the upper troposphere is related to the dissipation of non-WKBJ waves near a local maximum of the refractive index. The existence of the local maximum in refractive index at the tropopause is a result of the great gradient of the buoyancy frequency as well as the vertical shear of the zonal flow there.

Chen, Ping↗

NASA NTRS1992

Characteristics of 3-D transport simulations of the stratosphere and mesosphere

A 3D mechanistic, primitive-equation model of the stratosphere and mesosphere is coupled to an offline spectral transport model. The dynamics model is initialized with and forced by observations so that the coupled models may be used to study specific episodes. Results are compared with those obtained by transport online in the dynamics model. Although some differences are apparent, the results suggest that coupling of the models to a comprehensive photochemical package will provide a useful tool for studying the evolution of constituents in the middle atmosphere during specific episodes.

Fairlie, T. D. A.↗

NASA NTRS1993

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.↗

NASA NTRS1993

A diagnostic model to estimate winds and small-scale drag from Mars Observer PMIRR data

Theoretical and modeling studies indicate that small-scale drag due to breaking gravity waves is likely to be of considerable importance for the circulation in the middle atmospheric region (approximately 40-100 km altitude) on Mars. Recent earth-based spectroscopic observations have provided evidence for the existence of circulation features, in particular, a warm winter polar region, associated with gravity wave drag. Since the Mars Observer PMIRR experiment will obtain temperature profiles extending from the surface up to about 80 km altitude, it will be extensively sampling middle atmospheric regions in which gravity wave drag may play a dominant role. Estimating the drag then becomes crucial to the estimation of the atmospheric winds from the PMIRR-observed temperatures. An interative diagnostic model based upon one previously developed and tested with earth satellite temperature data will be applied to the PMIRR measurements to produce estimates of the small-scale zonal drag and three-dimensional wind fields in the Mars middle atmosphere. This model is based on the primitive equations, and can allow for time dependence (the time tendencies used may be based upon those computed in a Fast Fourier Mapping procedure). The small-scale zonal drag is estimated as the residual in the zonal momentum equation; the horizontal winds having first been estimated from the meridional momentum equation and the continuity equation. The scheme estimates the vertical motions from the thermodynamic equation, and thus needs estimates of the diabatic heating based upon the observed temperatures. The latter will be generated using a radiative model. It is hoped that the diagnostic scheme will be able to produce good estimates of the zonal gravity wave drag in the Mars middle atmosphere, estimates that can then be used in other diagnostic or assimilation efforts, as well as more theoretical studies.

Barnes, J. R.↗

NASA NTRS1993

The annual pressure cycle on Mars: Results from the LMD Martian atmospheric general circulation model

We have been developing a General Circulation Model (GCM) of the martian atmosphere since 1989. The model has been described rather extensively elsewhere and only the main characteristics are given here. The dynamical part of the model, adapted from the LMD terrestrial climate model, is based on a finite-difference formulation of the classical 'primitive equations of meteorology.' The radiative transfer code includes absorption and emission by CO2 (carefully validated by comparison to line-by-line calculations) and dust in the thermal range and absorption and scattering by dust in the visible range. Other physical parameterizations are included: modeling of vertical turbulent mixing, dry convective adjustment (in order to prevent vertical unstable temperature profiles), and a multilayer model of the thermal conduction in the soil. Finally, the condensation-sublimation of CO2 is introduced through specification of a pressure-dependent condensation temperature. The atmospheric and surface temperatures are prevented from falling below this critical temperature by condensation and direct precipitation onto the surface of atmospheric CO2. The only prespecified spatial fields are the surface thermal inertia, albedo, and topography.

Hourdin, Frederic↗

NASA NTRS1993

The influence of orography on the transport of atmospheric constituents

Given the fact that martian orography varies over a range of two scale heights (cf., about one scale height on Earth), it is thought that on Mars some aspects of the circulation should be affected to a large degree by the presence of orography. Initial results from two numerical models are presented that illustrate the effect of orography on cross-equatorial transport on Mars. On Earth, western boundary currents are found not only in the oceans as one would expect, but also in the troposphere, the most notable example being the East African Jet, which has an important effect on the monsoon circulation of the Eastern Indian Ocean. Two models were used for this investigation: a simple GCM (hereafter referred to as the SGCM); and a simple one-level barotropic model (BM). The SGCM solves the hydrostatic primitive equations on a sphere by representing model variables by spherical harmonics in the horizontal, with nonlinear terms being calculated in grid-point space. In the vertical, the model uses the sigma coordinate system in finite-difference form. The model represents the radiative processes in the atmosphere by a Newtonian relaxation to a specified zonal mean temperature state and surface drag by Rayleigh friction in the lowest model level.

Joshi, Manoj M.↗

NASA NTRS1993

An observing system simulation experiment for the Laser Atmospheric Wind Sounder (LAWS)

The present observing system-simulation experiments evaluate the potential of a Laser Atmospheric Wind Sounder (LAWS) instrument for 5-day forecasting, using a primitive-equation multilevel spectral global circulation model. A 55-deg-inclined and a 98-deg sun-synchronous orbit are examined, by adding LAWS wind profiles into a global 4D data-assimilation system, and comparing both the analyses and forecasts to a control experiment. The 4D data-assimilation system consists of a multivariate optimum interpolation analysis and a nonlinear, normal-mode intialization, using the aforementioned global circulation model.

Rohaly, G. D.↗

NASA NTRS1994

Experiments with the Mesoscale Atmospheric Simulation System (MASS) using the synthetic relative humidity

This study is intended to examine the impact of the synthetic relative humidity on the model simulation of mesoscale convective storm environment. The synthetic relative humidity is derived from the National Weather Services surface observations, and non-conventional sources including aircraft, radar, and satellite observations. The latter sources provide the mesoscale data of very high spatial and temporal resolution. The synthetic humidity data is used to complement the National Weather Services rawinsonde observations. It is believed that a realistic representation of initial moisture field in a mesoscale model is critical for the model simulation of thunderstorm development, and the formation of non-convective clouds as well as their effects on the surface energy budget. The impact will be investigated based on a real-data case study using the mesoscale atmospheric simulation system developed by Mesoscale Environmental Simulations Operations, Inc. The mesoscale atmospheric simulation system consists of objective analysis and initialization codes, and the coarse-mesh and fine-mesh dynamic prediction models. Both models are a three dimensional, primitive equation model containing the essential moist physics for simulating and forecasting mesoscale convective processes in the atmosphere. The modeling system is currently implemented at the Applied Meteorology Unit, Kennedy Space Center. Two procedures involving the synthetic relative humidity to define the model initial moisture fields are considered. It is proposed to perform several short-range (approximately 6 hours) comparative coarse-mesh simulation experiments with and without the synthetic data. They are aimed at revealing the model sensitivities should allow us both to refine the specification of the observational requirements, and to develop more accurate and efficient objective analysis schemes. The goal is to advance the MASS (Mesoscal Atmospheric Simulation System) modeling expertise so that the model output can provide reliable guidance for thunderstorm forecasting.

Chang, Chia-Bo↗

NASA NTRS1994

Evidence for eddy formation in the eastern Arabian Sea during the northeast monsoon

The seasonal formation of a large (500-800 km diameter) anticyclonic eddy in the upper 300-400 m of the eastern Arabian Sea during the northest monsoon period (December-April) is indicated fom hydrographic and satellite altimetry sea level observations, as well as from numerical model experiments. The center of the eddy circulation is approximately 10 deg N, 70 deg E, just to the west of the north-south Laccadive Island chain. In this paper the eddy is called the Laccadive High (LH). In some ways it is like a mirrorlike counterpart to the Great Whirl that develops during the southwest monsoon of the Somali coast (western Arabian Sea). The LH occurs at the same latitude but on the opposite side of the basin during the reversed monsoon. It is different from the Great Whirl, however, in its formation process, its intensity, and its decay. The hydrographic data obtained from surveys all during a single season give sufficiently close station spacing to allow reasonable contouring of the geopotential surfaces and of the properties within and around the LH region with minimum time aliasing. The Geostat altimeter record extends over 4 years, during which the seasonal variability of the LH indicates a dynamic relief of approximately 15-20 cm, which is in good agreement with the hydrographics observations. The altimetry time series also suggests a westward translation of the LH by January with a subsequent dissipation in midbasin. The model used is a wind-forced three-layer primitive equation model which depicts a LH agreement with the timing, position, and amplitude of both the hydrographic and altimetric measurements. The numerical simulation includes a passive tracer located in the Western Bay of Bengal; the western advection of the tracer around the south coasts of Sri Lanka and India in December and January is consistent with the appearance of low-salinity water observed to extend into the Arabian Sea during this period. The modeling studies suggest that both local and remote forcing are important in formation of the LH.

Bruce, John G.↗

NASA NTRS1994

Cyclogenesis in a saturated environment

The dynamics of baroclinic wave growth in a saturated environment is examined using linear and nonlinear models employing a parameterization of latent heat release that assumes all rising air is saturated, and saturation equivalent potential temperature is conserved on ascent. Piecewise potential vorticity (PV) diagnostics are used to interpret the results. When the stability to vertical displacements in saturated air is allowed to increase with height, as it must in an atmosphere with a constant, positive lapse rate of potential temperature, the growth rates of the most unstable modes of the Eady problem grow only marginally faster than the modes of the dry problem. The vertical variation of moist static stability produces a gradient of moist potential vorticity in the rising air, eliminating the short-wave cutoff present in the dry Eady problem. The destabilization of the short waves is shown to be associated with the interaction between surface potential temperature anomalies and diabatically generated lower-tropospheric potential vorticity anomalies. Nonlinear primitive equation simulations, starting from normal-mode initial conditions, show that while the dry wave grows at nearly the linear growth rate until maximum amplitude is reached, the moist wave grows significantly faster than the linear growth rate at finite amplitude. This enhanced growth is associated with the rapid amplification of a mesoscale PV anomaly generated by latent heat release at the warm front. The rapid amplification of the surface cyclone results from the superposition of the circulation associated with this misoscale PV anomaly upon the circulation associated with the surface and upper boundary potential temperature anomalies. Additional integrations with finite-amplitude initial conditions more typical of atmospheric conditions exhibit similar behavior. It is suggested that many of the rapid cyclogenesis events that occur as upper-troposheric PV anomalies cross the east coasts of continents may arise from the rapid generation of PV anomalies by condensational heating in the moist maritime lower troposphere.

Whitaker, Jeffrey S.↗

NASA NTRS1994

The simulated Indian monsoon: A GCM sensitivity study

A series of sensitivity experiments are conducted in an attempt to understand and correct deficiencies in the simulation of the seasonal mean Indian monsoon with a global atmospheric general circulation model. The seasonal mean precipitation is less than half that observed. This poor simulation in seasonal integrations is independent of the choice of initial conditions and global sea surface temperature data used. Experiments are performed to test the sensitivity of the Indian monsoon simulation to changes in orography, vegetation, soil, wetness, and cloudiness. The authors find that the deficiency of the model precipitation simulation may be attributed to the use of an enhanced orography in the integrations. Replacement of this orography with a mean orography results in a much more realistic simulation of Indian monsoon circulation and rainfall. Experiments with a linear primitive equation model on the sphere suggest that this striking improvement is due to modulations of the orographically forced waves in the lower troposphere. This improvement in the monsoon simulation is due to the kinematic and dynamical effects of changing the topography, rather than the thermal effects, which were minimal. The magnitude of the impact on the Indian monsoon of the other sensitivity experiments varied considerably, but was consistently less than the impact of using the mean orography. However, results from the soil moisture sensitivity experiments suggest a possibly important role for soil moisture in simulating tropical precipitation, including that associated with the Indian monsoon.

Fennessy, M. J.↗

NASA NTRS1994

Planetary-scale circulations in the presence of climatological and wave-induced heating

Interaction between the large-scale circulation and the convective pattern is investigated in a coupled system governed by the linearized primitive equations. Convection is represented in terms of two components of heating: A 'climatological component' is prescribed stochastically to represent convection that is maintained by fixed distributions of land and sea and sea surface temperature (SST). An 'induced component' is defined in terms of the column-integrated moisture flux convergence to represent convection that is produced through feedback with the circulation. Each component describes the envelope organizing mesoscale convective activity. As SST on the equator is increased, induced heating amplifies in the gravest zonal wavenumbers at eastward frequencies, where positive feedback offsets dissipation. Under barotropic stratification, a critical SST of 29.5 C results in positive feedback exactly cancelling dissipation in wavenumber 1 for an eastward phase speed of 6 m/s. Sympathetic interaction between the circulation and the induced heating is the basis for 'frictional wave-Conditional Instability of the Second Kind (CISK)', which is distinguished from classical wave-CISK by rendering the gravest zonal dimensions most unstable. Under baroclinic stratification, the coupled system exhibits similar behavior. The critical SST is only 26.5 C for conditions representative of equinox, but in excess of 30 C for conditions representative of solstice. Having the form of an unsteady Walker circulation, the disturbance produced by frictional wave-CISK compares favorably with the observed life cycle of the Madden-Julian oscillation (MJO). SST above the critical value produces an amplifying disturbance in which enhanced convection coincides with upper-tropospheric westerlies and is positively correlated with temperature and surface convergence. Conversely, SST below the critical value produces a decaying disturbance in which enhanced convection coincides with upper-tropospheric easterlies and is nearly in quadrature with temperature and surface convergence. While sharing essential features with the MJO in the Eastern Hemisphere, frictional wave-CISK does not explain observed behavior in the Western Hemisphere, where the convective signal is largely absent. Comprised of Kelvin structure with the same frequency, observed behavior in the Western Hemisphere can be understood as a propagating response that is excited in and radiates away from the fluctuation of convection in the Eastern Hemisphere.

Salby, Murry L↗

NASA NTRS1994

Equatorial wave activity derived from fluctuations in observed convection

The spectrum of equatorial wave activity propagating vertically into the stratosphere is calculated from high-resolution imagery of the global convective pattern. Synoptic Global Cloud Imagery (GCI), constructed from six satellites simultaneously observing the earth, is used to diabatically force the linearized primitive equations. Having resolution of 0.5 deg and 3 h, that imagery captures the dominant scales of organized convection, including several harmonics of the diurnal cycle. Its global coverage with high space-time resolution allows the GCI to represent heating variability and dynamical behavior excited by it over a wide range of scales. The dynamical response above the heating is evaluated globally in terms of a space-time spectrum of Hough modes, one which includes planetary-scale Kelvin waves, Rossby waves, and gravity waves down to the resolution of the GCI. The geopotential response, which is indicative of temperature fluctuations observed by satellite, is very red in frequency. Therefore, planetary-scale waves with periods longer than two days dominate the spectrum of geopotential, while high-frequency gravity waves make a comparatively small contribution. Some 80% of the geopotential variance is accounted for by the Kelvin and gravest-symmetric Rossby modes, while the Rossby-gravity mode is comparatively weak. In horizontal eddy motion, the excited wave spectrum is still dominated by planetary-scale components. However, meridional wind fluctuations associated with the Rossby-gravity mode have variance comparable to that of zonal wind fluctuations associated with the Kelvin mode, even though the Rossby-gravity mode is nearly invisible in the geopotential response. Estimates of tropospheric heating lead to amplitudes and propagation characteristics that are broadly consistent with satellite and radiosonde observations of wave activity in the lower stratosphere. The space-time spectrum of EP flux is significantly whiter than the response in either geopotential or motion. Gravity waves of small scale and high frequency carry a large fraction of the upward flux. Although it dominates eastward variance of geopotential and motion, the Kelvin mode carries only about 50% of the eastward EP flux at phase speeds of 20-40 m/s and only 35% of the total eastward flux transmitted to the stratosphere. The remainder is carried by the gravity wave spectrum, which carries nearly all of the westward flux at phase speeds greater than 20 m/s. The gravity wave spectrum also contributes significantly at phase speeds of 10-20 m/s, where only 25% of the flux is accounted for by zonal wavenumbers less than 20. The broad nature of the gravity wave spectrum suggests its absorption at critical levels will be distributed over a deep layer of the middle atmosphere.

Bergman, John W.↗

NASA NTRS1994

Simulations of the February 1979 stratospheric sudden warming: Model comparisons and three-dimensional evolution

The evolution of the stratopsheric flow during the major stratospheric sudden warming of February 1979 is studied using two primitive equation models of the stratosphere and mesosphere. The United Kingdom Meteorological Office Stratosphere-Mesosphere Model (SMM) uses log pressure as a vertical coordinate. A spectral, entropy coordinate version of the SMM (entropy coordinate model, or ECM) that has recently been developed is also used. Comparison of SMM simulations with forecasts performed using the University of California, Los Angeles general circulation model confirms the previously noted sensitivity of stratospheric forecasts to tropospheric forecasts and emphasizes the importance of adequate vertical resolution in modeling the stratosphere. The ECM simulations provide a schematic description of the three-dimensional evolution of the polar vortex and the motion of air through it. During the warming, the two cyclonic vortices tilt westward and equatorward with height. Strong upward velocities develop in the lower stratosphere on the west (cold) side of a baroclinic zone as it forms over Europe and Asia. Strong downward velocities appear in the upper stratosphere on the east (warm) side, strengthening the temperature gradients. After the peak of the warming, vertical velocities decrease, downward velocities move into the lower stratosphere, and upward velocities move into the upper stratosphere.

Manney, G. L.↗

NASA NTRS1993

A two-dimensional model with coupled dynamics, radiation, and photochemistry. 1: Simulation of the middle atmosphere

Spatial and seasonal distribution of ozone and other trace gases are simulated using a interactive two-dimensional model of the stratosphere updated to include full gas phase chemistry. The model consists of a primitive equation dynamics module, a full radiative transfer scheme, and a comprehensive gas phase chemistry module. The circulation is derived from heating rates in the stratosphere that are calculated using model-generated ozone. In the troposphere, parameterized heating rates are adopted. The eddy momentum flux divergence in the zonal mean momentum equation is given by the eddy fluxes of potential vorticity. Eddy fluxes of potential vorticity and tracers are parameterized using a set of predetermined diffusion coefficients. The adopted values for K(sub yy), show a hemispherical asymmetry in that the values in the lower stratosphere are consistently smaller in the southern hemisphere. The asymmetry in K(sub yy) and in the parameterization of the tropospheric heating rate results in an asymmetry in the circulation giving rise to unique signatures in the trace gas distributions. The model successfully simulates the observed asymmetry in the column abundance of the springtime ozone maxima between northern and southern hemisphere. Results for other trace gases are in agreement with the gross observed features although specific differences exist.

Ko, Malcolm K. W.↗

NASA NTRS1995

Short-range precipitation forecasts using assimilation of simulated satellite water vapor profiles and column cloud liquid water amounts

These observing system simulation experiments investigate the assimilation of satellite-observed water vapor and cloud liquid water data in the initialization of a limited-area primitive equations model with the goal of improving short-range precipitation forecasts. The assimilation procedure presented includes two aspects: specification of an initial cloud liquid water vertical distribution and diabatic initialization. The satellite data is simulated for the next generation of polar-orbiting satellite instruments, the Advanced Microwave Sounding Unit (AMSU) and the High-Resolution Infrared Sounder (HIRS), which are scheduled to be launched on the NOAA-K satellite in the mid-1990s. Based on cloud-top height and total column cloud liquid water amounts simulated for satellite data a diagnostic method is used to specify an initial cloud water vertical distribution and to modify the initial moisture distribution in cloudy areas. Using a diabatic initialization procedure, the associated latent heating profiles are directly assimilated into the numerical model. The initial heating is estimated by time averaging the latent heat release from convective and large-scale condensation during the early forecast stage after insertion of satellite-observed temperature, water vapor, and cloud water formation. The assimilation of satellite-observed moisture and cloud water, together withy three-mode diabatic initialization, significantly alleviates the model precipitation spinup problem, especially in the first 3 h of the forecast. Experimental forecasts indicate that the impact of satellite-observed temperature and water vapor profiles and cloud water alone in the initialization procedure shortens the spinup time for precipitation rates by 1-2 h and for regeneration of the areal coverage by 3 h. The diabatic initialization further reduces the precipitation spinup time (compared to adiabatic initialization) by 1 h.

Wu, Xiaohua↗

NASA NTRS1995

A steady-state Kalman filter for assimilating data from a single polar orbiting satellite

A steady-state scheme for data assimilation in the context of a single, short period (relative to a day), sun-synchronous, polar-orbiting satellite is examined. If the satellite takes observations continuously, the gains, which are the weights for blending observations and predictions together, are steady in time. For a linear system forced by random noise, the optimal steady-state gains (Wiener gains) are equivalent to those of a Kalman filter. Computing the Kalman gains increases the computational cost of the model by a large factor, but computing the Wiener gains does not. The latter are computed by iteration using prior estimates of the gains to assimilate simulated observations of one run of the model, termed 'truth' into another run termed 'prediction'. At each stage, the prediction errors form the basis for the next estimate of the gains. Steady state is achieved after three or four iterations. Further simplification is achieved by making the gains depend on longitudinal distance from the observation point, not on absolute longitude. For a single-layer primitive equation model, the scheme works well even if only the mass field is observed but not the velocity field. Although the scheme was developed for Mars Observer, it should be applicable to data retrieved from Earth atmosphere satellites, for example, UARS.

Banfield, Don↗