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 145 records · Page 8

A semianalytic radiance model of ocean color

A semianalytical radiance model is developed which predicts the upwelled spectral radiance at the sea surface as a function of the phytoplankton pigment concentration for Morel Case 1 waters. The model is in good agreement with experimental measurements carried out in waters which were not included in the data base used to derive it. It suggests that the observed variability in the radiance is due to variations in the backscattering of plankton and the associated detrital material. The model is extended to include other material in the water, such as dissolved organic material, referred to as yellow substances, and detached coccoliths from coccolithophorids, e.g., Emiliana huxleyi. Potential applications include an improved biooptical algorithm for the retrieval of pigment concentrations from satellite imagery in the presence of interference from detached coccoliths and an improved atmospheric correction for satellite imagery. The model also serves to identify and to interpret deviations from Case 1 waters.

Gordon, Howard R.↗

Understanding the Transport of Patagonian Dust and Its Influence on Marine Biological Activity in the South Atlantic Ocean

Modeling and remote sensing techniques were applied to examine the horizontal and vertical transport pathways of Patagonian dust and quantify the effect of soluble-iron- laden mineral dust deposition on marine primary productivity in the South Atlantic Ocean (SAO) surface waters. The global chemistry transport model GEOS-Chem, implemented with an iron dissolution scheme, was applied to evaluate the atmospheric transport and deposition of mineral dust and bioavailable iron during two dust outbreaks originating in the source regions of Patagonia. In addition to this "rapidly released" iron, offline calculations were also carried out to estimate the amount of bioavailable iron leached during the residence time of dust in the ocean mixed layer. Model simulations showed that the horizontal and vertical transport pathways of Patagonian dust plumes were largely influenced by the synoptic meteorological patterns of high and low pressure systems. Model-predicted horizontal and vertical transport pathways of Patagonian dust over the SAO were in reasonable agreement with remotely-sensed data. Comparison between remotely-sensed and offline calculated ocean surface chlorophyll-a concentrations indicated that, for the two dust outbreaks examined in this study, the deposition of bioavailable iron in the SAO through atmospheric pathways was insignificant. As the two dust transport episodes examined here represent typical outflows of mineral dust from South American sources, our study suggests that the atmospheric deposition of mineral dust is unlikely to induce large scale marine primary productivity and carbon sequestration in the South Atlantic sector of the Southern Ocean.

Johnson, Matthew↗

A coupled dynamic-thermodynamic model of an ice-ocean system in the marginal ice zone

Thermodynamics are incorporated into a coupled ice-ocean model in order to investigate wind-driven ice-ocean processes in the marginal zone. Upswelling at the ice edge which is generated by the difference in the ice-air and air-water surface stresses is found to give rise to a strong entrainment by drawing the pycnocline closer to the surface. Entrainment is shown to be negligible outside the areas affected by the ice edge upswelling. If cooling at the top is included in the model, the heat and salt exchanges are further enhanced in the upswelling areas. It is noted that new ice formation occurs in the region not affected by ice edge upswelling, and it is suggested that the high-salinity mixed layer regions (with a scale of a few Rossby radii of deformation) will overturn due to cooling, possibly contributing to the formation of deep water.

Hakkinen, Sirpa↗

Assimilation of Sea Ice Thickness Derived from Cryosat-2 Along-Track Freeboard Measurements into the Met Office's Forecast Ocean Assimilation Model (FOAM)

The feasibility of assimilating sea ice thickness (SIT) observations derived from CryoSat-2 along-track measurements of sea ice freeboard is successfully demonstrated using a 3D-Var assimilation scheme, NEMOVAR, within the Met Office's global, coupled ocean–sea-ice model, Forecast Ocean Assimilation Model (FOAM). The CryoSat-2 Arctic freeboard measurements are produced by the Centre for Polar Observation and Modelling (CPOM) and are converted to SIT within FOAM using modelled snow depth. This is the first time along-track observations of SIT have been used in this way, with other centres assimilating gridded and temporally averaged observations. The assimilation leads to improvements in the SIT analysis and forecast fields generated by FOAM, particularly in the Canadian Arctic. Arctic-wide observation-minus-background assimilation statistics for 2015–2017 show improvements of 0.75 m mean difference and 0.41 m root-mean-square difference (RMSD) in the freeze-up period and 0.46 m mean difference and 0.33 m RMSD in the ice break-up period. Validation of the SIT analysis against independent springtime in situ SIT observations from NASA Operation IceBridge (OIB) shows improvement in the SIT analysis of 0.61 m mean difference (0.42 m RMSD) compared to a control without SIT assimilation. Similar improvements are seen in the FOAM 5 d SIT forecast. Validation of the SIT assimilation with independent Beaufort Gyre Exploration Project (BGEP) sea ice draft observations does not show an improvement, since the assimilated CryoSat-2 observations compare similarly to the model without assimilation in this region. Comparison with airborne electromagnetic induction (Air-EM) combined measurements of SIT and snow depth shows poorer results for the assimilation compared to the control, despite covering similar locations to the OIB and BGEP datasets. This may be evidence of sampling uncertainty in the matchups with the Air-EM validation dataset, owing to the limited number of observations available over the time period of interest. This may also be evidence of noise in the SIT analysis or uncertainties in the modelled snow depth, in the assimilated SIT observations, or in the data used for validation. The SIT analysis could be improved by upgrading the observation uncertainties used in the assimilation. Despite the lack of CryoSat-2 SIT observations available for assimilation over the summer due to the detrimental effect of melt ponds on retrievals, it is shown that the model is able to retain improvements to the SIT field throughout the summer months due to prior, wintertime SIT assimilation. This also results in regional improvements to the July modelled sea ice concentration (SIC) of 5 % RMSD in the European sector, due to slower melt of the thicker sea ice.

