Search NASA⌕ Search

SEARCH · Search NASA

Results for “Geophysical fluid dynamics”

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 37 records · Page 2

Vertical structure of the wind field during the Special Observing Period I of the Global Weather Experiment

The vertical structure of the global atmosphere is analyzed for selected periods of the Special Observing Period I (SOP-I) for the Global Weather Experiment (GWE). The analysis consists of projection of the stream-function and velocity potential at 200 and 850 mb on spherical harmonics and of the wind and height fields on the normal modes of a linearized form of the primitive equations for a basic state at rest. The kinematic vertical structure is discussed in terms of correlation coefficients of the 200 mb and 850 mb winds and analysis of the internal and external normal modes of the primitive equations. The reliability of the results is checked by applying the same analysis methods to data sets obtained from three different institutions: Geophysical Fluid Dynamics Laboratory (GFDL), European Center for Medium Range Weather Forecasting (ECMWF), and Goddard Laboratory for the Atmospheres (GLA). It is found that, on a global basis, vertically reversing circulations are as important as the equivalent barotropic structures. For the verticaly reversing components, the gravity and mixed Rossby-gravity modes have contributions of the same order of magnitude as those of the Rossby modes in tropical latitudes.

Paegle, J. N.↗

Energetics analysis of the observed and simulated general circulation using three-dimensional normal mode expansions

The energetics characteristics of the observed and simulated general circulation are analyzed using three-dimensional normal mode expansions. The data sets involved are the Goddard Laboratory for Atmospheres (GLA) analysis and simulation data and the Geophysical Fluid Dynamics Laboratory (GFDL) analysis data. The spectral energy properties of the Rossby and gravity modes and energy transformations are presented. Significant influences of model characteristics and the assimilation techniques are observed in the barotropic energy spectrum, particularly for the gravity mode. Energy transformations of the zonal mean field in the GLA analysis and simulation are similar, but distinctly different from that in the GFDL analysis. However, overall, the energy generation in the baroclinic mode is largely balanced by the sink in the barotropic mode. The present study may demonstrate utilities of the three-dimensional normal mode energetics in the analysis of the general circulation.

Tanaka, Hiroshi↗

Energetics examination of winter blocking simulations in the Northern Hemisphere

Four numerical simulations of the global atmosphere for January 1979 are analyzed to study the formation of blocking in terms of Northern Hemisphere energetics. The Goddard Laboratory for Atmospheres (GLA) 4 x 5 deg latitude-longitude grid general circulation model (GCM) and 2 x 2.5 deg grid GCM are employed with the GLA and Geophysical Fluid Dynamics Laboratory initial datasets. Among four simulations, the simulations by the high resolution GCM produce realistically strong blockings with compatible spectral energetics as in the observed blockings episodes. The latitude-height cross sections of the energy variables of wavenumber 1 is presented to describe the dipole structure of blockings. Blocking development is also examined in time series of barotropic and baroclinic components of energy and associated conversions.

Kung, Ernest C.↗

A Pacific Ocean general circulation model for satellite data assimilation

A tropical Pacific Ocean General Circulation Model (OGCM) to be used in satellite data assimilation studies is described. The transfer of the OGCM from a CYBER-205 at NOAA's Geophysical Fluid Dynamics Laboratory to a CRAY-2 at NASA's Ames Research Center is documented. Two 3-year model integrations from identical initial conditions but performed on those two computers are compared. The model simulations are very similar to each other, as expected, but the simulations performed with the higher-precision CRAY-2 is smoother than that with the lower-precision CYBER-205. The CYBER-205 and CRAY-2 use 32 and 64-bit mantissa arithmetic, respectively. The major features of the oceanic circulation in the tropical Pacific, namely the North Equatorial Current, the North Equatorial Countercurrent, the South Equatorial Current, and the Equatorial Undercurrent, are realistically produced and their seasonal cycles are described. The OGCM provides a powerful tool for study of tropical oceans and for the assimilation of satellite altimetry data.

