Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “GCM”

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 433 records · Page 24
All sourcesNASA NTRSDOE OSTI
NASA NTRS1994

A class of the van Leer-type transport schemes and its application to the moisture transport in a general circulation model

A generalized form of the second-order van Leer transport scheme is derived. Several constraints to the implied subgrid linear distribution are discussed. A very simple positive-definite scheme can be derived directly from the generalized form. A monotonic version of the scheme is applied to the Goddard Laboratory for Atmospheres (GLA) general circulation model (GCM) for the moisture transport calculations, replacing the original fourth-order center-differencing scheme. Comparisons with the original scheme are made in idealized tests as well as in a summer climate simulation using the full GLA GCM. A distinct advantage of the monotonic transport scheme is its ability to transport sharp gradients without producing spurious oscillations and unphysical negative mixing ratio. Within the context of low-resolution climate simulations, the aforementioned characteristics are demonstrated to be very beneficial in regions where cumulus convection is active. The model-produced precipitation pattern using the new transport scheme is more coherently organized both in time and in space, and correlates better with observations. The side effect of the filling algorithm used in conjunction with the original scheme is also discussed, in the context of idealized tests. The major weakness of the proposed transport scheme with a local monotonic constraint is its substantial implicit diffusion at low resolution. Alternative constraints are discussed to counter this problem.

Lin, Shian-Jiann↗
NASA NTRS1994

Modeling global lightning distributions in a general circulation model

A general circulation model (GCM) is used to model global lightning distributions and frequencies. Both total and cloud-to-ground lightning frequencies are modeled using parameterizations that relate the depth of convective clouds to lightning frequencies. The model's simulations of lightning distributions in time and space show good agreement with available observations. The model's annual mean climatology shows a global lightning frequency of 77 flashes per second, with cloud-to-ground lightning making up 25% of the total. The maximum lightning activity in the GCM occurs during the Northern Hemisphere summer, with approximately 91% of all lightning occurring over continental and coastal regions.

Price, Colin↗
NASA NTRS1993

Transport and residence times of tropospheric aerosols inferred from a global three-dimensional simulation of Pb-210

A global three-dimensional model is used to investigate the transport and tropospheric residence time of Pb-210, an aerosol tracer produced in the atmosphere by radioactive decay of Rn-222 emitted from soils. The model uses meteorological input with 4 deg x 5 deg horizontal resolution and 4-hour temporal resolution from the Goddard Institute for Space Studies general circulation model (GCM). It computes aerosol scavenging by convective precipitation as part of the wet convective mass transport operator in order to capture the coupling between vertical transport and rainout. Scavenging in convective precipitation accounts for 74% of the global Pb-210 sink in the model; scavenging in large-scale precipitation accounts for 12%, and scavenging in dry deposition accounts for 14%. The model captures 63% of the variance of yearly mean Pb-210 concentrations measured at 85 sites around the world with negligible mean bias, lending support to the computation of aerosol scavenging. There are, however, a number of regional and seasonal discrepancies that reflect in part anomalies in GCM precipitation. Computed residence times with respect to deposition for Pb-210 aerosol in the tropospheric column are about 5 days at southern midlatitudes and 10-15 days in the tropics; values at northern midlatitudes vary from about 5 days in winter to 10 days in summer. The residence time of Pb-210 produced in the lowest 0.5 km of atmosphere is on average four times shorter than that of Pb-210 produced in the upper atmosphere. Both model and observations indicate a weaker decrease of Pb-210 concentrations between the continental mixed layer and the free troposphere than is observed for total aerosol concentrations; an explanation is that Rn-222 is transported to high altitudes in wet convective updrafts, while aerosols and soluble precursors of aerosols are scavenged by precipitation in the updrafts. Thus Pb-210 is not simply a tracer of aerosols produced in the continental boundary layer, but also of aerosols derived from insoluble precursors emitted from the surface of continents. One may draw an analogy between Pb-210 and nitrate, whose precursor NO(sub x) is sparingly soluble, and explain in this manner the strong correlation observed between nitrate and Pb-210 concentrations over the oceans.