Emma K. Fiedler↗

Dynamic and Regression Modeling of Ocean Variability in the Tide-Gauge Record at Seasonal and Longer Periods

Comparison of monthly mean tide-gauge time series to corresponding model time series based on a static inverted barometer (IB) for pressure-driven fluctuations and a ocean general circulation model (OM) reveals that the combined model successfully reproduces seasonal and interannual changes in relative sea level at many stations. Removal of the OM and IB from the tide-gauge record produces residual time series with a mean global variance reduction of 53%. The OM is mis-scaled for certain regions, and 68% of the residual time series contain a significant seasonal variability after removal of the OM and IB from the tide-gauge data. Including OM admittance parameters and seasonal coefficients in a regression model for each station, with IB also removed, produces residual time series with mean global variance reduction of 71%. Examination of the regional improvement in variance caused by scaling the OM, including seasonal terms, or both, indicates weakness in the model at predicting sea-level variation for constricted ocean regions. The model is particularly effective at reproducing sea-level variation for stations in North America, Europe, and Japan. The RMS residual for many stations in these areas is 25-35 mm. The production of "cleaner" tide-gauge time series, with oceanographic variability removed, is important for future analysis of nonsecular and regionally differing sea-level variations. Understanding the ocean model's strengths and weaknesses will allow for future improvements of the model.

Hill, Emma M.↗

Turbulent Diffusivities for Momentum, Heat, Salt and Passive Scalars