Chao, Y.↗

Studies of the intermediate and deep circulation in the western equatorial Atlantic

This proposal concerns the preparation and design of an experiment, the objective of which is to improve our knowledge of the intermediate and deep circulation in the western equatorial Atlantic Ocean. We shall focus on the description of the western boundary currents, of their crossing with the equator, on the estimation of their mass and heat fluxes, and their seasonal and interannual variations. We will use satellite altimetric data, tomographic measurements, and in situ observations (current measurements, hydrology, and floaters). We propose a feasibility study and the definition of a strategy based on a high-resolution Geophysical Fluid Dynamics Laboratory (GFDL) numerical model to define which in situ measurements are necessary to optimally complete the altimetric observations.

Desaubies, Yves↗

The impact of tropical wind data on the analysis and forcasts of the GLA GCM for the global weather experiment

It is well-known that divergent wind estimates are much more dependent upon the analysis system than are estimates of the rotational wind. This conclusion is supported in recent analyses of FGGE SOP1 data produced by the Goddard Laboratory for Atmospheres (GLA), the Geophysical Fluid Dynamics Laboratory (GFDL) and the European Center for Medium Range Weather Forecasting (ECMWF). These analyses differ in the forecast models that are used for the four-dimensional assimilation, in the data rejection criteria, and, to a certain extent, in the data density. Because the final divergent wind is a product of both model constraints and observation, it is relevant to inquire how much of each goes into the final product. We presently investigate this question through a systematic analysis of tropical data that are sampled at different densities by the GLA GCM.

Paegle, Jan↗

Investigation of finite element: ABC methods for electromagnetic field simulation

The mechanics of wave propagation in the presence of obstacles is of great interest in many branches of engineering and applied mathematics like electromagnetics, fluid dynamics, geophysics, seismology, etc. Such problems can be broadly classified into two categories: the bounded domain or the closed problem and the unbounded domain or the open problem. Analytical techniques have been derived for the simpler problems; however, the need to model complicated geometrical features, complex material coatings and fillings, and to adapt the model to changing design parameters have inevitably tilted the balance in favor of numerical techniques. The modeling of closed problems presents difficulties primarily in proper meshing of the interior region. However, problems in unbounded domains pose a unique challenge to computation, since the exterior region is inappropriate for direct implementation of numerical techniques. A large number of solutions have been proposed but only a few have stood the test of time and experiment. The goal of this thesis is to develop an efficient and reliable partial differential equation technique to model large three dimensional scattering problems in electromagnetics.

Chatterjee, A.↗

Precipitation recycling in the Amazon basin

Precipitation recycling is the contribution of evaporation within a region to precipitation in that same region. The recycling rate is a diagnostic measure of the potential for interactions between land surface hydrology and regional climate. In this paper we present a model for describing the seasonal and spatial variability of the recycling process. The precipitation recycling ratio, rho, is the basic variable in describing the recycling process. Rho is the fraction of precipitation at a certain location and time which is contributed by evaporation within the region under study. The recycling model is applied in studyiing the hydrologic cycle in the Amazon basin. It is estimated that about 25% of all the rain that falls in the Amazon basin is contributed by evaporation within the basin. This estimate is based on analysis of a data set supplied by the European Centre for Medium-range Weather Forecasts (ECMWF). The same analysis is repeated using a different data set from the Geophysical Fluid Dynamics Laboratory (GFDL). Based on this data set, the recycling ratio is estimated to be 35%. The seasonal variability of the recycling ratio is small compared with the yearly average. The new estimates of the recycling ratio are compared with results of previous studies, and the differences are explained.

Eltahir, E. A. B.↗

A simple hydrologically based model of land surface water and energy fluxes for general circulation models