Balkanski, Yves J.↗
NASA NTRS1994

The interpretation of remotely sensed cloud properties from a model paramterization perspective

A study has been made of the relationship between mean cloud radiative properties and cloud fraction in stratocumulus cloud systems. The analysis is of several Land Resources Satellite System (LANDSAT) images and three hourly International Satellite Cloud Climatology Project (ISCCP) C-1 data during daylight hours for two grid boxes covering an area typical of a general circulation model (GCM) grid increment. Cloud properties were inferred from the LANDSAT images using two thresholds and several pixel resolutions ranging from roughly 0.0625 km to 8 km. At the finest resolution, the analysis shows that mean cloud optical depth (or liquid water path) increases somewhat with increasing cloud fraction up to 20% cloud coverage. More striking, however, is the lack of correlation between the two quantities for cloud fractions between roughly 0.2 and 0.8. When the scene is essentially overcast, the mean cloud optical tends to be higher. Coarse resolution LANDSAT analysis and the ISCCP 8-km data show lack of correlation between mean cloud optical depth and cloud fraction for coverage less than about 90%. This study shows that there is perhaps a local mean liquid water path (LWP) associated with partly cloudy areas of stratocumulus clouds. A method has been suggested to use this property to construct the cloud fraction paramterization in a GCM when the model computes a grid-box-mean LWP.

HARSHVARDHAN↗
NASA NTRS1994

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↗
NASA NTRS1995

Geology and landscape evolution of the Hellas region of Mars

Hellas basin on Mars has been the site of volcanism, tectonism, and modification by fluvial, mass-wasting, and eolian processes over its more than 4-b.y. existence. Our detailed geologic mapping and related studies have resulted in the following new interpretations. The asymmetric distribution of highland massifs and other structures that define the uplifted basin rim suggest a formation of the basin by the impact of a low-angle bolide having a trajectory heading S60E. During the Late Noachian, the basin was infilled, perhaps by lava flows, that were sufficiently thick (greater than 1 km) to produce wrinkle ridges on the fill material and extensional faulting along the west rim of the basin. At about the same time, deposits buried northern Malea Planum, which are interpreted to be pyroclastic flows from Amphitrites and Peneus Paterae on the basis of their degraded morphology, topology, and the application of a previous model for pyroclastic volcanism on Mars. Peneus forms a distinctive caldera structure that indicates eruption of massive volumes of magma, whereas Amphitrites is a less distinct circular feature surrounded by a broad, low, dissected shield that suggests generally smaller volume eruptions. During the Early Hesperian, an approximately 1-to 2km-thick sequence of primarily fined-grained, eolian material was deposited on the floor of Hellas basin. Subsequently, the deposit was deeply eroded, except where armored by crater ejecta, and it retreated as much as 200-300 km along its western margin, leaving behind pedestal craters and knobby outliers of the deposit. Local debris flows within the deposit attest to concentrations of groundwater, perhaps in part brought in by outflow floods along the east rim of the basin. These floods may have deposited approximately 100-200m of sediment, subduing wrinkle ridges in the eastern part of the basin floor. During the Late Hesperian and Amazonian, eolian mantles were emplaced on the basin rim and floor and surrounding highlands. Their subsequent erosion resulted in pitted and etched plains and crater fill, irregular mesas, and pedestal craters. Local evidence occurs for the possible former presence of ground ice or ice sheets approximately 100 km across; however, we disagree with a hypothesis that suggest that the entire south rim and much of the floor of Hellas have been glaciated. Orientations of dune fields and yardangs in lower parts of Hellas basin follow directions of the strongest winds predicted by a recently published general circulation model (GCM). Transient frost and dust splotches in the region are, by contrast, related to the GCM prediction for the season in which the images they appear in were taken.

Tanaka, Kenneth L.↗
NASA NTRS1994

