Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “ocean model”

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 325 records · Page 18

Petrology and Geochemistry of Unbrecciated Harzburgitic Diogenite MIL 07001: A Window Into Vestan Geological Evolution

There is a strong case that asteroid 4 Vesta is the parent of the howardite, eucrite and diogenite (HED) meteorites. Models developed for the geological evolution of Vesta can satisfy the compositions of basaltic eucrites that dominate in the upper crust. The bulk compositional characteristics of diogenites - cumulate harzburgites and orthopyroxenites from the lower crust - do not fit into global magma ocean models that can describe the compositions of basaltic and cumulate eucrites. Recent more detailed formation models do make provision for a more complicated origin for diogenites, but this model has yet to be completely vetted. Compositional studies of bulk samples has led to the hypothesis that many diogenites were formed late by interaction of their parent melts with a eucritic crust, but those observations may alternatively be explained by subsolidus equilibration of trace elements between orthopyroxene and plagioclase and Ca-phosphate in the rocks. Differences in radiogenic Mg-26 content between diogenites and eucrites favors early formation of the former, not later formation. Understanding the origin of diogenites is crucial for understanding the petrologic evolution of Vesta. We have been doing coordinated studies of a suite of diogenites including petrologic investigations, bulk rock major and trace element studies, and in situ trace element analyses of orthopyroxene. Here we will focus on an especially unusual, and potentially key, diogenite, MIL 07001.

Mittlefehldt, D. W.↗

High-Resolution Global Ocean Simulation in the SWOT ERA: Introducing ECCO Llc4320v2

The Surface Water and Ocean Topography (SWOT) mission is opening up new opportunities to evaluate and drive improvement in km-scale global ocean and atmosphere-ocean models via its simultaneous swath observations of sea-surface height, wind speed, and significant wave height down to 2 km scales globally. Motivated by this, we have initiated a collaborative effort at NASA Ames to run a new global ocean-sea-ice simulation in MITgcm at a nominal 1/48-deg resolution of the year 2023 that enables direct comparisons with SWOT observations. New features relative to a prior simulation at a similar resolution of the year 2012 include much more realistic tidal forcing, thermodynamic coupling in ice-shelf cavities, high-frequency hourly atmospheric forcing, much improved river runoff, double the vertical resolution in the deep ocean (173 vs 90 levels), updated bathymetry including improvements in finescale bathymetry derived from SWOT, and updated vertical mixing scheme parameters that improve the equatorial Pacific circulation.

ECCO↗

GEOS Atmospheric Model: Challenges at Exascale

The Goddard Earth Observing System (GEOS) model at NASA's Global Modeling and Assimilation Office (GMAO) is used to simulate the multi-scale variability of the Earth's weather and climate, and is used primarily to assimilate conventional and satellite-based observations for weather forecasting and reanalysis. In addition, assimilations coupled to an ocean model are used for longer-term forecasting (e.g., El Nino) on seasonal to interannual times-scales. The GMAO's research activities, including system development, focus on numerous time and space scales, as detailed on the GMAO website, where they are tabbed under five major themes: Weather Analysis and Prediction; Seasonal-Decadal Analysis and Prediction; Reanalysis; Global Mesoscale Modeling, and Observing System Science. A brief description of the GEOS systems can also be found at the GMAO website. GEOS executes as a collection of earth system components connected through the Earth System Modeling Framework (ESMF). The ESMF layer is supplemented with the MAPL (Modeling, Analysis, and Prediction Layer) software toolkit developed at the GMAO, which facilitates the organization of the computational components into a hierarchical architecture. GEOS systems run in parallel using a horizontal decomposition of the Earth's sphere into processing elements (PEs). Communication between PEs is primarily through a message passing framework, using the message passing interface (MPI), and through explicit use of node-level shared memory access via the SHMEM (Symmetric Hierarchical Memory access) protocol. Production GEOS weather prediction systems currently run at 12.5-kilometer horizontal resolution with 72 vertical levels decomposed into PEs associated with 5,400 MPI processes. Research GEOS systems run at resolutions as fine as 1.5 kilometers globally using as many as 30,000 MPI processes. Looking forward, these systems can be expected to see a 2 times increase in horizontal resolution every two to three years, as well as less frequent increases in vertical resolution. Coupling these resolution changes with increases in complexity, the computational demands on the GEOS production and research systems should easily increase 100-fold over the next five years. Currently, our 12.5 kilometer weather prediction system narrowly meets the time-to-solution demands of a near-real-time production system. Work is now in progress to take advantage of a hybrid MPI-OpenMP parallelism strategy, in an attempt to achieve a modest two-fold speed-up to accommodate an immediate demand due to increased scientific complexity and an increase in vertical resolution. Pursuing demands that require a 10- to 100-fold increases or more, however, would require a detailed exploration of the computational profile of GEOS, as well as targeted solutions using more advanced high-performance computing technologies. Increased computing demands of 100-fold will be required within five years based on anticipated changes in the GEOS production systems, increases of 1000-fold can be anticipated over the next ten years.

