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

CMIP5 Historical Simulations (1850-2012) with GISS ModelE2

Observations of climate change during the CMIP5 extended historical period (1850-2012) are compared to trends simulated by six versions of the NASA Goddard Institute for Space Studies ModelE2 Earth System Model. The six models are constructed from three versions of the ModelE2 atmospheric general circulation model, distinguished by their treatment of atmospheric composition and the aerosol indirect effect, combined with two ocean general circulation models, HYCOM and Russell. Forcings that perturb the model climate during the historical period are described. Five-member ensemble averages from each of the six versions of ModelE2 simulate trends of surface air temperature, atmospheric temperature, sea ice and ocean heat content that are in general agreement with observed trends, although simulated warming is slightly excessive within the past decade. Only simulations that include increasing concentrations of long-lived greenhouse gases match the warming observed during the twentieth century. Differences in twentieth-century warming among the six model versions can be attributed to differences in climate sensitivity, aerosol and ozone forcing, and heat uptake by the deep ocean. Coupled models with HYCOM export less heat to the deep ocean, associated with reduced surface warming in regions of deepwater formation, but greater warming elsewhere at high latitudes along with reduced sea ice. All ensembles show twentieth-century annular trends toward reduced surface pressure at southern high latitudes and a poleward shift of the midlatitude westerlies, consistent with observations.

ocean heat content↗

The Impact of Satellite Sea Surface Salinity for Prediction of the Coupled Indo-Pacific System

We assess the impact of satellite sea surface salinity (SSS) observations on seasonal to interannual variability of tropical Indo-Pacific Ocean dynamics as well as on dynamical ENSO forecasts. Our coupled model is composed of a primitive equation ocean model for the tropical Indo-Pacific region that is coupled with the global SPEEDY atmospheric model (Molteni, 2003). The Ensemble Reduced Order Kalman Filter is used to assimilate observations to constrain dynamics and thermodynamics for initialization of the coupled model. The baseline experiment assimilates satellite sea level, SST, and in situ subsurface temperature and salinity observations. This baseline is then compared with experiments that additionally assimilate Aquarius (version 4.0) and SMAP (version 2.0) SSS. Twelve-month forecasts are initialized for each month from Sep. 2011 to Dec. 2016. We find that including satellite SSS significantly improves NINO 3.4 sea surface temperature anomaly validation after 1 out to 12 month forecast lead times. For initialization of the coupled forecast, the positive impact of SSS assimilation is brought about by surface freshening near the eastern edge of the western Pacific warm pool and density changes that lead to shallower mixed layer between 10 degrees South latitude-5 degrees North latitude. SST differences at initialization force wide-spread downwelling favorable curl over most of the tropical Pacific. Over an average forecast, SST remains warmer with SSS assimilation at the eastern edge of the warm pool. This warm SST propagates into the eastern Pacific and drags westerly wind anomalies eastward into the NINO 3.4 region. In addition, salting near the ITCZ (Intertropical Convergence Zone) leads to a deepening of the mixed layer and thermocline near 8 degrees North latitude. These patterns together lead to a funneling effect that provides the background state to amplify equatorial Kelvin waves. We show that the downwelling Kelvin waves are amplified by assimilating satellite SSS and lead to significantly improved forecasts particularly for the 2015 El Nino.

El Nino↗

The Impact of Satellite Sea Surface Salinity for Prediction of the Coupled Indo-Pacific System

We assess the impact of satellite sea surface salinity (SSS) observations on seasonal to interannual variability of tropical Indo-Pacific Ocean dynamics as well as on dynamical ENSO forecasts. Our coupled model is composed of a primitive equation ocean model for the tropical Indo-Pacific region that is coupled with the global SPEEDY atmospheric model (Molteni, 2003). The Ensemble Reduced Order Kalman Filter is used to assimilate observations to constrain dynamics and thermodynamics for initialization of the coupled model. The baseline experiment assimilates satellite sea level, SST, and in situ subsurface temperature and salinity observations. This baseline is then compared with experiments that additionally assimilate Aquarius (version 4.0) and SMAP (version 2.0) SSS. Twelve-month forecasts are initialized for each month from Sep. 2011 to Dec. 2016. We find that including satellite SSS significantly improves NINO3.4 sea surface temperature anomaly validation after 1 out to 12 month forecast lead times. For initialization of the coupled forecast, the positive impact of SSS assimilation is brought about by surface freshening near the eastern edge of the western Pacific warm pool and density changes that lead to shallower mixed layer between 10S-5N. SST differences at initialization force wide-spread downwelling favorable curl over most of the tropical Pacific. Over an average forecast, SST remains warmer with SSS assimilation at the eastern edge of the warm pool. This warm SST propagates into the eastern Pacific and drags westerly wind anomalies eastward into the NINO3.4 region. In addition, salting near the ITCZ leads to a deepening of the mixed layer and thermocline near 8N. These patterns together lead to a funneling effect that provides the background state to amplify equatorial Kelvin waves. We show that the downwelling Kelvin waves are amplified by assimilating satellite SSS and lead to significantly improved forecasts particularly for the 2015 El Nino.

