Search NASA⌕ Search

SEARCH · Search NASA

Results for “Ocean Modeling”

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 973 records · Page 54

Ocean Turbulence I: One-Point Closure Model Momentum and Heat Vertical Diffusivities

Since the early forties, one-point turbulence closure models have been the canonical tools used to describe turbulent flows in many fields. In geophysics, Mellor and Yamada applied such models using the 1980 state-of-the art. Since then, no improvements were introduced to alleviate two major difficulties: 1) closure of the pressure correlations, which affects the correct determination of the critical Richardson number Ri(sub cr) above which turbulent mixing is no longer possible and 2) the need to express the non-local third-order moments (TOM) in terms of lower order moments rather than via the down-gradient approximation as done thus far, since the latter seriously underestimates the TOMs. Since 1) and 2) are still being dealt with adjustable parameters which weaken the credibility of the models, alternative models, not based on turbulence modeling, have been suggested. The aim of this paper is to show that new information, partly derived from the newest 2-point closure model discussed, can be used to solve these shortcomings. The new one-point closure model, which in its simplest form is algebraic and thus simple to implement, is first shown to reproduce a variety of data. Then, it is used in a Ocean-General Circulation Model (O-GCM) where it reproduces well a large variety of ocean data. While phenomenological models are specifically tuned to ocean turbulence, the present model is not. It is first tested against laboratory data on stably stratified flows and then used in an O-GCM. It is more general, more predictive and more resilient, e.g., it can incorporate phenomena like wave-breaking at the surface, salinity diffusivity, non-locality, etc. One important feature that naturally comes out of the new model is that the predicted Richardson critical value Ri(sub cr) is Ri (sub cr approx. = 1) in agreement with both Large Eddy Simulations (LES) and empirical evidence while all previous models predicted Ri (sub cr approx. = 0.2) which led to a considerable underestimate of the extent of turbulent mixing and thus to an incorrect mixed layer depth. The predicted temperature and salinity profiles (vs. depth) are presented and compared with those of the Kolmogorov-Petruvsky-Piskunuv (KPP) model and Levitus data.

Canuto, V. M.↗

Exploring and Analyzing Climate Variations Online by Using NASA MERRA-2 Data at GES DISC

NASA Giovanni (Goddard Interactive Online Visualization ANd aNalysis Infrastructure) (http:giovanni.sci.gsfc.nasa.govgiovanni) is a web-based data visualization and analysis system developed by the Goddard Earth Sciences Data and Information Services Center (GES DISC). Current data analysis functions include Lat-Lon map, time series, scatter plot, correlation map, difference, cross-section, vertical profile, and animation etc. The system enables basic statistical analysis and comparisons of multiple variables. This web-based tool facilitates data discovery, exploration and analysis of large amount of global and regional remote sensing and model data sets from a number of NASA data centers. Long term global assimilated atmospheric, land, and ocean data have been integrated into the system that enables quick exploration and analysis of climate data without downloading, preprocessing, and learning data. Example data include climate reanalysis data from NASA Modern-Era Retrospective analysis for Research and Applications, Version 2 (MERRA-2) which provides data beginning in 1980 to present; land data from NASA Global Land Data Assimilation System (GLDAS), which assimilates data from 1948 to 2012; as well as ocean biological data from NASA Ocean Biogeochemical Model (NOBM), which provides data from 1998 to 2012. This presentation, using surface air temperature, precipitation, ozone, and aerosol, etc. from MERRA-2, demonstrates climate variation analysis with Giovanni at selected regions.

knowledge base↗

Symmetric Equations on the Surface of a Sphere as Used by Model GISS:IB

Standard vector calculus formulas of Cartesian three space are projected onto the surface of a sphere. This produces symmetric equations with three nonindependent horizontal velocity components. Each orthogonal axis has a velocity component that rotates around its axis (eastward velocity rotates around the north–south axis) and a specific angular momentum component that is the product of the velocity component multiplied by the cosine of axis’ latitude. Angular momentum components align with the fixed axes and simplify several formulas, whereas the rotating velocity components are not orthogonal and vary with location. Three symmetric coordinates allow vector resolution and calculus operations continuously over the whole spherical surface, which is not possible with only two coordinates. The symmetric equations are applied to one-layer shallow water models on cubed-sphere and icosahedral grids, the latter being computationally simple and applicable to an ocean domain. Model results are presented for three different initial conditions and five different resolutions.

Icosahedral grid↗

