Engineering Papers⌕ Search

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.

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

At least 217 records · Page 12

An approximate Kalman filter for ocean data assimilation: An example with an idealized Gulf Stream model

A practical method of data assimilation for use with large, nonlinear, ocean general circulation models is explored. A Kalman filter based on approximation of the state error covariance matrix is presented, employing a reduction of the effective model dimension, the error's asymptotic steady state limit, and a time-invariant linearization of the dynamic model for the error integration. The approximations lead to dramatic computational savings in applying estimation theory to large complex systems. We examine the utility of the approximate filter in assimilating different measurement types using a twin experiment of an idealized Gulf Stream. A nonlinear primitive equation model of an unstable east-west jet is studied with a state dimension exceeding 170,000 elements. Assimilation of various pseudomeasurements are examined, including velocity, density, and volume transport at localized arrays and realistic distributions of satellite altimetry and acoustic tomography observations. Results are compared in terms of their effects on the accuracies of the estimation. The approximate filter is shown to outperform an empirical nudging scheme used in a previous study. The examples demonstrate that useful approximate estimation errors can be computed in a practical manner for general circulation models.

Fukumori, Ichiro↗

Observational and numerical studies of extreme frontal scale contraction

The general objective of this effort is to increase understanding of how frontal scale contraction processes may create and sustain intense mesoscale precipitation along intensifying cold fronts. The five-part project (an expansion of the originally proposed two-part project) employed conventional meteorological data, special mesoscale data, remote sensing measurements, and various numerical models. First an idealized hydrostatic modeling study of the scale contraction effects of differential cloud cover on low-level frontal structure and dynamics was completed and published in a peer-reviewed journal. The second objective was to complete and publish the results from a three dimensional numerical model simulation of a cold front in which differential sensible heating related to cloud coverage patterns was apparently crucial in the formation of a severe frontal squall line. The third objective was to use a nonhydrostatic model to examine the nonlinear interactions between the transverse circulation arising from inhomogeneous cloud cover, the adiabatic frontal circulation related to semi-geostrophic forcing, and diabatic effects related to precipitation processes, in the development of a density current-like microstructure at the leading edge of cold fronts. Although the development of a frontal model that could be used to initialize such a primitive equation model was begun, we decided to focus our efforts instead on a project that could be successfully completed in this short time, due to the lack of prospects for continued NASA funding beyond this first year (our proposal was not accepted for future funding). Thus, a fourth task was added, which was to use the nonhydrostatic model to test tentative hypotheses developed from the most detailed observations ever obtained on a density current (primarily sodar and wind profiler data). These simulations were successfully completed, the findings were reported at a scientific conference, and the results have recently been submitted to a peer-reviewed journal. The fifth objective was to complete the analysis of data collected during the Cooperative Oklahoma Profiler Studies (COPS-91) field project, which was supported by NASA. The analysis of the mesoscale surface and sounding data, Doppler radar imagery, and other remote sensing data from multi frequency wind profiler, microwave radiometer, and the Radio Acoustic Sounding System has been completed. This study is a unique investigation of processes that caused the contraction of a cold front to a microscale zone exhibiting an undular bore-like structure. Results were reported at a scientific conference and are being prepared for publication. In summary, considerable progress has been achieved under NASA funding in furthering our understanding of frontal scale contraction and density current - gravity wave interaction processes, and in utilizing models and remotely sensed data in such studies.

Koch, Steven E.↗

Dynamics of the Antarctic Circumpolar Current. Evidence for Topographic Effects from Altimeter Data and Numerical Model Output

Geosat altimeter data and numerical model output are used to examine the circulation and dynamics of the Antarctic Circumpolar Current (ACC). The mean sea surface height across the ACC has been reconstructed from height variability measured by the altimeter, without assuming prior knowledge of the geoid. The results indicate locations for the Subantarctic and Polar Fronts which are consistent with in situ observations and indicate that the fronts are substantially steered by bathymetry. Detailed examination of spatial and temporal variability indicates a spatial decorrelation scale of 85 km and a temporal e-folding scale of 34 days. Empirical Orthogonal Function analysis suggests that the scales of motion are relatively short, occuring on 1000 km length-scales rather than basin or global scales. The momentum balance of the ACC has been investigated using output from the high resolution primitive equation model in combination with altimeter data. In the Semtner-Chervin quarter-degree general circulation model topographic form stress is the dominant process balancing the surface wind forcing. In stream coordinates, the dominant effect transporting momentum across the ACC is bibarmonic friction. Potential vorticity is considered on Montgomery streamlines in the model output and along surface streamlines in model and altimeter data. (AN)