ESMF↗

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↗

Geosat-derived sea level and surface current anomalies in the equatorial Pacific during the 1986-1989 El Nino and La Nina

Equatorial wave dynamics are essential in most oceanic models for reproducing aspects of the El Nino-Southern Oscillation (ENSO) phenomenon. Observational evidence of first baroclinic equatorial Kelvin and first symmetric meridional Rossby waves is found in Geosat-derived sea level anomalies (SLA) and surface zonal current anomalies (ZCA) in the equatorial Pacific ocean during the 1986-1987 El Nino and the ensuing 1988-1989 La Nina. This was made possible after extensive quality control and specific processing of the recently improved Geosat geophysical data records pertaining to the 17-day Exact Repeat Mission. In particular, the processing was made so that the Geosat-derived ZCA best fit near-surface zonal currents from three equatorial moorings at 165 deg E, 140 deg W, and 110 deg W. The Geosat-derived SLA and ZCA are decomposed into first baroclinic equatorial Kelvin and gravest Rossby modes. The emphasis is then put on the chronology of the ZCA at the equator where the currents are the most energetic and where Kelvin and first symmetric Rossby waves explain most of the variance in ZCA, in similar proportion. The 1986-1987 El Nino is mostly characterized by a strong downwelling Kelvin wave in December 1986, a series of downwelling Kelvin waves in March-October 1987, and a strong upwelling Rossby wave in March-September 1987. These waves are consistent with wind forcing, and all give rise to notable eastward ZCA for almost an entire year. During El Nino-La Nina transition period there is the occurrence of two downwelling Rossby waves originating from the eastern Pacific at times of favorable wind forcing. These Rossby waves switch the basin-wide ZCA from eastward to westward, terminating the warm event. Then, the 1988-1989 La Nina shows the generation of a series of upwelling Kelvin waves which are clearly identified in March/April and September/October 1988. These upwelling Kelvin waves are apparently related to the enhencement of the trade wind in the western half of the basin. These waves maintain the basin-wide westward ZCA for more than a year, although the ZCA is slightly reduced ba front of upwelling Rossby waves in phase with the normal seasonal cycle. The Kelvin wave forcing and sign of kelvin wave contribution versus Rossby wave contribution to ZCA in the western Pacific are discussed, and we conclude that neither pleads in favor of the delayed action oscillator mechanism.

Delcroix, Thierry↗

Global Climate Impacts of Greenland and Antarctic Glacial Melt: A Comparative Study

