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 631 records · Page 35

A model of the general ocean circulation determined from a joint solution for the Earth's gravity field

If the geoid and the satellite position are known accurately, satellite altimetry can be used to determine the geostrophic velocity of the surface ocean currents. The purpose of this investigation is to simultaneously estimate the sea surface topography, zeta, the model for the gravity field, and the satellite orbit. Satellite tracking data from fourteen satellites were used; along with Seasat and Geosat altimeter data as well as surface gravity data for the solution. The estimated model of zeta compares well at long wavelengths with the hydrographic model of zeta. Covariance studies show that the geoid is separable from zeta up to degree 9, at which point geoid error becomes comparable to the signal of zeta.

Nerem, R. S.↗

Numerical simulation of the world ocean circulation

A multi-level model, based on the primitive equations, is developed for simulating the temperature and velocity fields produced in the world ocean by differential heating and surface wind stress. The model ocean has constant depth, free slip at the lower boundary, and neglects momentum advection; so that there is no energy exchange between the barotropic and baroclinic components of the motion, although the former influences the latter through temperature advection. The ocean model was designed to be coupled to the UCLA atmospheric general circulation model, for the study of the dynamics of climate and climate changes. But here, the model is tested by prescribing the observed seasonally varying surface wind stress and the incident solar radiation, the surface air temperature and humidity, cloudiness and the surface wind speed, which, together with the predicted ocean surface temperature, determine the surface flux of radiant energy, sensible heat and latent heat.

Takano, K.↗

Global Ocean Primary Production Trends in the Modern Ocean Color Satellite Record (1998-2015)

Ocean primary production (PP), representing the uptake of inorganic carbon through photosynthesis, supports marine life and affects carbon exchange with the atmosphere. It is difficult to ascertain its magnitude, variability, and trends due to our inability to measure it directly at large scales. Yet it is paramount for understanding changes in marine health, fisheries, and the global carbon cycle. Using assimilation of ocean color satellite data into an ocean biogeochemical model, we estimate that global net ocean PP has experienced a small but significant decline −0.8 PgC/y (−2.1%) per decade (P < 0.05) in the 18-year satellite record from 1998 to 2015. This decline is associated with shallowing surface mixed layer depth (−2.4% per decade) and decreasing nitrate concentrations (−3.2% per decade). Relative contributions to PP by various types of ocean phytoplankton have changed, with decreases in production by intermediate-sized phytoplankton represented by chlorophytes (−14.3% per decade). This is partially compensated by increases from the unique, more nutrient-efficient, coccolithophores (8.4% per decade). Geographically, the North and Equatorial Indian Oceans are responsible for much of the decline in PP, falling 0.16 and 0.69 PgC/y per decade, respectively. Reduced production by large, fast-growing diatoms along with chlorophytes characterizes the decline here. In contrast, increases in PP are found in the North and North Central Pacific. The increases here are led by chlorophytes in the North Pacific and the small cyanobacteria in the North Central Pacific. These results suggest that the multi-decadal satellite observational record, coupled with an underlying representation of marine biodiversity in a model, can monitor the uptake of carbon by phytoplankton and that changes, although small, are occurring in the global oceans.

Watson W Gregg↗

Polar Motion Constraints on Models of the Fortnightly Tide

Estimates of the near-fortnightly Mf ocean tide from Topex/Poseidon satellite altimetry and from numerical solutions to the shallow water equations agree reasonably well, at least in their basin-scale features. For example, both show that the Pacific Ocean tide lags the Atlantic tide by roughly 30 degrees. There are hints of finer scale agreements in the elevation fields, but noise levels are high. In contrast, estimates of Mf currents are only weakly constrained by the TP data, because high-wavenumber Rossby waves (with intense currents) are associated with relatively small perturbations in surface elevation. As a result, a wide range of Mf current fields are consistent with both the TP data and the hydrodynamic equations within a priori plausible misfit bounds. We find that a useful constraint on the Mf currents is provided by independent estimates of the Earth's polar motion. At the Mf period polar motion shows a weak signal (both prograde and retrograde) which must be almost entirely caused by the ocean tide. We have estimated this signal from the SPACE2000 time series, after applying a broad-band correction for atmospheric angular momentum. Although the polar motion estimates have relatively large uncertainties, they are sufficiently precise to fix optimum data weights in a global ocean inverse model of Mf. These weights control the tradeoff between fitting a prior hydrodynamic model of Mf and fitting the relatively noisy T/P measurements of Mf. The predicted polar motion from the final inverse model agrees remarkably well with the Mf polar motion observations. The preferred model is also consistent with noise levels suggested by island gauges, and it is marginally consistent with differences observed by subsetting the altimetry (to the small extent that this is possible). In turn, this new model of the Mf ocean tide allows the ocean component to be removed from Mf estimates of length of day, thus yielding estimates of complex Love numbers less contaminated by oceanic effects than has hitherto been possible.