Orbital, Rotational and Climatic Interactions: Energy Dissipation and Angular Momentum Exchange in the Earth-Moon System

A numerical ocean tide model has been developed and tested using highly accurate TOPEX/Poseidon (T/P) tidal solutions. The hydrodynamic model is based on time stepping a finite difference approximation to the non-linear shallow water equations. Two novel features of our implementation are a rigorous treatment of self attraction and loading (SAL), and a physically based parameterization for internal tide (IT) radiation drag. The model was run for a range of grid resolutions, and with variations in model parameters and bathymetry. For a rational treatment of SAL and IT drag, the model run at high resolution (1/12 degree) fits the T/P solutions to within 5 cm RMS in the open ocean. Both the rigorous SAL treatment and the IT drag parameterization are required to obtain solutions of this quality. The sensitivity of the solution to perturbations in bathymetry suggest that the fit to T/P is probably now limited by errors in this critical input. Since the model is not constrained by any data, we can test the effect of dropping sea-level to match estimated bathymetry from the last glacial maximum (LGM). Our results suggest that the 100 m drop in sea-level in the LGM would have significantly increased tidal amplitudes in the North Atlantic, and increased overall tidal dissipation by about 40%. However, details in tidal solutions for the past 20 ka are sensitive to the assumed stratification. IT drag accounts for a significant fraction of dissipation, especially in the LGM when large areas of present day shallow sea were exposed, and this parameter is poorly constrained at present.

Egbert, Gary D.↗

Report on activities and findings under DOE grant “Collaborative research. An Interactive Multi-Model for Consensus on Climate Change”

The project takes a hierarchical approach. The supermodeling scheme was first studied exhaustively with simple systems of ordinary differential equations. Results were described in detail in the previous report. The principal findings were that 1) for highly non-linear systems, such as Lorenz-63, including systems which describe phenomena on very different (atmosphere/ocean) times scales, supermodeling is far superior to any form of output-averaging; 2) negative coefficients can be used to advantage in situations where all models err in the same way, but to different degrees; 3) an interesting variant of supermodeling, “weighted supermodeling”, is the limiting case where inter-model nudging coefficients in the originally conceived “connected supermodel” become infinite, but with fixed ratios, corresponding to a direct combination of the tendencies that appear in corresponding equations for the alternative models; 4) noise is useful for avoiding local optima in training the inter-model coefficients in the supermodel. The supermodeling scheme was then investigated with simple quasigeostrophic (QG) models. As described in the previous report, it was found that QG models on a sphere can be coupled most efficaciously by working in a basis which captures the most variance, rather than the most instability, a somewhat unexpected result that still deserves scrutiny in a broader context. Further studies (since the last report) with QG channel models addressed the central question of when supermodeling is superior to output averaging in situations where nonlinearites are less extreme than with the ODEs initially studied. It was found that for realistic variations in a parameter in the QG model, output averaging is sufficient to capture all but the most subtle quantitative and qualitative behavior. Supermodeling helps when qualitative differences between the models result from unrealistically large parameter differences, or when very detailed spatial structure of the modes of variability are of interest. Therefore, the scheme may still be useful in the case of full climate models with qualitatively different parametrization schemes. A supermodel was constructed from the intermediate-complexity SPEEDO model, a primitive equation model with ocean and land. Versions defined by different parameter choices, in a realistic range, were connected and the coefficients trained. Some improvement was found as compared to output averaging. The learning algorithm used thus far gives sub-optimal, but still useful results when the CO 2 level and other parameters are varied. Spatial structure remains to be studied. The first use of supermodeling with full climate models has been with variants of the ECHAM model that use different convection schemes. As yet the models are only connected at the ocean-atmosphere interface, where weighted combinations of fluxes from the two atmospheres are passed to a common ocean, and the weights adapted during a training period. The supermodel was surprisingly successful at avoiding unrealistic features such as the double-ITCZ (Intertropical Convergence Zone), a problem that arises in both of the two models run separately. The supermodels constructed thus far have not identified dynamical regime shifts in future climate. Thus the planned connection with the work of Tsonis on the relationship between regime shifts and synchronization/de-synchronization among the major climate modes (see U. Wisconsin report) has not yet been made. However the network analysis of the climate system, in observations and models, that was done in conjunction with that study, shows that models differ strongly from one another and from observations in regard to the dynamical structure described by correlation networks [Steinhaeuser and Tsonis 2013], providing a further justification for supermodeling. Toward a general software framework for supermodeling, three versions of CAM (the Community Atmosphere Model) at NCAR were configured for inter-model nudging using the DART (Data Assimilation Research Testbed) capability to stop and re-start models in synchrony. It was clearly established that the inter-model nudging adds almost no computational burden to the runs, but there appears to be a problem with the re-initialization software that is still being debugged. Publications: Several papers were published on the basic idea of the interactive multi-model (supermodel) including demonstrations with low-order ODEs. The last of these, a semi-philosophical review paper on the relevance of synchronization generally, encountered considerable resistance but was finally published in Entropy [Duane 2015]. A paper on the ECHAM/COSMOS supermodel, containing the most promising results so far [Shen et al. 2015] is presently under review.