Gille, Sarah T.↗

Intraseasonal and Interannual Variability of Mars Present Climate

This is a Final Report for a Joint Research Interchange (JRI) between NASA Ames Research Center and San Jose State University, Department of Meteorology. The focus of this JRI has been to investigate the nature of intraseasonal and interannual variability of Mars'present climate. We have applied a three-dimensional climate model based on the full hydrostatic primitive equations to determine the spatial, but primarily, the temporal structures of the planet's large-scale circulation as it evolves during a given seasonal advance, and, over multi-annual cycles. The particular climate model applies simplified physical parameterizations and is computationally efficient. It could thus easily be integrated in a perpetual season or advancing season configuration, as well as over many Mars years. We have assessed both high and low-frequency components of the circulation (i.e., motions having periods of Omicron(2-10 days) or greater than Omicron(10 days), respectively). Results from this investigation have explored the basic issue whether Mars' climate system is naturally 'chaotic' associated with nonlinear interactions of the large-scale circulation-regardless of any allowance for year-to-year variations in external forcing mechanisms. Titles of papers presented at scientific conferences and a manuscript to be submitted to the scientific literature are provided. An overview of a areas for further investigation is also presented.

Hollingsworth, Jeffery L.↗

Utilization of STRAT Data

One of the goals of the STRAT mission was to improve our understanding of transport in the stratosphere. We have used our coupled 2-D model to analyse transport between the tropics and midlatitudes using aircraft measurements of (CO2) and (N2O) made during the SPADE and ASHOE/MAESA campaigns (Boering et al. 1994) for model validation. We have further used the model to estimate errors in the experimental determination of the age of stratospheric air that arise from nonlinear increases of real tracers such as (CO2) and (SF6). The interpretation of measurements made during STRAT is still in progress. Most of the work performed under this contract has been described in the second year progress report, submitted in April of 1996. In the following, the research is summarized and some additional work that has become relevant in the meantime is described. We have developed an interactive model of the dynamics, chemistry and radiation of the stratosphere with partial support from this contract. The dynamics module integrates the primitive equations on a sphere. Lower boundary zonal winds and eddy geopotential heights are specified from observations. Chemistry and tracer transport are done in two dimensions. The chemistry model is a family model. Long lived species are transported and short lived species are determined in intervals of 5 to 10 days by assuming photochemical steady state. Full diurnal integrations with the fast species are performed periodically and diurnal coefficients, derived from the explicit integration are used to determine the steady state solution. Heating rates in the stratosphere are calculated from model ozone and the model temperatures and circulation fields are used in the chemistry transport module.

Source record↗

Observational and Modeling Studies of Radiative, Chemical, and Dynamical Interactions in the Earth''s Atmosphere

A 3-dimensional model was developed to support mechanistic studies. The model solves the global primitive equations in isentropic coordinates, which directly characterize diabatic processes forcing the Brewer-Dobson circulation of the middle atmosphere. It's numerical formulation is based on Hough harmonics, which partition horizontal motion into its rotational and divergent components. These computational features, along with others, enable 3D integrations to be performed practically on RISC computer architecture, on which they can be iterated to support mechanistic studies. The model conserves potential vorticity quite accurately under adiabatic conditions. Forced by observed tropospheric structure, in which integrations are anchored, the model generates a diabatic circulation that is consistent with satellite observations of tracer behavior and diabatic cooling rates. The model includes a basic but fairly complete treatment of gas-phase photochemistry that represents some 20 chemical species and 50 governing reactions with diurnally-varying shortwave absorption. The model thus provides a reliable framework to study transport and underlying diabatic processes, which can then be compared against chemical and dynamical structure observed and in GCM integrations. Integrations with the Langley GCM were performed to diagnose feedback between simulated convection and the tropical circulation. These were studied in relation to tropospheric properties controlling moisture convergence and environmental conditions supporting deep convection, for comparison against mechanistic integrations of wave CISK that successfully reproduce the Madden-Julian Oscillation (MJO) of the tropical circulation. These comparisons were aimed at identifying and ultimately improving aspects of the convective simulation, with the objective of recovering a successful simulation of the MJO in the Langley GCM, behavior that should be important to budgets of upper-tropospheric water vapor and chemical species.

Salby, Murry↗

Modeling the Gulf Stream System: How Far from Reality?