The scale and persistence of soil moisture anomalies as simulated in a global model

Short term variability of climate is intimately connected with soil moisture variability. Soil moisture provides the storage and subsequent return to the atmosphere, through evaporation and transpiration, of precipitation anomalies over land. Global Circulation Model (GCM) simulations enable consistent identification of correlations and dynamical connections between the hydrologic variables, many of which are incompletely observed. One way to facilitate understanding with these increasingly intricate models is to perform sensitivity studies in which a boundary condition or process is prescribed. In this study we will report on a sensitivity study in which a GCM with a sophisticated land surface representation is used to investigate soil moisture variability in the model climate. The simulations to be used in this study were made at R15 resolution (approximately 4.5 deg latitude x 7.5 deg longitude) with prescribed sea surface temperatures (SST) in the GENESIS model (Thompson and Pollard, 1994), which is coupled to a Land Surface Transfer model (LSX) at 2 deg x 2 deg resolution (Pollard and Thompson, 1994). All the results represented here were taken from the monthly averages of the model results. The LSX model accounts for the physical effects of vegetation with two layers specified at each grid point. Vegetation attributes such as leaf area indices, fractional cover, leaf albedos, etc., were taken from the global dataset in Dorman and Sellers (1989). A six-layer soil model extends from the surface to 4.25 m depth. SST's were prescribed in two ten year experiments using monthly SST values with the daily value being interpolated from the nearest two months. In the first experiment monthly climatological values were used, and in the second, the Atmospheric Model Intercomparison Project (AMIP) observed SST's for the years 1979 through 1988 were used (Gates, 1992). Thus, the former experiment gives a measure of the intrinsic model variability, to be compared with that of the latter experiment, which includes month-to-month variability due to ocean forcing.

Fitzjarrald, D.↗
NASA NTRS1994

Calculation of clear-sky outgoing longwave radiation using ECMWF gridded fields and ISCCP C1 cloud data

Clear-Sky OLR (outgoing longwave radiation) has a variety of definitions. Clear-sky OLR can be measured from satellites by examining the data to detect clouds, eliminating those observations that are contaminated with clouds, and averaging the remaining data. The disadvantage of this method is that satellites do not provide time continuous coverage, requiring many observations and careful analysis to remove sampling biases. Alternatively, clear-sky OLR can be found continuously using a radiative flux code, combined with analyzed fields of temperature, humidity and clouds, either from assimilated observations or a GCM simulation. The disadvantage is that the cloud field is typically very low resolution, on the order of 200 km, so they lack important small scale cloud detail. Even if fractional clouds data is available, it is not obvious how GCM and ERBE results should be compared. The purpose of this study is to assess clear-sky OLR calculated from these methods using available data. Cloud data from ISCCP provides the fractional cloud coverage data necessary to test these methods against calculations based on observed data. Clear-sky OLR is calculated using the ECMWF/TOGA archive. Monthly averages are made in the manner of the several methods listed above and compared to each other and to ERBE.

Ramsey, Perry G.↗
NASA NTRS1995

Remote Sensing and Holocene Vegetation: History of Global Change

Predictions of the future evolution of the earth's atmospheric chemistry and its impact on global circulation patterns are based on Global Climate Models (GCM's) that integrate the complex interactions of the biosphere, atmosphere and the oceans. Most of the available records of climate and environment are shortterm records (from decades to a few hundred years) with convolved information of real trends and short-term fluctuations. GCM's must be tested beyond the short-term record of climate and environment to insure that predictions are based on trends and therefore are appropriate to support long term policy making. An appropriate timeframe should extend over the Holocene period (the last 10,000 years) when most contemporary climate and environmental processes began. Since its inception in 1916, pollen analysis has successfully reconstructed the paleoecology of the last 10,000 years for many sites around the world, thus providing a powerful time-link between short- and long-term processes in the biosphere. However, pollen analytic results cannot be used in physiological models driven by remotely sensed data. Further, modern ecology and climate data are necessary to calibrate pollen analytical models. These are available for extensive regions in the northern hemisphere, particularly for eastern United States and Canada, and western Europe. In other parts of the world, weather stations are scattered, records extend over a period of only few years, and there are no systematic climate records for large portions of the globe. This is the case of Patagonia in Argentina where a few weather stations are located close to the Atlantic seaports, fewer stations are in towns located near the eastern Andean foothill, and fewer still are scattered on the extensive Patagonian plateau. This problem became evident after completion of the Argentine-German Program of Palynology (PROPAL), a cooperative effort of National University of Mar del Plata (Argentina) and University Bamberg (Germany) to produce a modern pollen database for the Pampa and Patagonia regions.