58 GEOSCIENCES↗

A real two-phase submarine debris flow and tsunami

The general two-phase debris flow model proposed by Pudasaini is employed to study subaerial and submarine debris flows, and the tsunami generated by the debris impact at lakes and oceans. The model, which includes three fundamentally new and dominant physical aspects such as enhanced viscous stress, virtual mass, and generalized drag (in addition to buoyancy), constitutes the most generalized two-phase flow model to date. The advantage of this two-phase debris flow model over classical single-phase, or quasi-two-phase models, is that the initial mass can be divided into several parts by appropriately considering the solid volume fraction. These parts include a dry (landslide or rock slide), a fluid (water or muddy water; e.g., dams, rivers), and a general debris mixture material as needed in real flow simulations. This innovative formulation provides an opportunity, within a single framework, to simultaneously simulate the sliding debris (or landslide), the water lake or ocean, the debris impact at the lake or ocean, the tsunami generation and propagation, the mixing and separation between the solid and fluid phases, and the sediment transport and deposition process in the bathymetric surface. Applications of this model include (a) sediment transport on hill slopes, river streams, hydraulic channels (e.g., hydropower dams and plants); lakes, fjords, coastal lines, and aquatic ecology; and (b) submarine debris impact and the rupture of fiber optic, submarine cables and pipelines along the ocean floor, and damage to offshore drilling platforms. Numerical simulations reveal that the dynamics of debris impact induced tsunamis in mountain lakes or oceans are fundamentally different than the tsunami generated by pure rock avalanches and landslides. The analysis includes the generation, amplification and propagation of super tsunami waves and run-ups along coastlines, debris slide and deposition at the bottom floor, and debris shock waves. It is observed that the submarine debris speed can be faster than the tsunami speed. This information can be useful for early warning strategies in the coastal regions. These findings substantially increase our understanding of complex multi-phase systems and multi-physics and flows, and allows for the proper modeling of landslide and debris induced tsunami, the dynamics of turbidity currents and sediment transport, and the associated applications to hazard mitigation, geomorphology and sedimentology.

42 ENGINEERING↗

Impact of Pacific Ocean Heatwaves on Phytoplankton Community Composition

Since 2013, marine heatwaves have become recurrent throughout the equatorial and northeastern Pacific Ocean and are expected to increase in intensity relative to historic norms. Among the ecological ramifications associated with these high temperature anomalies are increased mortality of higher trophic organisms such as marine mammals and seabirds, which are likely triggered by changes in the composition of phytoplankton, the base of the marine trophic food web. Here, we assimilated satellite ocean color data into an ocean biogeochemical model to describe changes in the abundance of phytoplankton functional types (PFTs) during the last decade's (2010s) warm anomalies in the equatorial and northeastern Pacific Ocean. We find important changes associated with the "Blob" warm anomaly in the Gulf of Alaska, where reduced silica supply led to a switch in community composition from diatoms to dinoagellates, resulting in an increase in surface ocean chlorophyll during the Summer-Fall of 2014. A more dramatic change was observed in the equatorial Pacific, where the extreme warm conditions of the 2016 El Niño resulted in a major decline of about 40 % in surface chlorophyll, which was associated with a nearly total collapse in diatoms.

Pacific Ocean↗

Thermo-Hydro-Chemical Simulation of Mid-Ocean Ridge Hydrothermal Systems: Static 2D Models and Effects of Paleo-Seawater Chemistry