Analyses of a primitive equation ocean model simulation of the Atlantic Ocean circulation at 1/6 deg horizontal resolution are presented with a focus on the Gulf Stream region. Among many successful features of this simulation, this letter describes the Gulf Stream separation from the coast of North America near Cape Hatteras, meandering of the Gulf Stream between Cape Hatteras and the Grand Banks, and the vertical structure of temperature and velocity associated with the Gulf Stream. These results demonstrate significant improvement in modeling the Gulf Stream system using basin- to global scale ocean general circulation models. Possible reasons responsible for the realistic Gulf Stream simulation are discussed, contrasting the major differences between the present model configuration and those of previous eddy resolving studies.

Choa, Yi↗

A Preliminary Model Study of the Large-Scale Seasonal Cycle in Bottom Pressure Over the Global Ocean

Output from the primitive equation model of Semtner and Chervin is used to examine the seasonal cycle in bottom pressure (Pb) over the global ocean. Effects of the volume-conserving formulation of the model on the calculation Of Pb are considered. The estimated seasonal, large-scale Pb signals have amplitudes ranging from less than 1 cm over most of the deep ocean to several centimeters over shallow, boundary regions. Variability generally increases toward the western sides of the basins, and is also larger in some Southern Ocean regions. An oscillation between subtropical and higher latitudes in the North Pacific is clear. Comparison with barotropic simulations indicates that, on basin scales, seasonal Pb variability is related to barotropic dynamics and the seasonal cycle in Ekman pumping, and results from a small, net residual in mass divergence from the balance between Ekman and Sverdrup flows.

Ponte, Rui M.↗

Responses of the Tropical Pacific to Wind Forcing as Observed by Spaceborne Sensors and Simulated by an Ocean General Circulation Model

In this study, satellite observations, in situ measurements, and model simulations are combined to assess the oceanic response to surface wind forcing in the equatorial Pacific. The surface wind fields derived from observations by the spaceborne special sensor microwave imager (SSM/I) and from the operational products of the European Centre for Medium-Range Weather Forecasts (ECMWF) are compared. When SSM/I winds are used to force a primitive-equation ocean general circulation model (OGCM), they produce 3 C more surface cooling than ECMWF winds for the eastern equatorial Pacific during the cool phase of an El Nino-Southern Oscillation event. The stronger cooling by SSM/I winds is in good agreement with measurements at the moored buoys and observations by the advanced very high resolution radiometer, indicating that SSM/I winds are superior to ECMWF winds in forcing the tropical ocean. In comparison with measurements from buoys, tide gauges, and the Geosat altimeter, the OGCM simulates the temporal variations of temperature, steric, and sea level changes with reasonable realism when forced with the satellite winds. There are discrepancies between model simulations and observations that are common to both wind forcing fields, one of which is the simulation of zonal currents; they could be attributed to model deficiencies. By examining model simulations under two winds, vertical heat advection and uplifting of the thermocline are found to be the dominant factors in the anomalous cooling of the ocean mixed layer.

Liu, W. Timothy↗

Ocean Thermal and Color Evolution During the 1997/1998 ENSO Event

A reduced gravity primitive equation modeling and assimilation system is used to study the evolution of the tropical Pacific during the 1997/1998 ENSO cycle. The modeling/assimilation scheme ingests satellite altimeter data and TAO temperature profiles and uses SSM/I satellite derived winds as surface boundary forcing. The four-dimensional structure of the upper ocean circulation structure will be compared against available in situ observations across the Pacific basin. In particular, variability near the Galapagos Islands will be highlighted during the spring of 1998 when phytoplankton concentrations were observed to increase a hundred-fold over a two week period.

Rienecker, Michele↗

An Extended Kalman Filter to Assimilate Altimetric Data into a Non-Linear Model of the Tropical Pacific

A new implementation of the extended Kaman filter is developed for the purpose of assimilating altimetric observations into a primitive equation model of the tropical Pacific. Its specificity consists in defining the errors into a reduced basis that evolves in time with the model dynamic. Validation by twin experiments is conducted and the method is shown to be efficient in quasi real conditions. Data from the first 2 years of the Topex/Poseidon mission are assimilated into the Gent & Cane [1989] model. Assimilation results are evaluated against independent in situ data, namely TAO mooring observations.

Gourdeau, L.↗

Can TOPEX/Poseidon Data be Used for Monitoring the Indonesian Throughflow?

The effects of the Indonesian throughflow(ITF) are studied in a reduced gravity, primitive equation, sigma coordinate model. The surface heat fluxes are provided by coupling the ocean GCAF to an advective atmospheric mixed layer model. There is an ENSO related signal in the ITF but the correlation with the Southern Oscillation Index (SOI) is only -0.31. When the winds over the Indian Ocean are held to climatology, this correlation jumps to -0.65 indicating that the non-EASO signal in the ITF is caused by the downstream winds. On interannual time-scales the ITF can be explained in terms of sea level differences between the western Pacific and eastern Indian Oceans when appropriate representative locations are chosen as demonstrated in both model and TOPEX data. It is shown that the main climatological effect of the ITF is to warm the Indian Ocean and cool the Pacific. Mile the main ENSO indices, NIN03 and NIN04, are almost identical with and without the ITF, the total SSTs show significant ENSO dependence.