SST↗

Small-scale physics of the ocean

Observations and theoretical models of small-scale phenomena in the oceans are reviewed, with a focus on progress during the period 1983-1986. Topics examined include surface layers, equatorial turbulence, off-equator mixed layers, the scaling of mixing, turbulence concepts, laboratory results, internal waves and mixing, rings, the nature of the bottom layer, double diffusion and intrusions, salt fingers, and biological interactions. Also discussed are developments in instrumentation (fast sampling profilers with upward-profiling capability, deep profilers, ship-motion correction, horizontal samplers, small submersibles, submarines, towed packages, conductivity sensors, dissolved-oxygen sensors, and acoustic Doppler current profilers) and goals for future research.

Caldwell, Douglas R.↗

Efficient Flowline Simulations of Ice Shelf-Ocean Interactions: Sensitivity Studies with a Fully Coupled Model

Thermodynamic flowline and plume models for the ice shelf-ocean system simplify the ice and ocean dynamics sufficiently to allow extensive exploration of parameters affecting ice-sheet stability while including key physical processes. Comparison between geophysically and laboratory-based treatments of ice-ocean interface thermodynamics shows reasonable agreement between calculated melt rates, except where steep basal slopes and relatively high ocean temperatures are present. Results are especially sensitive to the poorly known drag coefficient, highlighting the need for additional field experiments to constrain its value. These experiments also suggest that if the ice-ocean interface near the grounding line is steeper than some threshold, further steepening of the slope may drive higher entrainment that limits buoyancy, slowing the plume and reducing melting; if confirmed, this will provide a stabilizing feedback on ice sheets under some circumstances.

polar oceanography↗

Evolution of major mineral compositions and trace element abundances during fractional crystallization of a model lunar composition

The evolution of major mineral compositions and trace element abundances during perfect fractional crystallization of a model lunar magma ocean was calculated. The minerals in the model lunar composition were olivine, orthopyroxene, clinopyroxene, and plagioclase. Lunar bulk composition data, major mineral/melt equilibria data, and trace element partition data were taken from published sources. The results show olivine beginning to crystallize at 1380 C. Approximately 50% of the system crystallizes as olivine. From 50 to 60% solidification, orthopyroxene crystallizes alone. During the final 40% solidification, Ca-rich clinopyroxene and plagioclase crystallize together. Various changes in composition of all these minerals are also noted during the process. Concomitant evolution of major element abundances in the melt is followed along with that of trace element abundances. Consequences of the results for constraints on some aspects of the composition of the lunar magma ocean and of the primitive moon are discussed.

Drake, M. J.↗

Free oscillations in a climate model with ice-sheet dynamics

A study of stable periodic solutions to a simple nonlinear model of the ocean-atmosphere-ice system is presented. The model has two dependent variables: ocean-atmosphere temperature and latitudinal extent of the ice cover. No explicit dependence on latitude is considered in the model. Hence all variables depend only on time and the model consists of a coupled set of nonlinear ordinary differential equations. The globally averaged ocean-atmosphere temperature in the model is governed by the radiation balance. The reflectivity to incoming solar radiation, i.e., the planetary albedo, includes separate contributions from sea ice and from continental ice sheets. The major physical mechanisms active in the model are (1) albedo-temperature feedback, (2) continental ice-sheet dynamics and (3) precipitation-rate variations. The model has three-equilibrium solutions, two of which are linearly unstable, while one is linearly stable. For some choices of parameters, the stability picture changes and sustained, finite-amplitude oscillations obtain around the previously stable equilibrium solution. The physical interpretation of these oscillations points to the possibility of internal mechanisms playing a role in glaciation cycles.

Kallen, E.↗

Drift current in a thin surface layer of the ocean

A theoretical model is proposed for the drift current in the surface layer of the sea that considers molecular transfer exchange in the thermal conductivity coefficient for the thin layers of water and air near the ocean/atmosphere boundary as well as a layer of depth in the sea in which the coefficient rises from the molecular to turbulent values.

Shumilov, A. V.↗

A satellite technique for quantitatively mapping rainfall rates over the oceans