Decades of research have resulted in characterization of the ocean floor manifestations of mid-ocean ridge (MOR) hydrothermal systems, yet numerical models accounting for the connections between heat transfer, hydrology and geochemistry have been slow to develop. The Thermo-hydro-chemical code ToughReact can be used to describe the coupled effects of fluid flow, heat transfer, and fluid-rock chemical interactions that occur in MOR systems. We describe the results of 2-dimensional simulations of steady state flow in fractured diabase with mineral-fluid chemical reactions. Basal heating and specified permeability yield maximum temperature of 400°C. Total fluid flux and high fracture flow velocities are in accord with observations. Fluid chemistry, mineralogical changes and 87 Sr/ 86 Sr ratios can be compared to observations to assess and calibrate models. Simulated high temperature fracture fluids have Mg and SO 4 near zero, elevated Ca and 87 Sr/ 86 Sr of about 0.7040. Total alteration is 10%–50% for simple models of spreading. Anhydrite forms mainly near the base of the upwelling zone and results in substantial local fracture porosity reduction. A calibrated model is used to predict how Sr isotopes and other features of altered oceanic crust would be different in the Cretaceous (95 Ma) early Proterozoic (1,800 Ma) and Archean (3,800 Ma), when seawater may have had high Ca and Sr concentrations, lower pH, higher temperature, and lower Na, Mg, and SO 4 . The simulations are offered as a start on what ultimately may require a longer-term community effort to better understand the role of MOR thermo-hydro-chemical systems in Earth evolution.

58 GEOSCIENCES↗

The Dependence of Climate Sensitivity on the Meridional Distribution of Radiative Forcing

Abstract This study investigates how climate sensitivity depends upon the spatial pattern of radiative forcing. Sensitivity experiments using a coupled ocean‐atmosphere model were conducted by adding anomalous incoming solar radiation over the entire globe, Northern Hemisphere mid‐latitudes, Southern Ocean, and tropics. The varied forcing patterns led to highly divergent climate sensitivities. Specifically, the climate is nearly twice as sensitive to Southern Ocean forcing as tropical forcing. Strong coupling between the surface and free troposphere in the tropics increases the inversion strength, leading to smaller cloud feedback in the tropical forcing experiments. In contrast, the extratropics exhibit weaker coupling, a decrease or near‐zero change in the inversion strength, and strong positive cloud feedback. These results contrast with the conventional SST‐pattern effect in which tropical surface temperature changes regulate climate sensitivity. They also have important implications for other potentially asymmetric forcings, such as those from geoengineering, volcanic eruptions, and paleoclimatic changes.

54 ENVIRONMENTAL SCIENCES↗

Using the NASA Giovanni DICCE Portal to Investigate Land-Ocean Linkages with Satellite and Model Data

Data-enhanced Investigations for Climate Change Education (DICCE), a NASA climate change education project, employs the NASA Giovanni data system to enable teachers to create climate-related classroom projects using selected satellite and assimilated model data. The easy-to-use DICCE Giovanni portal (DICCE-G) provides data parameters relevant to oceanic, terrestrial, and atmospheric processes. Participants will explore land-ocean linkages using the available data in the DICCE-G portal, in particular focusing on temperature, ocean biology, and precipitation variability related to El Ni?o and La Ni?a events. The demonstration includes the enhanced information for educators developed for the DICCE-G portal. The prototype DICCE Learning Environment (DICCE-LE) for classroom project development will also be demonstrated.

Acker, James G.↗

A simple object-oriented and open-source model for scientific and policy analyses of the global climate system – Hector v1.0

Simple climate models play an integral role in the policy and scientific communities. They are used for climate mitigation scenarios within integrated assessment models, complex climate model emulation, and uncertainty analyses. Here we describe Hector v1.0, an open source, object-oriented, simple global climate carbon-cycle model. This model runs essentially instantaneously while still representing the most critical global-scale earth system processes. Hector has a three-part main carbon cycle: a one-pool atmosphere, land, and ocean. The model's terrestrial carbon cycle includes primary production and respiration fluxes, accommodating arbitrary geographic divisions into, e.g., ecological biomes or political units. Hector actively solves the inorganic carbon system in the surface ocean, directly calculating air–sea fluxes of carbon and ocean pH. Hector reproduces the global historical trends of atmospheric [CO 2 ], radiative forcing, and surface temperatures. The model simulates all four Representative Concentration Pathways (RCPs) with equivalent rates of change of key variables over time compared to current observations, MAGICC (a well-known simple climate model), and models from the 5th Coupled Model Intercomparison Project. Hector's flexibility, open-source nature, and modular design will facilitate a broad range of research in various areas.

54 ENVIRONMENTAL SCIENCES↗