A generalization of the single soil layer variable infiltration capacity (VIC) land surface hydrological model previously implemented in the Geophysical Fluid Dynamics Laboratory (GFDL) general circulation model (GCM) is described. The new model is comprised of a two-layer characterization of the soil column, and uses an aerodynamic representation of the latent and sensible heat fluxes at the land surface. The infiltration algorithm for the upper layer is essentially the same as for the single layer VIC model, while the lower layer drainage formulation is of the form previously implemented in the Max-Planck-Institut GCM. The model partitions the area of interest (e.g., grid cell) into multiple land surface cover types; for each land cover type the fraction of roots in the upper and lower zone is specified. Evapotranspiration consists of three components: canopy evaporation, evaporation from bare soils, and transpiration, which is represented using a canopy and architectural resistance formulation. Once the latent heat flux has been computed, the surface energy balance is iterated to solve for the land surface temperature at each time step. The model was tested using long-term hydrologic and climatological data for Kings Creek, Kansas to estimate and validate the hydrological parameters, and surface flux data from three First International Satellite Land Surface Climatology Project Field Experiment (FIFE) intensive field campaigns in the summer-fall of 1987 to validate the surface energy fluxes.

Liang, XU↗

Decadal variability of the Tropical Atlantic Ocean Surface Temperature in shipboard measurements and in a Global Ocean-Atmosphere model

Sea surface temperature (SST) variability was investigated in a 200-yr integration of a global model of the coupled oceanic and atmospheric general circulations developed at the Geophysical Fluid Dynamics Laboratory (GFDL). The second 100 yr of SST in the coupled model's tropical Atlantic region were analyzed with a variety of techniques. Analyses of SST time series, averaged over approximately the same subregions as the Global Ocean Surface Temperature Atlas (GOSTA) time series, showed that the GFDL SST anomalies also undergo pronounced quasi-oscillatory decadal and multidecadal variability but at somewhat shorter timescales than the GOSTA SST anomalies. Further analyses of the horizontal structures of the decadal timescale variability in the GFDL coupled model showed the existence of two types of variability in general agreement with results of the GOSTA SST time series analyses. One type, characterized by timescales between 8 and 11 yr, has high spatial coherence within each hemisphere but not between the two hemispheres of the tropical Atlantic. A second type, characterized by timescales between 12 and 20 yr, has high spatial coherence between the two hemispheres. The second type of variability is considerably weaker than the first. As in the GOSTA time series, the multidecadal variability in the GFDL SST time series has approximately opposite phases between the tropical North and South Atlantic Oceans. Empirical orthogonal function analyses of the tropical Atlantic SST anomalies revealed a north-south bipolar pattern as the dominant pattern of decadal variability. It is suggested that the bipolar pattern can be interpreted as decadal variability of the interhemispheric gradient of SST anomalies. The decadal and multidecadal timescale variability of the tropical Atlantic SST, both in the actual and in the GFDL model, stands out significantly above the background 'red noise' and is coherent within each of the time series, suggesting that specific sets of processes may be responsible for the choice of the decadal and multidecadal timescales. Finally, it must be emphasized that the GFDL coupled ocean-atmosphere model generates the decadal and multidecadal timescale variability without any externally applied force, solar or lunar, at those timescales.

Mehta, Vikram M.↗

Comparison of tropical pacific temperature and current simulations with two vertical mixing schemes embedded in an ocean general circulation model and reference to observations

The Pacanowski-Philander (PP) and Mellor-Yamada (MY) parameterization models of vertical mixing by turbulent processes were embedded in the Geophysical Fluid Dynamics Laboratory high-resolution ocean general circulation model of the tropical Pacific Ocean. All other facets of the numerical simulations were the same. Simulations were made for the 1987-1988 period. At the equator the MY simulation produced near-surface temperatures more uniform with depth, a deeper thermocline, a deeper core speed of the Equatorial Undercurrent, and a South Equatorial Current with greater vertical thickness compared with that computed with the PP method. Along 140 deg W, between 5 deg N and 10 deg N, both simulations were the same. Moored buoy current and temperature observations had been recorded by the Pacific Marine Environmental Laboratory at three sites (165 deg E, 140 deg W, 110 deg W) along the equator and at three sites (5 deg N, 7 deg N, 9 deg N) along 140 deg W. Simulated temperatures were lower than those observed in the near-surface layer and higher than those observed in the thermocline. Temperature simulations were in better agreement with observations compared to current simulations. At the equator, PP current and temperature simulations were more representative of the observations than MY simulations.