A 10-member ensemble simulation using the NASA GISS-E2-1-G climate model shows a clear bifurcation in the Atlantic Meridional Overturning Circulation (AMOC) strength under the SSP2-4.5 future scenario, leading to 8 stable AMOC-on and 2 stable AMOC-off climate states after 400 years of integration, despite identical forcing of each ensemble member. Ocean models have previously shown a tipping point where sufficient forcing (e.g. by freshwater inputs) can shut down the AMOC. A variety of fully coupled models have demonstrated this, either through hosing experiments or increased precipitation and warming at high latitudes due to increased concentrations of greenhouse gases. In the GISS simulations, there are no external freshwater perturbations. The bifurcation arises freely in the coupled system and is the result of stochastic variability associated with sea-ice transport and melting in the Irminger Sea following a slowing of the GHG concentration increase at the end of the twenty-first century. We believe this is the first time that a CMIP-class model has shown such a bifurcation across an initial condition ensemble.

GISS ModelE2.1↗

On a generating mechanism for Yanai waves and the 25-day oscillation

A spectral Chebyshev-collocation method applied to the linear, 1.5 layer reduced-gravity ocean model equations is used to study the dynamics of Yanai (or mixed Rossby-gravity) wave packets. These are of interest because of the observations of equatorial instability waves (which have the characteristics of Yanai waves) and their role in the momentum and heat budgets in the tropics. A series of experiments is performed to investigate the generation of the waves by simple cross-equatorial wind stress forcings in various configurations and the influence of a western boundary on the waves. They may be generated in the interior ocean as well as from a western boundary. The observations from all the oceans indicate that the waves have a preferential period and wavelength of around 25 days and 1000 km respectively. These properties are also seen in the model results and a plausible explanation is provided as being due to the dispersive properties of Yanai waves.

Kelly, Brian G.↗

Sensitivity of Ocean Reflectance Inversion Models for Identifying and Discriminating Between Phytoplankton Functional Groups

The daily, synoptic images provided by satellite ocean color instruments provide viable data streams for observing changes in the biogeochemistrY of marine ecosystems. Ocean reflectance inversion models (ORMs) provide a common mechanism for inverting the "color" of the water observed a satellite into marine inherent optical properties (lOPs) through a combination of empiricism and radiative transfer theory. lOPs, namely the spectral absorption and scattering characteristics of ocean water and its dissolved and particulate constituents, describe the contents of the upper ocean, information critical for furthering scientific understanding of biogeochemical oceanic processes. Many recent studies inferred marine particle sizes and discriminated between phytoplankton functional groups using remotely-sensed lOPs. While all demonstrated the viability of their approaches, few described the vertical distributions of the water column constituents under consideration and, thus, failed to report the biophysical conditions under which their model performed (e.g., the depth and thickness of the phytoplankton bloom(s)). We developed an ORM to remotely identifY Noctiluca miliaris and other phytoplankton functional types using satellite ocean color data records collected in the northern Arabian Sea. Here, we present results from analyses designed to evaluate the applicability and sensitivity of the ORM to varied biophysical conditions. Specifically, we: (1) synthesized a series of vertical profiles of spectral inherent optical properties that represent a wide variety of bio-optical conditions for the northern Arabian Sea under aN Miliaris bloom; (2) generated spectral remote-sensing reflectances from these profiles using Hydrolight; and, (3) applied the ORM to the synthesized reflectances to estimate the relative concentrations of diatoms and N Miliaris for each example. By comparing the estimates from the inversion model to those from synthesized vertical profiles, we were able to identifY those bio-optical conditions under which the inversion model performs both well and poorly.

Werdell, P. Jeremy↗

Net diffusivity in ocean general circulation models with nonuniform grids

The numerical vertical diffusivity K(num), embedded in a numerical ocean general circulation model with nonuniform vertical grid, is estimated. It is shown that in a downwelling region, K(num) is negative for a grid with grid size increasing with depth. When the grid size increment, or the downward vertical velocity, is large, K(num) may exceed the vertical diffusivity specified and may result in a negative effective vertical diffusivity. Therefore care needs to be taken to specify the vertical diffusivity in a numerical model with nonuniform grid, and a lower bound is generally imposed in order to avoid an unphysical negative value. Some possible effects of the negative effective diffusivity are discussed.