Ray, Richard D.↗

A New Theoretical Framework for Parameterizing Nonequilibrium Fractionation During Evaporation From the Ocean

The evaporation model for water isotopes proposed by Craig and Gordon (1965, https://books.google.co.in/books?id=6wIKAQAAIAAJ) is used in most isotope-enabled atmospheric models for the parameterization of nonequilibrium fractionation during evaporation from the ocean. In this model, one of the most uncertain parameters is the nonequilibrium fractionation factor 𝛼 𝑘 . Many isotope models use the formulation of Merlivat and Jouzel (1979, https://doi.org/10.1029/jc084ic08p05029), which parameterizes 𝛼 𝑘 as a function of wind speed and distinguishes between a smooth and a rough regime to account for the effect of ocean waves. The resulting discontinuity in 𝛼 𝑘 between smooth and rough regimes has been disputed by several empirical studies. Here, we present a new approach to parameterizing 𝛼𝑘 by explicitly accounting for the influence of wave drag on the momentum flux near the surface. Following Cifuentes-Lorenzen et al. (2018, https://doi.org/10.1007/s10546-018-0376-0), we add a third wave-induced component to the total momentum flux, in addition to the viscous and turbulent components, and extend the definition of the eddy viscosity to account for the momentum flux due to waves and turbulent dissipation near the surface. The new scheme predicts a slight decrease of 𝛼 𝑘 with wind speed, similar to the smooth-regime parameterization of Merlivat and Jouzel (1979, https://doi.org/10.1029/jc084ic08p05029). This new parameterization is incorporated into the isotope-enabled Community Atmosphere Model, where it improves the correlation of simulated and measured vapor deuterium excess relative to the default version and a version with constant 𝛼 𝑘 , suggesting that it may be used as a valid representation of fractionation during evaporation from the ocean in future isotope models.

54 ENVIRONMENTAL SCIENCES↗

LANL Institutional Computing Close-out Report for Project t22_ocean_time_step

This year, my team made a very productive use of LANL IC time. Four papers were published that used IC resources, and two more are under review. These publications fall into three categories: 1.) Improving tide modeling in the global ocean by adding self attraction and loading (Barton et. al. 2022), ice shelf cavities (Pal et al. 2023), and local time stepping (Lilly et al. 2023). 2.) A new sea ice numerical formulation (Capodaglio 2023). 3.) Performance comparisons, methods, and test cases for ocean model development, This includes a verification suite for ocean models (Bishnu, submitted) and a comparison between Julia and Fortran (Strauss, submitted). A new time-stepping method was introduced in Calandrini et al. (2022). These publications, and the use of LANL IC resources, go hand-in-hand with our mentoring efforts to train young scientists. Two of the lead authors are graduate students who are conducting their PhD research: Kristin Barton at the University of Michigan and Jeremy Lilly at Oregon State. These are both their first publications, and they both have DOE funding and DOE mentors. In addition, Bishnu is a post-doctoral researcher at LANL; Strauss conducted his research as a senior in high school; and Capodaglio, Calandrini and Pal are all early-career scientists who were converted to staff in 2021 or 2022. Here we highlight two publications on improvements in tidal modeling: Barton et al. 2022 and Lilly et al. 2023.

58 GEOSCIENCES↗

Mid-latitude wind forced ocean circulation studies