A theoretical model for calculating microwave radiative transfer in raining atmospheres is developed. These calculations are compared with microwave brightness temperatures at a wavelength of 1.55 cm measured on the Nimbus-5 satellite and rain rates derived from WSR-57 meteorological radar measurements. A specially designed ground based verification experiment was also performed wherein upward viewing microwave brightness temperature measurements at wavelengths of 1.55 cm and 0.81 cm were compared with directly measured rain rates.

Wilheit, T. T.↗

Evaluation of the Large-Scale and Regional Climatic Response Across North Africa to Natural Variability in Oceanic Modes and Terrestrial Vegetation Among the CMIP5 Models

Hydrologic variability is a serious threat to the poverty-stricken regions of North Africa. Meanwhile, the scientific community struggles to attribute these extreme climatic episodes to specific oceanic and land drivers. Prior modeling studies have not assessed simulated feedbacks over North Africa against an established observational benchmark, so their results are considered as untested, model-specific findings. The Coupled Model Intercomparison Project Phase Five (CMIP5) archive represents the state-of-the-art in climate projections, with most CMIP5 models now containing interactive vegetation phenology; however, there have been no studies to date of their simulated vegetation feedbacks. The capability of CMIP5 models at accurately simulating the forcing of both sea-surface temperature (SST) and regional leaf area index (LAI) anomalies on North African climate needs to be a key consideration in assessing the models’ overall credibility and determining appropriate model weighting for developing climate projections. The Generalized Equilibrium Feedback Assessment (GEFA) is a promising statistical method, which can be applied to either model output or observations, for isolating the local and remote impacts of individual oceanic or terrestrial forcings on regional climate. We performed a combined observational and modeling assessment of land-ocean-atmosphere interactions across the distinct ecological and moisture gradients of North Africa. We evaluated and demonstrated the reliability of the GEFA statistical method over North Africa using the Community Earth System Model (CESM), applied GEFA to observational data to quantify the observed forcing of ocean basin SST anomalies and regional LAI anomalies on North African climate, evaluated the CMIP5 models’ performance in terms of representing these key observed feedbacks, and formulated CMIP5 feedback performance metrics for weighting North African climate projections. Our study represented the first attempt to separate the observed roles of oceanic and vegetation feedbacks across North Africa, the first systematic assessment and intercomparison of land-ocean-atmosphere feedbacks in CMIP5, and the first exploration of vegetation feedbacks among CMIP5 models. The following work was accomplished. (1) The team generated 14 scientific publications and 30 presentations based on the DOE-funded research. (2) A stepwise version of the GEFA statistical method was developed in which unimportant forcings were dropped to increase the reliability of the results. (3) SGEFA was successfully validated through experiments with CESM in terms of its ability to isolate the atmospheric responses to individual oceanic or land forcings. (4) Observational evidence was revealed for the Sahel’s positive vegetation-rainfall feedback on the seasonal to interannual time scale, and it was attributed to a moisture recycling mechanism rather than an albedo mechanism. (5) An approach with SGEFA was developed in which the individual contributions of soil moisture versus vegetation forcings could be separated, revealing that the former forcing outweighed the latter for sub-Saharan Africa. (6) Tropical ocean temperatures were found to be key regulators of pan-tropical vegetation variability, especially for arid and semi-arid regions, including sub-Saharan Africa. (7) The CMIP5 models largely underestimated the importance of land feedbacks across the Sahel. (8) The general consensus among CMIP5 models indicates, for the late 21st century, a diminished seasonal predictability of sub-Saharan African regional climate and an elevated role of the land surface compared to oceanic drivers in regulating regional climate variability. (9) Land-ocean-atmosphere interactions were demonstrated to be key contributors to the seasonal predictability of African wildfire activity. (10) The El Djouf was determined to be of greater importance than the Bodélé depression in terms of providing trans-Atlantic dust transport to the Americas.

54 ENVIRONMENTAL SCIENCES↗

A Gradient Based Subgrid-Scale Parameterization for Ocean Mesoscale Eddies