Yin, F. L.↗

Effects of High-frequency Wind Sampling on Simulated Mixed Layer Depth and Upper Ocean Temperature

Effects of high-frequency wind sampling on a near-global ocean model are studied by forcing the model with a 12 hourly averaged wind product and its 24 hourly subsamples in separate experiments. The differences in mixed layer depth and sea surface temperature resulting from these experiments are examined, and the underlying physical processes are investigated. The 24 hourly subsampling not only reduces the high-frequency variability of the wind but also affects the annual mean wind because of aliasing. While the former effect largely impacts mid- to high-latitude oceans, the latter primarily affects tropical and coastal oceans. At mid- to high-latitude regions the subsampled wind results in a shallower mixed layer and higher sea surface temperature because of reduced vertical mixing associated with weaker high-frequency wind. In tropical and coastal regions, however, the change in upper ocean structure due to the wind subsampling is primarily caused by the difference in advection resulting from aliased annual mean wind, which varies with the subsampling time. The results of the study indicate a need for more frequent sampling of satellite wind measurement and have implications for data assimilation in terms of identifying the nature of model errors.

mixed layered depth↗

Diagnosing Seasonal to Multi-decadal Phytoplankton Group Dynamics in a Highly Productive Coastal Ecosystem

The Santa Barbara Channel, CA (SBC) is a biodiverse marine ecosystem fueled largely by phytoplankton productivity, and the composition of the phytoplankton community influences the magnitude and fates of this productivity. Here, we create a 22-year monthly time series of phytoplankton biomarker pigment concentrations in the SBC by combining 12 years of high performance liquid chromatography phytoplankton pigment concentrations with bio-optical models and 10 additional years of bio-optical observations. The bio-optical models skillfully predict biomarker pigment concentrations representative of five distinct phytoplankton groups (PGs; diatoms, dinoflagellates, chlorophytes, prymnesiophytes, and picophytoplankton) and resolve seasonal responses to the annual upwelling-relaxation cycle for all PGs except the dinoflagellates. Our observations indicate that nanophytoplankton groups respond most rapidly to seasonal upwelling, followed by diatoms, and then by picophytoplankton as the water column stratifies in the summer. A Regional Ocean Modeling System (ROMS) solution is used to relate advection of different source waters to the observed PG dynamics. The ROMS simulation results suggest that, on seasonal time scales, pronounced cross-SBC differences in PG seasonality are related to cross-SBC differences in source waters. El Niño Southern Oscillation events drive interannual variability in the upwelling response of most PGs. On decadal time scales, dinoflagellate blooms are associated with the warm phase of the North Pacific Gyre Oscillation and anomalous advection of Southern California Bight source waters into the SBC. Taken together, our results provide a novel view of phytoplankton community succession in response to seasonal upwelling by considering the dynamics of pico- and nano-phytoplankton and suggest that regional surface ocean advection plays a substantial role in driving phytoplankton composition in the SBC.

Phytoplankton↗

Geographical representation of radial orbit perturbations due to ocean tides: Implications for satellite altimetry

In analogy to the geographical representation of the zeroth-order radial orbit perturbations due to the static geopotential, similar relationships have been derived for radial orbit perturbations due to the ocean tides. At each location these perturbations are seen to be coherent with the tide height variations. The study of this singularity is of obvious importance to the estimation of ocean tides from satellite altimeter data. We derive analytical expressions for the sensitivity of altimeter derived ocean tide models to the ocean tide force model induced errors in the orbits of the altimeter satellite. In particular, we focus on characterizing and quantifying the nonresonant tidal orbit perturbations, which cannot be adjusted into the empirical accelerations or radial perturbation adjustments commonly used during orbit determination and in altimeter data processing. As an illustration of the utility of this technique, we study the differences between a TOPEX/POSEIDON-derived ocean tide model and the Cartwright and Ray 1991 Geosat model. This analysis shows that nearly 60% of the variance of this difference for M(sub 2) can be explained by the Geosat radial orbit eror due to the omission of coefficients from the GEM-T2 background ocean tide model. For O(sub 1), K(sub 1), S(sub 2), and K(sub 2) the orbital effects account for approximately 10 to 40% of the variances of these differences. The utility of this technique to assessment of the ocean tide induced errors in the TOPEX/POSEIDON-derived tide models is also discussed.