A simple barotropic vorticity equation model was developed to study some of the various modeling factors that affect the characteristics of strong western boundary currents like the Gulf Stream and Kuroshio. Successful prediction of sea surface temperature, both in the climatological mean and over periods as short as 1 month requires that the heating tendency, due to horizontal advection of heat by these currents, be accurately modeled. Conventional, coarse resolution ocean models do not satisfactorily reproduce the dominant features of these currents. It is concluded that it is important to understand why they do not and what must be done to do so in the future.

Harrison, D. E.↗

Forecasting Ocean Chlorophyll in the Equatorial Pacific

Using a global ocean biogeochemical model combined with a forecast of physical oceanic and atmospheric variables from the NASA Global Modeling and Assimilation Office, we assess the skill of a chlorophyll concentrations forecast in the Equatorial Pacific for the period 2012-2015 with a focus on the forecast of the onset of the 2015 El Niño event. Using a series of retrospective 9-month hindcasts, we assess the uncertainties of the forecasted chlorophyll by comparing the monthly total chlorophyll concentration from the forecast with the corresponding monthly ocean chlorophyll data from the Suomi-National Polar-orbiting Partnership Visible Infrared Imaging Radiometer Suite (S-NPP VIIRS) satellite. The forecast was able to reproduce the phasing of the variability in chlorophyll concentration in the Equatorial Pacific, including the beginning of the 2015-2016 El Niño. The anomaly correlation coefficient (ACC) was significant (p less than 0.05) for forecast at 1-month (R=0.33), 8-month (R=0.42) and 9-month (R=0.41) lead times. The root mean square error (RMSE) increased from 0.0399 microgram ch1 L(exp -1) for the 1-month lead forecast to a maximum of 0.0472 microgram ch1 L(exp -1) for the 9-month lead forecast indicating that the forecast of the amplitude of chlorophyll concentration variability was getting worse. Forecasts with a 3-month lead time were on average the closest to the S-NPP VIIRS data (23% or 0.033 microgram ch1 L(exp -1)) while the forecast with a 9-month lead time were the furthest (31% or 0.042 microgram ch1 L(exp -1)). These results indicate the potential for forecasting chlorophyll concentration in this region but also highlights various deficiencies and suggestions for improvements to the current biogeochemical forecasting system. This system provides an initial basis for future applications including the effects of El Niño events on fisheries and other ocean resources given improvements identified in the analysis of these results.

Anomaly Correlation Coefficient (ACC)↗

Assessing the Skill of Chlorophyll Forecasts: Latest Development and Challenges Ahead Using the Case of the Equatorial Pacific

Using a global ocean biogeochemical model combined with a forecast of physical oceanic and atmospheric variables from the NASA Global Modeling and Assimilation Office, we assess the skill of a chlorophyll concentrations forecast in the Equatorial Pacific for the period 2012-2015 with a focus on the forecast of the onset of the 2015 El Nino event. Using a series of retrospective 9-month hindcasts, we assess the uncertainties of the forecasted chlorophyll by comparing the monthly total chlorophyll concentration from the forecast with the corresponding monthly ocean chlorophyll data from the Suomi-National Polar-orbiting Partnership Visible Infrared Imaging Radiometer Suite (S-NPP VIIRS) satellite. The forecast was able to reproduce the phasing of the variability in chlorophyll concentration in the Equatorial Pacific, including the beginning of the 2015-2016 El Nino. The anomaly correlation coefficient (ACC) was significant (p less than 0.05) for forecast at 1-month (R=0.33), 8-month (R=0.42) and 9-month (R=0.41) lead times. The root mean square error (RMSE) increased from 0.0399 microgram chl L(exp -1) for the 1-month lead forecast to a maximum of 0.0472 microgram chl L(exp -1) for the 9-month lead forecast indicating that the forecast of the amplitude of chlorophyll concentration variability was getting worse. Forecasts with a 3-month lead time were on average the closest to the S-NPP VIIRS data (23% or 0.033 microgram chl L(exp -1)) while the forecast with a 9-month lead time were the furthest (31% or 0.042 microgram chl L(exp -1)). These results indicate the potential for forecasting chlorophyll concentration in this region but also highlights various deficiencies and suggestions for improvements to the current biogeochemical forecasting system. This system provides an initial basis for future applications including the effects of El Nino events on fisheries and other ocean resources given improvements identified in the analysis of these results.

El Nin↗