Halpern, David↗

Wavenumber Vacillation in Weakly-Stratified Baroclinic Flows

The wavenumber selection in an unstable baroclinic flow has long been a subject of investigation in geophysical fluid dynamics. For a weakly unstable flow, the linearly most unstable wave is always the final realized wave. As the flow becomes more unstable, the dominant wave gradually shifts from the most unstable wave to a longer wave. However, when the supercriticality is reduced from the state with a dominant longwave, it may remain in that state even beyond the threshold where it first developed. Hysteresis in a quasigeostrophic system also displays a complex behavior: it occurs not only between the states of different dominant wavenumbers, but also between the states of identical dominant wavenumber but of different dynamic characteristics. In a recent numerical experiment using a high-resolution spectral Eady model with asymmetric Ekman dissipation, it is found that in a weakly stratified flow with strong meridional forcing the flow does not show a preference for any single wave. Instead, the dominant wave periodically varies among the linearly most unstable wave and the longer waves. The mechanism of wavenumber vacillation will be examined in terms of meridional heat flux, the wave PV gradient, and wave saturation. The importance of wave-wave and wave-mean flow interactions will also be investigated.

Chou, Shih-Hung↗

Turbulent Diffusivities for Momentum, Heat, Salt and Passive Scalars

A program began ten years ago to build a turbulence model to describe high Reynolds numbers flows. Specifically, the aim was to devise a model that would satisfy two basic conditions: reproduce available turbulence data (laboratory, numerical simulations, etc.) concerning shear driven flows, buoyancy driven flows, 2D turbulence, freely decaying turbulence etc., and yet be manageable so as to be used, for example, in General Circulation Models (GCM'S). The model was presented in a series of papers that have appeared in Physics of Fluids since 1996. A total of about 80 turbulence statistics were reproduced. The model has no adjustable parameters. The next step was to apply the model to construct the vertical diffusities K for moment K(sub m), temperature K(sub h), salt K(sub s). and passive scalars K(sub c). First, we took K(sub s) = K(sub h) and tested the model using the GFDL ocean model. The results for the profiles of T and S vs. depth are indistinguishable from those derived using the latest model, the KPP model by the NCAR group. Presently, we are running the same Geophysical Fluid Dynamics Laboratory (GFDL) code relaxing the assumptive K(sub s) = K(sub h). Indeed, the turbulence model yields a salt diffusivity that depends on Ri and R rho (= Beta delta S/delta z/alpha delta T/delta z) in such a way that K(sub s) may be quite different from K(sub h). Salt fingers and double diffusivity laboratory data are reproduced. Results from the ocean model will be available shortly. Finally, we are trying to derive the horizontal diffusivities with the goal of providing a physically acceptable representation of mesoscale eddies. The recently suggested GMW parameterization has improved several O-GCM results and the goal here is to try to derive/justify it from a turbulence model and/or propose improvements/modifications. Theoretical work is in progress.

Canuto, V. M.↗

Turbulent Diffusivities for Momentum, Heat, Salt and Passive Scalars

A program began ten years ago to build a turbulence model to describe high Reynolds numbers flows. Specifically, the aim was to devise a model that would satisfy two basic conditions: reproduce available turbulence data (laboratory, numerical simulations, etc.) concerning shear driven flows, buoyancy driven flows, 2D turbulence, freely decaying turbulence etc., and yet be manageable so as to be used, for example, in GCM'S. The model was presented in a series of papers that have appeared in Physics of FLuids since 1996. A total of about 80 turbulence statistics were reproduced. The model has no adjustable parameters. The next step was to apply the model to construct the vertical diffusivities K for moment K(sub m), temperature K(sub h), salt K(sub g) and passive scalars K(sub c). First, we took K(sub g) = K(sub h) and tested the model using the Geophysical Fluid Dynamics Laboratory (GFDL) ocean model. The results for the profiles of T and S vs. depth are indistinguishable from those derided using the latest model, the KPP model by the NCAR group. Presently, we are running the same GFDL code relaxing the assumptive K(sub g) = K(sub h). Indeed, the turbulence model yields a salt diffusivity that depends on Ri and R Rho (= Beta Sigma S/Sigma z/Alpha Sigma T/Sigma z) in such a way that K(sub g) may be quite different from K(sub h). Salt fingers and double diffusivity laboratory data are reproduced. Results from the ocean model will be available shortly. Finally, we are trying to derive the horizontal diffusivities with the goal of providing a physically acceptable representation of mesoscale eddies. The recently suggested GMW parameterization has improved several O-GCM results and the goal here is to try to derive/justify it from a turbulence model and/or propose improvements/modifications. Theoretical work is in progress.

Canuto, Vittorio↗

Sensitivity of Age-of-Air Calculations to the Choice of Advection Scheme

The age of air has recently emerged as a diagnostic of atmospheric transport unaffected by chemical parameterizations, and the features in the age distributions computed in models have been interpreted in terms of the models' large-scale circulation field. This study shows, however, that in addition to the simulated large-scale circulation, three-dimensional age calculations can also be affected by the choice of advection scheme employed in solving the tracer continuity equation, Specifically, using the 3.0deg latitude X 3.6deg longitude and 40 vertical level version of the Geophysical Fluid Dynamics Laboratory SKYHI GCM and six online transport schemes ranging from Eulerian through semi-Lagrangian to fully Lagrangian, it will be demonstrated that the oldest ages are obtained using the nondiffusive centered-difference schemes while the youngest ages are computed with a semi-Lagrangian transport (SLT) scheme. The centered- difference schemes are capable of producing ages older than 10 years in the mesosphere, thus eliminating the "young bias" found in previous age-of-air calculations. At this stage, only limited intuitive explanations can be advanced for this sensitivity of age-of-air calculations to the choice of advection scheme, In particular, age distributions computed online with the National Center for Atmospheric Research Community Climate Model (MACCM3) using different varieties of the SLT scheme are substantially older than the SKYHI SLT distribution. The different varieties, including a noninterpolating-in-the-vertical version (which is essentially centered-difference in the vertical), also produce a narrower range of age distributions than the suite of advection schemes employed in the SKYHI model. While additional MACCM3 experiments with a wider range of schemes would be necessary to provide more definitive insights, the older and less variable MACCM3 age distributions can plausibly be interpreted as being due to the semi-implicit semi-Lagrangian dynamics employed in the MACCM3. This type of dynamical core (employed with a 60-min time step) is likely to reduce SLT's interpolation errors that are compounded by the short-term variability characteristic of the explicit centered-difference dynamics employed in the SKYHI model (time step of 3 min). In the extreme case of a very slowly varying circulation, the choice of advection scheme has no effect on two-dimensional (latitude-height) age-of-air calculations, owing to the smooth nature of the transport circulation in 2D models. These results suggest that nondiffusive schemes may be the preferred choice for multiyear simulations of tracers not overly sensitive to the requirement of monotonicity (this category includes many greenhouse gases). At the same time, age-of-air calculations offer a simple quantitative diagnostic of a scheme's long-term diffusive properties and may help in the evaluation of dynamical cores in multiyear integrations. On the other hand, the sensitivity of the computed ages to the model numerics calls for caution in using age of air as a diagnostic of a GCM's large-scale circulation field.

Eluszkiewicz, Janusz↗

Mixing in the Extratropical Stratosphere: Model-measurements Comparisons using MLM Diagnostics

We evaluate transport processes in the extratropical lower stratosphere for both models and measurements with the help of equivalent length diagnostic from the modified Lagrangian-mean (MLM) analysis. This diagnostic is used to compare measurements of long-lived tracers made by the Cryogenic Limb Array Etalon Spectrometer (CLAES) on the Upper Atmosphere Research Satellite (UARS) with simulated tracers. Simulations are produced in Chemical and Transport Models (CTMs), in which meteorological fields are taken from the Goddard Earth Observing System Data Assimilation System (GEOS DAS), the Middle Atmosphere Community Climate Model (MACCM2), and the Geophysical Fluid Dynamics Laboratory (GFDL) "SKYHI" model, respectively. Time series of isentropic equivalent length show that these models are able to capture major mixing and transport properties observed by CLAES, such as the formation and destruction of polar barriers, the presence of surf zones in both hemispheres. Differences between each model simulation and the observation are examined in light of model performance. Among these differences, only the simulation driven by GEOS DAS shows one case of the "top-down" destruction of the Antarctic polar vortex, as observed in the CLAES data. Additional experiments of isentropic advection of artificial tracer by GEOS DAS winds suggest that diabatic movement might have considerable contribution to the equivalent length field in the 3D CTM diagnostics.

Ma, Jun↗

The Impact of Land Initialization and Assimilation on Climate Predictability and Prediction

Analysis will be presented which explores the impact of land conditions on monthly to seasonal climate simulations in a variety of atmospheric general circulation models (AGCMs). In one set of experiments, the Geophysical Fluid Dynamics Laboratory (GDFL) AGCM is used to explore the nature of soil-moisture predictability and associated climate predictability as an initial value problem. For another set of experiments, the Center for Ocean Land Atmosphere (COLA) and the Goddard Earth Observing System 2 (GEOS-2) AGCMs are used to investigate the impact of realistic snow initialization and assimilation in retrospective climate forecasts for the northern hemisphere spring (March-June).

Schlosser, C. Adam↗

Assessing Impacts of Global Warming on Tropical Cyclone Tracks

A new approach is proposed to assess the possible impacts of the global climate change on tropical cyclone (TC) tracks in the western North Pacific (WNP) basin. The idea is based on the premise that the future change of TC track characteristics is primarily determined by changes in large-scale environmental steering flows. It is demonstrated that the main characteristics of the current climatology of TC tracks can be derived from the climatological mean velocity field of TC motion by using a trajectory model. The climatological mean velocity of TC motion, which is composed of the large-scale steering and beta drift, is determined on each grid of the basin. The mean beta drift is estimated from the best track data, and the mean large-scale steering flow is computed from the NCEP/NCAR reanalysis for the current climate state. The derived mean beta drift agrees well with the results of previous observational and numerical studies in terms of its direction and magnitude. The outputs of experiments A2 and B2 of the Geophysical Fluid Dynamics Laboratory (GFDL) R30 climate model suggest that the subtropical high will be persistently weak over the western part of the WNP or shift eastward during July-September in response to the future climate change. By assuming that the mean beta drift in the future climate state is unchanged, the change in the general circulation by 2059 will decrease the TC activities in the WNP, but favor a northward shift of typical TC tracks. As a result, the storm activities in the South China Sea will decrease by about 12%, while the Japan region will experience an increase of TCs by 12-15%. During the period of 2000-2029, the tropical storms that affect the China region will increase by 5-6%, but return to the current level during 2030-2059. It is also suggested that, during the period of 2030-2059 tropical storms will more frequently affect Japan and the middle latitude region of China given that the formation locations remain the same as in the current climate state.

Wu, Li-Guang↗