DAntoni, Hector↗
NASA NTRS1997

Parameterizing Grid-Averaged Longwave Fluxes for Inhomogeneous Marine Boundary Layer Clouds

This paper examines the relative impacts on grid-averaged longwave flux transmittance (emittance) for Marine Boundary Layer (MBL) cloud fields arising from horizontal variability of optical depth tau and cloud sides, First, using fields of Landsat-inferred tau and a Monte Carlo photon transport algorithm, it is demonstrated that mean all-sky transmittances for 3D variable MBL clouds can be computed accurately by the conventional method of linearly weighting clear and cloudy transmittances by their respective sky fractions. Then, the approximations of decoupling cloud and radiative properties and assuming independent columns are shown to be adequate for computation of mean flux transmittance. Since real clouds have nonzero geometric thicknesses, cloud fractions A'(sub c) presented to isotropic beams usually exceed the more familiar vertically projected cloud fractions A(sub c). It is shown, however, that when A(sub c)less than or equal to 0.9, biases for all-sky transmittance stemming from use of A(sub c) as opposed to A'(sub c) are roughly 2-5 times smaller than, and opposite in sign to, biases due to neglect of horizontal variability of tau. By neglecting variable tau, all-sky transmittances are underestimated often by more than 0.1 for A(sub c) near 0.75 and this translates into relative errors that can exceed 40% (corresponding errors for all-sky emittance are about 20% for most values of A(sub c). Thus, priority should be given to development of General Circulation Model (GCM) parameterizations that account for the effects of horizontal variations in unresolved tau, effects of cloud sides are of secondary importance. On this note, an efficient stochastic model for computing grid-averaged cloudy-sky flux transmittances is furnished that assumes that distributions of tau, for regions comparable in size to GCM grid cells, can be described adequately by gamma distribution functions. While the plane-parallel, homogeneous model underestimates cloud transmittance by about an order of magnitude when 3D variable cloud transmittances are less than or equal to 0.2 and by approx. 20% to 100% otherwise, the stochastic model reduces these biases often by more than 80%.

Barker, Howard W.↗
NASA NTRS1994

Two-Layer Variable Infiltration Capacity Land Surface Representation for General Circulation Models

A simple two-layer variable infiltration capacity (VIC-2L) land surface model suitable for incorporation in general circulation models (GCMs) is described. The model consists of a two-layer characterization of the soil within a GCM grid cell, and uses an aerodynamic representation of latent and sensible heat fluxes at the land surface. The effects of GCM spatial subgrid variability of soil moisture and a hydrologically realistic runoff mechanism are represented in the soil layers. The model was tested using long-term hydrologic and climatalogical data for Kings Creek, Kansas to estimate and validate the hydrological parameters. Surface flux data from three First International Satellite Land Surface Climatology Project Field Experiments (FIFE) intensive field compaigns in the summer and fall of 1987 in central Kansas, and from the Anglo-Brazilian Amazonian Climate Observation Study (ABRACOS) in Brazil were used to validate the mode-simulated surface energy fluxes and surface temperature.

GENERAL CIRCULATION MODELS↗
NASA NTRS1998

Incorporating Parallel Computing into the Goddard Earth Observing System Data Assimilation System (GEOS DAS)

Atmospheric data assimilation is a method of combining actual observations with model forecasts to produce a more accurate description of the earth system than the observations or forecast alone can provide. The output of data assimilation, sometimes called the analysis, are regular, gridded datasets of observed and unobserved variables. Analysis plays a key role in numerical weather prediction and is becoming increasingly important for climate research. These applications, and the need for timely validation of scientific enhancements to the data assimilation system pose computational demands that are best met by distributed parallel software. The mission of the NASA Data Assimilation Office (DAO) is to provide datasets for climate research and to support NASA satellite and aircraft missions. The system used to create these datasets is the Goddard Earth Observing System Data Assimilation System (GEOS DAS). The core components of the the GEOS DAS are: the GEOS General Circulation Model (GCM), the Physical-space Statistical Analysis System (PSAS), the Observer, the on-line Quality Control (QC) system, the Coupler (which feeds analysis increments back to the GCM), and an I/O package for processing the large amounts of data the system produces (which will be described in another presentation in this session). The discussion will center on the following issues: the computational complexity for the whole GEOS DAS, assessment of the performance of the individual elements of GEOS DAS, and parallelization strategy for some of the components of the system.

Larson, Jay W.↗
NASA NTRS1995

Geology and landscape evolution of the Hellas region of Mars

Hellas basin on Mars has been the site of volcanism, tectonism, and modification by fluvial, mass-wasting, and eolian processes over its more than 4-b.y. existence. Our detailed geologic mapping and related studies have resulted in the following new interpretations. The asymmetric distribution of highland massifs and other structures that define the uplifted basin rim suggest a formation of the basin by the impact of a low-angle bolide having a trajectory heading S 60 deg E. During the Late Noachian, the basin was infilled, perhaps by lava flows, that were sufficiently thick (>1 km) to produce wrinkle ridges on the fill material and extensional faulting along the west rim of the basin. At about the same time, deposits buried northern Malea Planum, which are interpreted to be pyroclastic flows from Amphitrites and Peneus Paterae on the basis of their degraded morphology, topography, and the application of a previous model for pyroclastic volcanism on Mars. Peneus forms a distinctive caldera structure that indicates eruption of massive volumes of magma, whereas Amphitrites is a less distinct circular feature surrounded by a broad, low, dissected shield that suggests generally smaller volume eruptions. During the Early Hesperian, an approx. 1- to 2-km-thick sequence of primarily fined-grained, eolian material was deposited on the floor of Hellas basin. Subsequently, the deposit was deeply eroded, except where armored by crater ejecta, and it retreated as much as 200-300 km along its western margin, leaving behind pedestal craters and knobby outliers of the deposit. Local debris flows within the deposit attest to concentrations of groundwater, perhaps in part brought in by outflow floods along the east rim of the basin. These floods may have deposited approx. 100-200 m of sediment, subduing wrinkle ridges in the eastern part of the basin floor. During the Late Hesperian and Amazonian, eolian mantles were emplaced on the basin rim and floor and surrounding highlands. Their subsequent erosion resulted in pitted and etched plains and crater fill, irregular mesas, and pedestal craters. Local evidence occurs for the possible former presence of ground ice or ice sheets approx. 100 km across; however, we disagree with a hypothesis that suggests that the entire south rim and much of the floor of Hellas have been glaciated. Orientations of dune fields and yardangs in lower parts of Hellas basin follow directions of the strongest winds predicted by a recently published general circulation model (GCM). Transient frost and dust splotches in the region are, by contrast, related to the GCM prediction for the season in which the images they appear in were taken.

Tanaka, Kenneth L.↗
NASA NTRS1998

MAMI: Modeling of the Magnetosphere-Ionosphere-Atmosphere System

Major emphasis of the investigation was the development of theoretical and numerical models of the aurora and of high latitude ionospheric processes. In particular: (1) the NCAR TIGCM (Thermosphere-Ionosphere Global Circulation Model) was updated to include mid and low latitude electrodynamics (this version is called TIE-GCM); (2) the NCAR TIE-GCM was modified to include a more realistic representation of the aurora; (3) the UAF auroral electron transport model was modified to include local acceleration processes; (4) a local ionospheric and auroral model (AURORA) was developed to allow detailed studies of the aurora; (5) a proton-hydrogen transport code has been developed to model proton aurora; and (6) a theory for the production of suprathermal atoms and ions in the upper atmosphere was developed and applied to studies of atomic nitrogen transport and helium escape on open field lines. These models enable us to devise schemes for tile interpretation and quantitative analysis of data obtained by the POLAR spacecraft. Parameterizations were formulated and are available to the GGS investigators. The UVI and VIS teams have adopted these parameterizations and include them in their data analysis. We have developed software for the quantitative interpretation of UVI and VIS images. After the launch of the POLAR satellite we used the image data from UVI and VIS in combination with ground based data from SuperDARN and incoherent scatter radars, magnetometers, and in situ observations from DMSP and NOAA satellites to characterize the state of the ionosphere. A number of event studies have been carried out in cooperation with other CGS theory teams.

Source record↗
NASA NTRS1994

Atmospheric, Climatic, and Environmental Research

The climate and atmospheric modeling project involves analysis of basic climate processes, with special emphasis on studies of the atmospheric CO2 and H2O source/sink budgets and studies of the climatic role Of CO2, trace gases and aerosols. These studies are carried out, based in part on use of simplified climate models and climate process models developed at GISS. The principal models currently employed are a variable resolution 3-D general circulation model (GCM), and an associated "tracer" model which simulates the advection of trace constituents using the winds generated by the GCM.

Broecker, Wallace S.↗
NASA NTRS1998

The Origin of Antarctic Precipitation: A Modeling Approach

Isotope concentrations in polar ice cores have long been used to estimate paleotemperatures. Underlying the use of this "isotope paleothermometer" is the assumption that the relationship between surface temperature and isotope concentration over time at a single geographical point is the same as that observed over space during the present-day climate. The validity of this assumption may in fact be compromised by several factors related to climate change. The specific factor studied in this paper involves the evaporative sources for polar precipitation. Climatic changes in the relative strengths of these sources would imply a need for a recalibration of the paleothermometer. To quantify such changes, we performed two GCM simulations, one of present-day climate and the other of the climate during the Last Glacial Maximum (LGM), roughly 18000 years ago. Evaporative sources of Antarctic precipitation were established using special tracer diagnostics. Results suggest that polar precipitation during the LGM does indeed consist of (relatively) more water from tropical oceans, a direct reflection of the LGM's increased equator-to-pole temperature gradient and its increased sea ice extent, which reduces high latitude evaporation. This result implies that an uncalibrated ice core paleothermometer would produce LGM temperatures that are biased slightly low. Because LGM boundary conditions are still under debate, we performed a third GCM simulation using a modified set of LGM boundary conditions. Using this simulation gives some qualitatively similar results, though the tropical contribution is not quite as high. Uncertainties in the LGM boundary conditions does hamper success in calibrating the paleothermometer.

Delaygue, Gilles↗
NASA NTRS1999

Two-D Representation of Three-D Modeling of Atmospheric Effects of Aircraft

Our principal accomplishments during the period of this research grant were as follows: (1) development of 2-D Transport Coefficients from a 3-D GCM; (2) comparison of various 2-D Treatments of Temperature; Asymmetries for Use in 2-D Models; (3) development and Testing of Convective Flux Formulation from a GCM for Use in a 2-D Model; and (4) investigation of Bromine Balance in the Upper Troposphere and Lower Stratosphere.

Source record↗
NASA NTRS1999

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 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 c) = 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 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 c) = K(sub h). Indeed, the turbulence model yelds a salt diffusivity that depends on R(sub i) and R(sub rho) = ((beta)(delta)S/(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.↗