The 26 December 2004 tsunami source estimated from satellite radar altimetry and seismic waves

The 26 December 2004 Indian Ocean tsunami was the first earthquake tsunami of its magnitude to occur since the advent of both digital seismometry and satellite radar altimetry. Both have independently recorded the event from different physical aspects. The seismic data has then been used to estimate the earthquake fault parameters, and a three-dimensional ocean-general-circulation-model (OGCM) coupled with the fault information has been used to simulate the satellite-observed tsunami waves. Here we show that these two datasets consistently provide the tsunami source using independent methodologies of seismic waveform inversion and ocean modeling. Cross-examining the two independent results confirms that the slip function is the most important condition controlling the tsunami strength, while the geometry and the rupture velocity of the tectonic plane determine the spatial patterns of the tsunami.

ocean general circulation model (OGCM)↗

The Potential for Predicting Precipitation on Seasonal-to-Interannual Timescales

The ability to predict precipitation several months in advance would have a significant impact on water resource management. This talk provides an overview of a project aimed at developing this prediction capability. NASA's Seasonal-to-Interannual Prediction Project (NSIPP) will generate seasonal-to-interannual sea surface temperature predictions through detailed ocean circulation modeling and will then translate these SST forecasts into forecasts of continental precipitation through the application of an atmospheric general circulation model and a "SVAT"-type land surface model. As part of the process, ocean variables (e.g., height) and land variables (e.g., soil moisture) will be updated regularly via data assimilation. The overview will include a discussion of the variability inherent in such a modeling system and will provide some quantitative estimates of the absolute upper limits of seasonal-to-interannual precipitation predictability.

Koster, R. D.↗

Methods of testing parameterizations: Vertical ocean mixing

The ocean's velocity field is characterized by an exceptional variety of scales. While the small-scale oceanic turbulence responsible for the vertical mixing in the ocean is of scales a few centimeters and smaller, the oceanic general circulation is characterized by horizontal scales of thousands of kilometers. In oceanic general circulation models that are typically run today, the vertical structure of the ocean is represented by a few tens of discrete grid points. Such models cannot explicitly model the small-scale mixing processes, and must, therefore, find ways to parameterize them in terms of the larger-scale fields. Finding a parameterization that is both reliable and plausible to use in ocean models is not a simple task. Vertical mixing in the ocean is the combined result of many complex processes, and, in fact, mixing is one of the less known and less understood aspects of the oceanic circulation. In present models of the oceanic circulation, the many complex processes responsible for vertical mixing are often parameterized in an oversimplified manner. Yet, finding an adequate parameterization of vertical ocean mixing is crucial to the successful application of ocean models to climate studies. The results of general circulation models for quantities that are of particular interest to climate studies, such as the meridional heat flux carried by the ocean, are quite sensitive to the strength of the vertical mixing. We try to examine the difficulties in choosing an appropriate vertical mixing parameterization, and the methods that are available for validating different parameterizations by comparing model results to oceanographic data. First, some of the physical processes responsible for vertically mixing the ocean are briefly mentioned, and some possible approaches to the parameterization of these processes in oceanographic general circulation models are described in the following section. We then discuss the role of the vertical mixing in the physics of the large-scale ocean circulation, and examine methods of validating mixing parameterizations using large-scale ocean models.

Tziperman, Eli↗

Visible-infrared remote-sensing model and applications for ocean waters