Murtugudde, Ragu↗

Modeling of Dense Water Production and Salt Transport from Alaskan Coastal Polynyas

The main significance of this paper is that a realistic, three-dimensional, high-resolution primitive equation model has been developed to study the effects of dense water formation in Arctic coastal polynyas. The model includes realistic ambient stratification, realistic bottom topography, and is forced by time-variant surface heat flux, surface salt flux, and time-dependent coastal flow. The salt and heat fluxes, and the surface ice drift, are derived from satellite observations (SSM/I and NSCAT sensors). The model is used to study the stratification, salt transport, and circulation in the vicinity of Barrow Canyon during the 1996/97 winter season. The coastal flow (Alaska coastal current), which is an extension of the Bering Sea throughflow, is formulated in the model using the wind-transport regression. The results show that for the 1996/97 winter the northeastward coastal current exports 13% to 26% of the salt produced by coastal polynyas upstream of Barrow Canyon in 20 to 30 days. The salt export occurs more rapidly during less persistent polynyas. The inclusion of ice-water stress in the model makes the coastal current slightly weaker and much wider due to the combined effects of surface drag and offshore Ekman transport.

Signorini, Sergio R.↗

Data Assimilation and Transport Modeling in Terrestrial and Planetary Atmospheres

Data assimilation is a blanket term used to describe a number of techniques for retrieving important physical parameters from observational data, subject to constraints imposed by prior knowledge (such as, in the case of meteorology, the primitive equations that govern atmospheric motion). Since these newly developed methods make efficient use of computational resources, they are of great importance in the interpretation of the voluminous datasets that are now produced by satellite missions. As proposed, these techniques have been applied to the study of the Martian and terrestrial atmospheres based on available satellite observations. In addition, a sophisticated hydrodynamic model (non-hydrostatic, and therefore applicable to the study of the interiors of the giant planets) has also been developed and successfully applied to the study of tidally induced motions in Jupiter.

Houben, Howard C.↗

Effects of Convective Asymmetries on Hurricane Intensity: A Numerical Study

The influence of the uniform large-scale flow, beta effect, and vertical shear of the environmental flow on hurricane intensity is investigated in the context of the induced convective or potential vorticity asymmetries with a hydrostatic primitive equation hurricane model. In agreement with the previous studies, imposing of one of these environmental effects can substantially weaken the simulated tropical cyclones. In response t o the environmental influence, significant asymmetries develop with a structure similar to the spiral bands in real hurricanes, which are dominated by wavenumber-one components. The tendencies of the mean radial, azimuthal winds and temperature associated with the environment-induced convective asymmetries are evaluated respectively. The resulting asymmetries can effectively reduce hurricane intensity by directly producing the negative tendency of the mean tangential wind in the vicinity of the radius of maximum wind, and by weakening the mean radial circulation. The reduction effects are closely associated with the spiral structure of the induced asymmetries. The time lag observed between the imposition of the environmental influence and the resulting rise in the minimum central pressure is the time required for developing the spiral structure. This study also confirms the axisymmetrization process associated with the induced wavenumber-one components of potential vorticity asymmetries, but it exists only within the radius of maximum wind.

Wu, Liguang↗

Optimal Observations for Variational Data Assimilation

An important aspect of Ocean state estimation is the design of an observing system that allows the efficient study of climate aspects in the ocean. A solution of the design problem is presented here in terms of optimal observations that emerge as nondimensionalized singular vectors of the modified data resolution matrix. The actual computation is feasible only for scalar quantities in the limit of large observational errors. In the framework of a lo resolution North Atlantic primitive equation model it is demonstrated that such optimal observations when applied to determining the strength of the volume and heat transport across the Greenland-Scotland ridge, perform significantly better than traditional section data. On seasonal to inter-annual time-scales optimal observations are located primarily along the continental shelf and information about heat-transport, wind stress and stratification is being communicated via boundary waves and advective processes. On time-scales of about a month, sea surface height observations appear to be more efficient in reconstructing the cross-ridge heat transport than hydrographic observations. Optimal observations also provide a tool for understanding how the ocean state is effected by anomalies of integral quantities such as meridional heat transport.

Koehl, Armin↗