A program began ten years ago to build a turbulence model to describe high Reynolds numbers flows. Specifically, the aim was to devise a model that would satisfy two basic conditions: reproduce available turbulence data (laboratory, numerical simulations, etc.) concerning shear driven flows, buoyancy driven flows, 2D turbulence, freely decaying turbulence etc., and yet be manageable so as to be used, for example, in General Circulation Models (GCM'S). The model was presented in a series of papers that have appeared in Physics of Fluids since 1996. A total of about 80 turbulence statistics were reproduced. The model has no adjustable parameters. The next step was to apply the model to construct the vertical diffusities K for moment K(sub m), temperature K(sub h), salt K(sub s). and passive scalars K(sub c). First, we took K(sub s) = K(sub h) and tested the model using the GFDL ocean model. The results for the profiles of T and S vs. depth are indistinguishable from those derived using the latest model, the KPP model by the NCAR group. Presently, we are running the same Geophysical Fluid Dynamics Laboratory (GFDL) code relaxing the assumptive K(sub s) = K(sub h). Indeed, the turbulence model yields a salt diffusivity that depends on Ri and R rho (= Beta delta S/delta z/alpha delta T/delta z) in such a way that K(sub s) may be quite different from K(sub h). Salt fingers and double diffusivity laboratory data are reproduced. Results from the ocean model will be available shortly. Finally, we are trying to derive the horizontal diffusivities with the goal of providing a physically acceptable representation of mesoscale eddies. The recently suggested GMW parameterization has improved several O-GCM results and the goal here is to try to derive/justify it from a turbulence model and/or propose improvements/modifications. Theoretical work is in progress.

Canuto, V. M.↗

Turbulent Diffusivities for Momentum, Heat, Salt and Passive Scalars

A program began ten years ago to build a turbulence model to describe high Reynolds numbers flows. Specifically, the aim was to devise a model that would satisfy two basic conditions: reproduce available turbulence data (laboratory, numerical simulations, etc.) concerning shear driven flows, buoyancy driven flows, 2D turbulence, freely decaying turbulence etc., and yet be manageable so as to be used, for example, in GCM'S. The model was presented in a series of papers that have appeared in Physics of Fluids since 1996. A total of about 80 turbulence statistics were reproduced. The model has no adjustable parameters. The next step was to apply the model to construct the vertical diffusities K for moment K(sub m), temperature K(sub h), salt K(sub s), and passive scalars K(sub c). First, we took K(sub c) = K(sub h) and tested the model using the GFDL ocean model. The results for the profiles of T and S vs. depth are indistinguishable from those derided using the latest model, the KPP model by the NCAR group. Presently, we are running the same GFDL code relaxing the assumptive K(sub c) = K(sub h). Indeed, the turbulence model yelds a salt diffusivity that depends on R(sub i) and R(sub rho) = ((beta)(delta)S/(alpha)(delta)T/(delta)z) in such a way that K(sub s) may be quite different from K(sub h). Salt fingers and double diffusivity laboratory data are reproduced. Results from the ocean model will be available shortly. Finally, we are trying to derive the horizontal diffusivities with the goal of providing a physically acceptable representation of mesoscale eddies. The recently suggested GMW parameterization has improved several O-GCM results and the goal here is to try to derive/justify it from a turbulence model and/or propose improvements/modifications. Theoretical work is in progress.

Canuto, V. M.↗

An implementation of a barotropic quasigeostrophic model of ocean circulation on the MPP

The implementation on the Massively Parallel Processor (MPP) of a barotropic quasigeostrophic model of ocean circulation is discussed. The mathematical model, including scalings and boundary conditions is discussed. The numerical scheme, which uses compact differencing is also discussed. The implementation of this model on the MPP is then presented. Finally, some performance results are given and compared to results obtained using the VPS-32 and one processor of a CRAY-2.

Grosch, C. E.↗

Data Assimilation Into a Coupled Ocean Atmosphere Model: Application to the 1997-1998 El Nino

As part of JPL's ocean data assimilation effort to study ocean circulation and seasonal-interannual climate variability, sea level anomaly observed by TOPEX altimeter, together with sea surface temperature and wind stress data, are assimilated into a simple coupled ocean atmosphere model of the tropical Pacific. Model-data consistency is examined. Impact of the assimilation (as initialization) on El Nino Southern Oscillation (ENSO) forecasts is evaluated. The coupled model consists of a shallow water component with two baroclinic modes, an Ekman shear layer, a simplified mixed-layer temperature equation, and a statistical atmosphere based on dominant correlations between historical surface temperature and wind stress anomaly data. The adjoins method is used to fit the coupled model to the data over various six-month periods from late 1996 to early 1998 by optimally adjusting the initial state, model parameters, and basis functions of the statistical atmosphere. On average, the coupled model can be fitted to the data to approximately within the data and representation errors (5 cm, 0.5 C, and 10 sq m/sq m for sea level, surface temperature, and pseudo wind stress anomalies, respectively). The estimated fields resemble observed spatio-temporal structure reasonably well. Hindcasts/forecasts of the 1997/1998 El Nino initialized from forced estimated ocean states and parameters are much more realistic than those simply initialized from ocean states (see figure below). In particular, the ability of the model to produce significant warming beyond the initial state is dramatically improved. Parameter estimation, which compensates for some model errors, is found to be important to obtaining better fits of the model to data and to improving forecasts.

Lee, Tong↗

A scattering model for ocean surface

A surface-scattering model based on an approximate solution of the integral equations for the surface tangential fields is developed for non-Gaussian distributed finitely conducting surfaces. It is a function of both the surface spectrum and the surface bispectrum which are defined to be the Fourier transforms of the surface-correlation function and the surface-skewness function, respectively. By employing a sea-surface spectrum proposed by Pierson (1964) and modified by Lee and Fung (1982) it is shown that good agreement is obtained between model and field measurements in polarization, incidence angle, and azimuth angle without adjusting surface parameters.

Chen, K. S.↗

Establishing an IERS Sub-Center for Ocean Angular Momentum

The primary responsibilities of the NCAR component of this project are the following: (1) Acting as liaison with the international ocean modeling community; and (2) Developing standardized algorithms to compute desired ocean model products, and providing template source code for these algorithms in widely used global ocean models.

Bryan, Frank↗

Ice-shelf freshwater triggers for the Filchner–Ronne Ice Shelf melt tipping point in a global ocean–sea-ice model

Abstract. Some ocean modeling studies have identified a potential tipping point from a low to a high basal melt regime beneath the Filchner–Ronne Ice Shelf (FRIS), Antarctica, with significant implications for subsequent Antarctic ice sheet mass loss. To date, investigation of the climate drivers and impacts of this possible event have been limited because ice-shelf cavities and ice-shelf melting are only now starting to be included in global climate models. Using a global ocean–sea-ice configuration of the Energy Exascale Earth System Model (E3SM) that represents both ocean circulations and melting within ice-shelf cavities, we explore freshwater triggers (iceberg melt and ice-shelf basal melt) of a transition to a high-melt regime at FRIS in a low-resolution (30 km in the Southern Ocean) global ocean–sea-ice model. We find that a realistic spatial distribution of iceberg melt fluxes is necessary to prevent the FRIS melt regime change from unrealistically occurring under historical-reanalysis-based atmospheric forcing. Further, improvement of the default parameterization for mesoscale eddy mixing significantly reduces a large regional fresh bias and weak Antarctic Slope Front structure, both of which precondition the model to melt regime change. Using two different stable model configurations, we explore the sensitivity of FRIS melt regime change to regional ice-sheet freshwater fluxes. Through a series of sensitivity experiments prescribing incrementally increasing melt rates from the smaller, neighboring ice shelves in the eastern Weddell Sea, we demonstrate the potential for an ice-shelf melt “domino effect” should the upstream ice shelves experience increased melt rates. The experiments also reveal that modest ice-shelf melt biases in a model, especially at coarse ocean resolution where narrow continental shelf dynamics are not well resolved, can lead to an unrealistic melt regime change at downstream ice shelves. Thus, we find that remote connections between melt fluxes at different ice shelves are sensitive to baseline model conditions. Our results highlight both the potential and the peril of simulating prognostic Antarctic ice-shelf melt rates in a low-resolution global model.

54 ENVIRONMENTAL SCIENCES↗

Global Earth Response to Loading by Ocean Tide Models

Mathematical and programming techniques to numerically calculate Earth response to global semidiurnal and diurnal ocean tide models were developed. Global vertical crustal deformations were evaluated for M sub 2, S sub 2, N sub 2, K sub 2, K sub 1, O sub 1, and P sub 1 ocean tide loading, while horizontal deformations were evaluated for the M sub 2 tidal load. Tidal gravity calculations were performed for M sub 2 tidal loads, and strain tensor elements were evaluated for M sub 2 loads. The M sub 2 solution used for the ocean tide included the effects of self-gravitation and crustal loading.

Estes, R. H.↗

The sensitivity of latent heat flux to the air humidity approximations used in ocean circulation models

In deriving the surface latent heat flux with the bulk formula for the thermal forcing of some ocean circulation models, two approximations are commonly made to bypass the use of atmospheric humidity in the formula. The first assumes a constant relative humidity, and the second supposes that the sea-air humidity difference varies linearly with the saturation humidity at sea surface temperature. Using climatological fields derived from the Marine Deck and long time series from ocean weather stations, the errors introduced by these two assumptions are examined. It is shown that the errors reach above 100 W/sq m over western boundary currents and 50 W/sq m over the tropical ocean. The two approximations also introduce erroneous seasonal and spatial variabilities with magnitudes over 50 percent of the observed variabilities.

Liu, W. Timothy↗

Quantifying the Influence of Size, Shape, and Density of Microplastics on Their Transport Modes: A Modeling Approach

Microplastics (MPs) pose significant risks to marine ecosystems and human health, necessitating accurate predictions of their distributions in aquatic environments for effective risk mitigation. However, understanding MP transport dynamics is challenging because of the inadequate representation of MP characteristics such as size, shape, and density in numerical models. Further, the accuracy of the MP vertical profiles in existing models has not been thoroughly validated. Thus, we developed an MP transport model within the Finite Volume Community Ocean Model framework (FVCOM-MP) by integrating MP characteristics. We validated FVCOM-MP against experimental and analytical data, focusing on various MP transport modes and transitions. FVCOM-MP successfully replicates MP profiles in different transport modes, including the bedload, surface-load, suspended-load, and mixed-load modes. Additionally, we introduce phase diagrams for classifying MP transport modes based on particle characteristics, enhancing our understanding of MP dynamics in aquatic systems. The transport modes for a number of real-world MP particles, including fishing line, plastic bag/bottle fragments, synthetic fibers, tire wear particles, polyvinyl chloride and expanded polystyrene pellets, were analyzed with these phase diagrams.

Microplastic transport, Settling velocity, Rising ↗

An Improved Heat Budget Estimation Including Bottom Effects for General Ocean Circulation Models

This paper studies the effects of the underwater light field on heat-budget calculations of general ocean circulation models for shallow waters. The presence of a bottom significantly alters the estimated heat budget in shallow waters, which affects the corresponding thermal stratification and hence modifies the circulation. Based on the data collected during the COBOP field experiment near the Bahamas, we have used a one-dimensional turbulence closure model to show the influence of the bottom reflection and absorption on the sea surface temperature field. The water depth has an almost one-to-one correlation with the temperature rise. Effects of varying the bottom albedo by replacing the sea grass bed with a coral sand bottom, also has an appreciable effect on the heat budget of the shallow regions. We believe that the differences in the heat budget for the shallow areas will have an influence on the local circulation processes and especially on the evaporative and long-wave heat losses for these areas. The ultimate effects on humidity and cloudiness of the region are expected to be significant as well.

Carder, Kendall↗