Remote sensing has become important in the ocean sciences, especially for research involving large spatial scales. To estimate the in-water constituents through remote sensing, whether carried out by satellite or airplane, the signal emitted from beneath the sea surface, the so called water-leaving radiance (L(w)), is of prime importance. The magnitude of L(w) depends on two terms: one is the intensity of the solar input, and the other is the reflectance of the in-water constituents. The ratio of the water-leaving radiance to the downwelling irradiance (E(d)) above the sear surface (remote-sensing reflectance, R(sub rs)) is independent of the intensity of the irradiance input, and is largely a function of the optical properties of the in-water constituents. In this work, a model is developed to interpret r(sub rs) for ocean water in the visible-infrared range. In addition to terms for the radiance scattered from molecules and particles, the model includes terms that describe contributions from bottom reflectance, fluorescence of gelbstoff or colored dissolved organic matter (CDOM), and water Raman scattering. By using this model, the measured R(sub rs) of waters from the West Florida Shelf to the Mississippi River plume, which covered a (concentration of chlorophyll a) range of 0.07 - 50 mg/cu m, were well interpreted. The average percentage difference (a.p.d.) between the measured and modeled R(sub rs) is 3.4%, and, for the shallow waters, the model-required water depth is within 10% of the chart depth. Simple mathematical simulations for the phytoplankton pigment absorption coefficient (a(sub theta)) are suggested for using the R(sub rs) model. The inverse problem of R(sub rs), which is to analytically derive the in-water constituents from R(sub rs) data alone, can be solved using the a(sub theta) functions without prior knowledge of the in-water optical properties. More importantly, this method avoids problems associated with a need for knowledge of the shape and value of the chlorophyll-specific absorption coefficient. The simulation was tested for a wide range of water types, including waters from Monterey Bay, the West Florida Shelf, and the Mississippi River plume. Using the simulation, the R(sub rs)-derived in-water absorption coefficients were consistent with the values from in-water measurements (r(exp 2) greater than 0.94, slope approximately 1.0). In the remote-sensing applications, a new approach is suggested for the estimation of primary production based on remote sensing. Using this approach, the calculated primary production (PP) values based upon remotely sensed data were very close to the measured values for the euphotic zone (r(exp 2) = 0.95, slope 1.26, and 32% average difference), while traditional, pigment-based PP model provided values only one-third the size of the measured data. This indicates a potential to significantly improve the accuracy of the estimation of primary production based upon remote sensing.

Lee, Zhongping↗

Improved Upper Ocean/Sea Ice Modeling in the GISS GCM for Investigating Climate Change

This project built on our previous results in which we highlighted the importance of sea ice in overall climate sensitivity by determining that for both warming and cooling climates, when sea ice was not allowed to change, climate sensitivity was reduced by 35-40%. We also modified the GISS 8 deg x lO deg atmospheric GCM to include an upper-ocean/sea-ice model involving the Semtner three-layer ice/snow thermodynamic model, the Price et al. (1986) ocean mixed layer model and a general upper ocean vertical advection/diffusion scheme for maintaining and fluxing properties across the pycnocline. This effort, in addition to improving the sea ice representation in the AGCM, revealed a number of sensitive components of the sea ice/ocean system. For example, the ability to flux heat through the ice/snow properly is critical in order to resolve the surface temperature properly, since small errors in this lead to unrestrained climate drift. The present project, summarized in this report, had as its objectives: (1) introducing a series of sea ice and ocean improvements aimed at overcoming remaining weaknesses in the GCM sea ice/ocean representation, and (2) performing a series of sensitivity experiments designed to evaluate the climate sensitivity of the revised model to both Antarctic and Arctic sea ice, determine the sensitivity of the climate response to initial ice distribution, and investigate the transient response to doubling CO2.

Source record↗

Unforced decadal fluctuations in a coupled model of the atmosphere and ocean mixed layer

Global average temperature in a 100-year control run of a model used for greenhouse gas response simulations showed low-frequency natural variability comparable in magnitude to that observed over the last 100 years. The model variability was found to be barotropic in the atmosphere, and located in the tropical strip with largest values near the equator in the Pacific. The model variations were traced to complex, low-frequency interactions between the meridional sea surface temperature gradients in the eastern equatorial Pacific, clouds at both high and low levels, and features of the tropical atmospheric circulation. The variations in these and other model parameters appear to oscillate between two limiting climate states. The physical scenario accounting for the oscillations on decadal time scales is almost certainly not found in the real world on shorter time scales due to limited resolution and the omission of key physics (e.g., equatorial ocean dynamics) in the model. The real message is that models with dynamical limitations can still produce significant long-term variability. Only a thorough physical diagnosis of such simulations and comparisons with decadal-length data sets will allow one to decide if faith in the model results is, or is not, warranted.

Barnett, T. P.↗

Interactive Visualization of High-Dimensional Petascale Ocean Data