Mesoscale eddies play an important role in transport of heat and biogeochemical tracers in the global ocean circulation. Resolving these energetic eddies, however, is challenging in ocean general circulation models (OGCM) because it requires a horizontal grid spacing of ≲1/8° that is computationally expensive. As a result, we are required to parameterize mesoscale eddy effects on large-scale ocean flows. In this work, we introduce a new subgrid-scale (SGS) model that is developed based on a Taylor series expansion of resolved variables to parameterize subgrid mesoscale eddy transports and momentum fluxes in OGCM. We have performed an a priori study to evaluate the performance of our new gradient model using high-resolution ocean simulations. Our results show that the gradient model well predicts the actual SGS thickness fluxes in the zonal and meridional directions in coarse-resolution simulations with the grid spacing ≳1/4°. The unresolved kinetic energy at the ocean surface is also skillfully estimated. More importantly, unlike current mesoscale eddy parameterizations, which are mainly developed based on an assumption of flat bottom ocean, our new SGS model can capture the structure of unresolved standing meanders at the ocean surface. We have also developed a dynamic procedure for setting in non-dimensional parameters in our new parameterization through a non-ad hoc and tuning-free method. Overall, this work suggests that implementing the gradient model in OGCM can improve the model accuracy with an affordable computational cost in eddy-permitting and non-eddying simulations.

54 ENVIRONMENTAL SCIENCES↗

Synthetic-aperture-radar imaging of the ocean surface using the slightly-rough facet model and a full surface-wave spectrum

A new model of synthetic-aperture-radar (SAR) imaging of ocean waves is described. The model is based on mapping individual, slightly-rough surface facets through the SAR processor into the image and responses of the facets in the image domain are added together coherently to give the composite image. A windowing technique allows both the orbital motion and the phase velocity of the long waves to be included. It is determined that the azimuthal cut-off is due to a smearing of the response of the facets in the image induced by the random orbital motion of the intermediate large-scale waves and that the focus adjustment that gives the greatest image contrast is half the phase velocity of the dominant long wave. The optimal processing technique, however, may consist of spatially offsetting the multiple looks on the image domain to compensate the propagation of long waves during the integration time of the SAR.

West, James C.↗

Simultaneous Aerosol and Ocean Polarimeter Products Using Coupled Atmosphere-Ocean Vector Radiative Transfer and Neural Networks: The PACE-MAPP Algorithm

We describe the PACE-MAPP algorithm that simultaneously retrieves aerosol and ocean optical parameters using multiangle and multi-channel polarimeter measurements from the SPEXone, Hyper-Angular Rainbow Polarimeter 2 (HARP2), and Ocean Color Instrument (OCI) instruments onboard the NASA Plankton, Aerosol, Cloud, ocean Ecosystem (PACE) observing system PACE-MAPP is adapted from the Research Scanning Polarimeter (RSP) Microphysical Aerosol Properties from Polarimetry (RSP-MAPP) algorithm. A key feature of the MAPP family of algorithms is the use of a coupled vector radiative transfer model such that the atmosphere and ocean are always considered together as one system. Consequently, conservation of energy ensures that negative water-leaving radiances do not occur. PACE-MAPP uses optimal estimation to simultaneously characterize the optical and microphysical properties of aerosol and ocean constituents, find the optimal solution, and reliably account for the uncertainties of each parameter. This coupled approach, together with multiangle, multi-channel polarimeter measurements, will enable retrievals of aerosol and water properties across the Earth’s oceans. The PACE-MAPP algorithm provides aerosol and ocean products for both the open ocean and coastal areas and is designed to be accurate, modular, and efficient by using fast neural networks that replace the time-consuming vector radiative transfer calculations. We provide an overview of the PACE-MAPP framework and also describe its modular components including its aerosol and hydrosol models, ocean bio-optical models, and thin cirrus model.

Snorre Stamnes↗

Impact into the earth's ocean floor - Preliminary experiments, a planetary model, and possibilities for detection

Impact processes and plate tectonics are invoked in an experimental study of craters larger than 100 km in diameter on the ocean floor. Although the results obtained from 22-caliber (383 m/sec) ammunition experiments using dense, saturated sand as a target medium cannot be directly scaled to large events, the phenomenology exhibited is that expected of actual craters on the ocean floor: steep, mixed ejecta plume, gravitational adjustment of the crater to form a shallow basin, and extensive reworking of the ejecta, rim, and floor materials by violent collapse of the transient water cavity. Excavation into the mantle is predicted, although asthenospheric influence on outer ring formation is not. The clearest geophysical signature of such a crater is not topography; detection should instead be based on gravity and geoid anomalies due to uplift of the Moho, magnetic anomalies, and seismic resolution of the Moho uplift and crater formation fault planes.

Mckinnon, W. B.↗

Limits on the Range and Rate of Change in Power Take-Off Load in Ocean Wave Energy Conversion: A Study Using Model Predictive Control