The self-consistent dynamic pole tide in global oceans

The dynamic pole tide is characterized in a self-consistent manner by means of introducing a single nondifferential matrix equation compatible with the Liouville equation, modelling the ocean as global and of uniform depth. The deviations of the theory from the realistic ocean, associated with the nonglobality of the latter, are also given consideration, with an inference that in realistic oceans long-period modes of resonances would be increasingly likely to exist. The analysis of the nature of the pole tide and its effects on the Chandler wobble indicate that departures of the pole tide from the equilibrium may indeed be minimal.

Dickman, S. R.↗

Ocean Turbulence: Two-Point Closure Model Momentum, Heat and Salt Vertical Diffusivities in the Presence of Shear - Paper 3

In papers 1 and 2 we have presented the results of the most updated 1-point closure model for the turbulent vertical diffusivities of momentum, heat and salt, K(sub m,h,s). In this paper, we derive the analytic expressions for K(sub m,h,s) using a new 2-point closure model that has recently been developed and successfully tested against some approx. 80 turbulence statistics for different flows. The new model has no free parameters. The expressions for K(sub m, h. s) are analytical functions of two stability parameters: the Turner number R(sub rho) (salinity gradient/temperature gradient) and the Richardson number R(sub i) (temperature gradient/shear). The turbulent kinetic energy K and its rate of dissipation may be taken local or non-local (K-epsilon model). Contrary to all previous models that to describe turbulent mixing below the mixed layer (ML) have adopted three adjustable "background diffusivities" for momentum. heat and salt, we propose a model that avoids such adjustable diffusivities. We assume that below the ML, K(sub m,h,s) have the same functional dependence on R(sub i) and R(sub rho) derived from the turbulence model. However, in order to compute R(sub i) below the ML, we use data of vertical shear due to wave-breaking measured by Gargett et al. (1981). The procedure frees the model from adjustable background diffusivities and indeed we use the same model throughout the entire vertical extent of the ocean. Using the new K(sub m,h, s), we run an O-GCM and present a variety of results that we compare with Levitus and the KPP model. Since the traditional 1-point (used in papers 1 and 2) and the new 2-point closure models used here represent different modeling philosophies and procedures, testing them in an O-GCM is indispensable. The basic motivation is to show that the new 2-point closure model gives results that are overall superior to the 1-point closure in spite of the fact that the latter rely on several adjustable parameters while the new 2-point closure has none. After the extensive comparisons presented in papers 1 and 2, we conclude that the new model presented here is overall superior for it not only is parameter free but also 2 because is part of a more general turbulence model that has been previously successfully tested on a wide variety of other types of turbulent flows.

Canuto, V. M.↗

The black carbon cycle and its role in the Earth system

Black carbon (BC) is produced by incomplete combustion of biomass by wildfires and burning of fossil fuels. BC is environmentally persistent over centuries to millennia, sequestering carbon in marine and terrestrial environments. However, its production, storage and dynamics, and therefore its role in the broader carbon cycling during global change, are poorly understood. Here, in this Review, we discuss BC cycling across the land-to-ocean continuum. Wildfires are the main source of BC, producing 128 ± 84 teragrams per year. Negative climate–BC feedbacks could arise as wildfire increases with anthropogenic warming, producing more BC, which in turn will sequester carbon, but the magnitude of these effects are unknown. Most BC is stored in terrestrial systems with some transported to the ocean via rivers and the atmosphere. However, the oceanic BC budget is not balanced, with known BC removal fluxes exceeding BC inputs. We demonstrate these observed inconsistencies using a simple ocean box model, which highlights key areas of future research. Measurements of BC mineralization and export rates along the land-to-ocean continuum and quantification of previously unexplored sources of oceanic BC are needed to close the global BC budget.

aquatic↗

Diurnal variability of the hydrologic cycle in a general circulation model

In the present Colorado State University GCM simulation-based analysis of the diurnal and semidiurnal variability of precipitation, precipitable water, evaporation, cloudiness, horizontal moisture flux convergence, and cloud radiative forcing, a realistic afternoon precipitation maximum is obtained over land in warm rainy regions, as well as an early morning maximum over the oceans. The model has been further used to investigate the bases for the oceanic diurnal-precipitation cycle; the results thus obtained indicate that such an oceanic cycle occurs even in the absence of neighboring continents, and tends to have a morning maximum, although the observed phenomenon is generally stronger than the results indicate.

Randall, David A.↗