We describe an application for interactive visualization of 5 petabytes of time-varying multivariate data from a high-resolution global ocean circulation model. The input data are 10311 hourly (ocean time) time steps of various 2D and 3D fields from a 22-billion point 1/48- degree “lat-lon cap” configuration of the MIT General Circulation Model (MITgcm). We map the global horizontal model domain onto our 128-screen (8x16) tiled display wall to produce a canonical tiling with approximately one MITgcm grid point per display pixel, and using this tiling we encode the entire time series for multiple native and computed scalar quantities at a collection of ocean depths. We reduce disk bandwidth requirements by converting the model’s floating point data to 16-bit fixed point values, and compressing those values with a lossless video encoder, which together allow synchronized playback at 24 time steps per second across all 128 displays. The application allows dynamic assignment of any two encoded tiles to any display, and has multiple interfaces for quickly specifying various orderly arrangements of tiles. All subsequent rendering is done on the fly, with run time control of colormaps, transfer functions, histogram equalization, and labeling. The two data streams on each screen can be rendered independently and combined in various ways, including blending, differencing, horizontal/ vertical wipes, and checkerboarding. The two data streams on any screen can optionally be displayed as a scatterplot in their joint attribute space. All scatterplots and map-view plots from the same x/y location and depth are linked so they all show the current brushable selection. Ocean scientists have used the system, and have found previously unidentified features in the data.

Interactive↗

Pushing the Limits of Aquatic Remote Sensing: Synthetic Data and Deep Learning for Fast Inverse Emulation of A Coupled Ocean-Atmosphere Radiative Transfer Model

The inversion of electromagnetic information to physical and biological properties of the water column is a notoriously difficult problem, yet fundamental to our ability of understanding aquatic processes on large time and space scales. There is now a growing necessity to develop pragmatic approaches that allow timely and effective extrapolation of local processes, to spatially resolved global products, and to promote operational and sustainable resource policy management. This presentation will discuss research integrating advanced biological and radiative modeling, high-end computation, and machine learning to develop a portable global processor for simultaneous retrieval of atmosphere and water optics for diverse aquatic systems from the open and coastal ocean to optically extreme inland waters and harmful algal blooms. We will discuss some of the basic concepts behind the forward modeling approach including DEAP, the novel Distributed Equivalent Algal Populations model, for developing large spectral libraries of aquatic particle optics to aid in our ability to distinguish phytoplankton functional types (PFTs) and inorganic material, as well as other factors which enable comprehensive modeling from the benthos to top-of-atmosphere (TOA). This information is being used to understand how we can leverage next-generation deep learning methods for maximum information retrieval and rapid image processing, while also providing capabilities to identify minimum sensor spectral requirements necessary for certain aquatic applications. Further, I will touch on how we envision this research to enable the aquatic community for science discovery and how we are moving closer towards the capability for high-fidelity global analysis of aquatic ecosystems.

Jeremy Alan Kravitz↗

Year-Round Pack Ice in the Weddell Sea, Antarctica: Response and Sensitivity to Atmospheric and Oceanic Forcing

Using a dynamic-thermodynamic numerical sea-ice model, external oceanic and atmospheric forcings on sea ice in the Weddell Sea are examined to identify physical processes associated with the seasonal cycle of pack ice, and to identify further the parameters that coupled models need to consider in predicting the response of the pack ice to climate and ocean-circulation changes. In agreement with earlier studies, the primary influence on the winter ice-edge maximum extent is air temperature. Ocean heat flux has more impact on the minimum-ice-edge extent and in reducing pack-ice thickness, especially in the eastern-Weddell Sea. Low relative humidity enhances ice growth in thin ice and open-water regions, producing a more realistic ice edge along the coastal areas of the western-Weddell Sea where dry continental air has an impact. The modeled extent of the Weddell summer pack is equally sensitive to ocean heat flux and atmospheric relative humidity variations with the more dynamic responses being from the atmosphere. Since the atmospheric regime in the eastern Weddell is dominated by marine intrusions from lower latitudes, with high humidity already, it is unlikely that either the moisture trans- port could be further raised or that it could be significantly lowered because of its distance from the continent (the lower humidity source). Ocean heat-transport variability is shown to lead to overall ice thinning in the model response and is a known feature of the actual system, as evidenced by the occurrence of the Weddell Polynya in the mid 1970s.

Geiger, Cathleen A.↗