Bettadpur, Srinivas V.↗

On Fofonoff's Mode

The classical barotropic ocean model of Veronis (1966) is reexamined to verify the possible existence of Fofonoff's (1954) steady free inertial mode by considering the balance where J is the Jacobian, eta is the streamfunction of the horizontal velocity, and R is a parameter of nonlinearity. A linear functional relationship between the potential vorticity and the streamfunction and inertial boundary layers along the ocean basin was suggested by Fofonoff. However, the numerical integration showed that the system is barotropically unstable. The typical instantaneous field is given. The statistical equilibrium state eta is also given. Although it resembles Fofonoff's solution, it is maintained by the Reynolds stress field induced by the eddies in addition to the wind stress and dissipation. The plot of eta versus R delta sq +y is given. It is shown that the Fofonoff's linear relation between the potential vorticity and the streamfunction does not hold uniformly over the entire basin. The above results show that the Fofonoff's solution can not be realized physically.

Merkine, L. O.↗

Integration of space and in situ observations to study global climate change

The use of model-based global data sets of atmospheric circulation for studying fundamental dynamical and physical processes is discussed, focusing on limitations of the available model-based data sets. Data from the Global Weather Experiment in 1979 were analyzed by two authorized level IIIb data centers in 1980 and in 1981. The analyses led to difference in data-sparse regions such as the tropics. Study areas which can be addressed by an internally-consistent long-term multivariate data set for the atmospheric circulation are considered, including mean climate, forcing for the ocean models, global hydrological cycle, atmospheric energetics, intraseasonal variability, land surface processes, and structure and variability of vertical velocity, divergence, and diabatic heating. It is concluded that the most comprehensive technique for integrating space and in situ observations to produce this type of data set would be a four-dimensional data assimilation system with a realistic physical model of the type employed in operational numerical weather prediction.

Bengtsson, L.↗

Establishing an IERS Sub-Center for Ocean Angular Momentum

With the objective of establishing the Special Bureau for the Oceans (SBO), a new archival center for ocean angular momentum (OAM) products, we have computed and analyzed a number of OAM products from several ocean models, with and without data assimilation. All three components of OAM (axial term related to length of day variations and equatorial terms related to polar motion) have been examined in detail, in comparison to the respective Earth rotation parameters. An 11+ year time series of OAM given at 5-day intervals has been made publicly available. Other OAM products spanning longer periods and with higher temporal resolution, as well as products calculated from ocean model/data assimilation systems, have been prepared and should become part of the SBO archives in the near future.

Ponte, Rui M.↗

E-Theatre

Mesoscale/severe storms and results from Terra Fritz Hasler (NASA/Goddard) will show the latest spectacular animations of the hurricanes and tropical storms: Floyd, Georges, Mitch, Bonnie etc. from Geostationary Operational Environmental Satellites (GOES) and Tropical Rain Measuring Mission (TRMM) supported by MMS 3D nested numerical model results. See dust storms in Africa and smoke plumes from fires in Mexico. See how High-Definition Television (HDTV) is revolutionizing entertainment and the way we communicate scientific results. Fly in from space and take a tour of Africa. Drop in on Albuquerque using 1-m Ikonos 'Spy Satellite' technology. See data from the new Earth Observing System (EOS) Terra satellite as well as Landsat 7, and SeaWiFS. See climate change in action with Global Land and Ocean productivity changes over the last 20 years. Remote sensing observations of Sea Surface Temperature (SST), El Nino, height, winds and ocean color from GOES, National Oceanic and Atmosphere Administration Advanced Very High Resolution Radiometer (AVHRR), Special Sensor Microwave Imager (SSMI), and SeaWiFS put in context with ten year LANL global ocean model. Compare the Gulf Stream observed by Terra with the simulations. Additional information is contained in the original extended abstract.