Previous work comparing power take-off (PTO) architectures for ocean wave-powered reverse osmosis suggests that variable displacement in the wave energy converter (WEC)-driven pump does not offer a significant performance advantage. A limitation of that study is that the WEC was subject to a constant load within a given sea state (“Coulomb damping”) and did not account for controlled, moment-to-moment variation of the PTO load enabled by a variable displacement pump. This study explores the potential performance advantage of a variable PTO load over Coulomb damping. Model predictive control is used to provide optimal load control with constraints on the PTO load. The constraints include minimum and maximum loads and a limit on the rate of load adjustment. Parameter studies on these constraints enable conclusions about PTO design requirements in addition to providing an estimated performance advantage over Coulomb damping. Numerical simulation of the Oyster 1 WEC is carried out with performance weighted by historical sea state data from Humboldt Bay, CA. The results show a performance advantage of up to 20% higher yearly-average power absorption over Coulomb damping. Additionally, the parameter studies suggest that the PTO load should be adjustable down to at least 25% of the maximum load and should be adjustable between the minimum and maximum loads within a few seconds.

16 TIDAL AND WAVE POWER↗

Estimation of Non-Seasonal Oceanic Mass Redistribution and its Geodynamic Effects from the Topex/Poseidon Observation

Mass redistribution in terms of angular momentum will change the Earth's rotation under the conservation of angular momentum. It also changes the external gravitational field according to Newton's gravitational law. These small geodynamic variations can be measured very accurately by modem space geodetic techniques. We report on new estimates of the oceanic tidal angular momentum and hence the tidal variations in the Earth's rotation rate and polar motion. These are deduced from the new tidal solutions GOT99 and TPXO.3, plus M2 and Mf estimates (with error bars) from TPXO.4a. Comparisons with observed diurnal/semidiurnal tidal variations are generally similar to previous model estimates. Comparisons of models and observations of Mf polar motion are confused, with neither models nor observations agreeing among themselves or with each other. Part of this is caused by differences in Mf currents; model TPXO.4 shows current-intensive Rossby modes along western boundaries of major basins, which are indeed expected on theoretical grounds. Unfortunately, there are no reliable in situ observations of Mf currents useful for comparing models. On non-tidal oceanic mass redistribution we compute the principal components of the multi-year T/P sea surface height (SSH) field using EOF (empirical orthogonal functions) taking into consideration of the latitude-dependent area weighting. Emphasis is placed upon non-seasonal and interannual oscillations on regional scales, including ENSO and the North Atlantic Oscillation. Steric corrections to the observed SSH due to sea surface temperature changes (which has no geodynamic effects) is now under study using POCM ocean model output.

Chao, Benjamin F.↗

Porting the WAVEWATCH III (v6.07) wave action source terms to GPU

Abstract. Surface gravity waves play a critical role in several processes, including mixing, coastal inundation, and surface fluxes. Despite the growing literature on the importance of ocean surface waves, wind–wave processes have traditionally been excluded from Earth system models (ESMs) due to the high computational costs of running spectral wave models. The development of the Next Generation Ocean Model for the DOE’s (Department of Energy) E3SM (Energy Exascale Earth System Model) Project partly focuses on the inclusion of a wave model, WAVEWATCH III (WW3), into E3SM. WW3, which was originally developed for operational wave forecasting, needs to be computationally less expensive before it can be integrated into ESMs. To accomplish this, we take advantage of heterogeneous architectures at DOE leadership computing facilities and the increasing computing power of general-purpose graphics processing units (GPUs). This paper identifies the wave action source terms, W3SRCEMD, as the most computationally intensive module in WW3 and then accelerates them via GPU. Our experiments on two computing platforms, Kodiak (P100 GPU and Intel(R) Xeon(R) central processing unit, CPU, E5-2695 v4) and Summit (V100 GPU and IBM POWER9 CPU) show respective average speedups of 2× and 4× when mapping one Message Passing Interface (MPI) per GPU. An average speedup of 1.4× was achieved using all 42 CPU cores and 6 GPUs on a Summit node (with 7 MPI ranks per GPU). However, the GPU speedup over the 42 CPU cores remains relatively unchanged (∼ 1.3×) even when using 4 MPI ranks per GPU (24 ranks in total) and 3 MPI ranks per GPU (18 ranks in total). This corresponds to a 35 %–40 % decrease in both simulation time and usage of resources. Due to too many local scalars and arrays in the W3SRCEMD subroutine and the huge WW3 memory requirement, GPU performance is currently limited by the data transfer bandwidth between the CPU and the GPU. Ideally, OpenACC routine directives could be used to further improve performance. However, W3SRCEMD would require significant code refactoring to make this possible. We also discuss how the trade-off between the occupancy, register, and latency affects the GPU performance of WW3.

58 GEOSCIENCES↗