Hasler, Arthur F.↗

Cloud Property Retrieval and 3D Radiative Transfer

Cloud thickness and photon mean-free-path together determine the scale of "radiative smoothing" of cloud fluxes and radiances. This scale is observed as a change in the spatial spectrum of cloud radiances, and also as the "halo size" seen by off beam lidar such as THOR and WAIL. Such of beam lidar returns are now being used to retrieve cloud layer thickness and vertical scattering extinction profile. We illustrate with recent measurements taken at the Oklahoma ARM site, comparing these to the-dependent 3D simulations. These and other measurements sensitive to 3D transfer in clouds, coupled with Monte Carlo and other 3D transfer methods, are providing a better understanding of the dependence of radiation on cloud inhomogeneity, and to suggest new retrieval algorithms appropriate for inhomogeneous clouds. The international "Intercomparison of 3D Radiation Codes" or I3RC, program is coordinating and evaluating the variety of 3D radiative transfer methods now available, and to make them more widely available. Information is on the Web at: http://i3rc.gsfc.nasa.gov/. Input consists of selected cloud fields derived from data sources such as radar, microwave and satellite, and from models involved in the GEWEX Cloud Systems Studies. Output is selected radiative quantities that characterize the large-scale properties of the fields of radiative fluxes and heating. Several example cloud fields will be used to illustrate. I3RC is currently implementing an "open source" 3d code capable of solving the baseline cases. Maintenance of this effort is one of the goals of a new 3DRT Working Group under the International Radiation Commission. It is hoped that the 3DRT WG will include active participation by land and ocean modelers as well, such as 3D vegetation modelers participating in RAMI.

Cahalan, Robert F.↗

Role of Ocean Currents and Bottom Pressure Variability on Seasonal Polar Motion

Changes in the ocean angular momentum (OAM) components about the equatorial axes, either due to fluctuations in currents or bottom pressure (mass redistribution), can induce movements of the Earth's pole of rotation, commonly referred to as polar motion or wobble. Output from a 1 deg resolution ocean model is used to calculate the effective equatorial OAM functions chi(sub 1, sup O) and chi(sub 2, sup O), corresponding to polar motion excitation about the equatorial axis pointing to the Greenwich and 90 deg E meridians, respectively. Time series of chi(sup O) are combined with similar atmospheric series chi(sup A), computed from the National Centers for Environmental Prediction/National Center for Atmospheric Research reanalyses, to interpret the observed low-frequency polar motion excitation for the period 1985-1996. Results indicate that the oceans are a very important excitation source for the Chandler (approximately 433 days), annual, and semiannual wobbles, providing for much better amplitude and phase agreement with the observed excitation at these periods, in comparison with what is obtained when only the atmosphere is considered. Both oceanic mass and motion terms are found to be important but with mass signals having somewhat larger amplitudes. The role of regional variability in ocean currents and bottom pressure in contributing to chi(sup O) signals is quantified. Midlatitude regions (approximately 30 deg - 70 deg) figure prominently as places of strong local oceanic excitation signals. The North Pacific basin is found to be generally important for chi(sup O) excitation, while the Southern Ocean is important for both chi(sub 1, sup O) and chi(sub 2, sup O). The largest positive covariances of local with global chi(sup O) signals occur in the Kuroshio region near the western boundary of the North Pacific for chi(sub 1, sup O) and southwest of Australia for chi(sub 2, sup O).

Ponte, Rui M.↗