Tsunami Squares: Earthquake driven inundation mapping and validation by comparison to the Regional Ocean Modeling System
Explore the source record for details and available documents.
David P. Grzan↗
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.
Skill Assessment of a Spectral Ocean-Atmosphere Radiative Model
Ocean phytoplankton, detrital material, and water absorb and scatter light spectrally. The Ocean- Atmosphere Spectral Irradiance Model (OASIM) is intended to provide surface irradiance over the oceans with sufficient spectral resolution to support ocean ecology, biogeochemistry, and heat exchange investigations, and of sufficient duration to support inter-annual and decadal investigations. OASIM total surface irradiance (integrated 200 nm to 4 microns) was compared to in situ data and three publicly available global data products at monthly 1-degree resolution. OASIM spectrally-integrated surface irradiance had root mean square (RMS) difference= 20.1 W/sq m (about 11%), bias=1.6 W/sq m (about 0.8%), regression slope= 1.01 and correlation coefficient= 0.89, when compared to 2322 in situ observations. OASIM had the lowest bias of any of the global data products evaluated (ISCCP-FD, NCEP, and ISLSCP 11), and the best slope (nearest to unity). It had the second best RMS, and the third best correlation coefficient. OASIM total surface irradiance compared well with ISCCP-FD (RMS= 20.7 W/sq m; bias=-11.4 W/sq m, r=0.98) and ISLSCP II (RMS =25.2 W/sq m; bias= -13.8 W/sq m; r=0.97), but less well with NCEP (RMS =43.0 W/sq m ;bias=-22.6 W/sq m; x=0.91). Comparisons of OASIM photosynthetically available radiation (PAR) with PAR derived from SeaWiFS showed low bias (-1.8 mol photons /sq m/d, or about 5%), RMS (4.25 mol photons /sq m/d ' or about 12%), near unity slope (1.03) and high correlation coefficient (0.97). Coupled with previous estimates of clear sky spectral irradiance in OASIM (6.6% RMS at 1 nm resolution), these results suggest that OASIM provides reasonable estimates of surface broadband and spectral irradiance in the oceans, and can support studies on ocean ecosystems, carbon cycling, and heat exchange.
Gregg, Watson, W.↗
A Modeling Study of Oceanic Response to Daily and Monthly Surface Forcing
The goal of this study is to investigate the effect of high-frequency surface forcing (wind stresses and heat fluxes) on upper-ocean response. We use the reduced-gravity quasi-isopycnal ocean model by Schopf and Loughe (1995) for this study. Two experiments are performed: one with daily and the other with monthly surface forcing. The two experiments are referred to as DD and MM, respectively. The daily surface wind stress is produced from the SSM/I wind data (Atlas et al. 1991) using the drag coefficient of Large and Pond (1982). The surface latent and sensible heat fluxes are estimated using the atmospheric mixed layer model by Seager et al. (1995) with the time-varying air temperature and specific humidity from the NCEP-NCAR reanalysis (Kalnay et al. 1996). The radiation is based on climatological shortwave radiation from the Earth Radiation Budget Experiment (ERBE) [Harrison et al. 1993] and the daily GEWEX SRB data. The ocean model domain is restricted to the Pacific Ocean with realistic land boundaries. At the southern boundary the model temperature and salinity are relaxed to the Levitus (1994) climatology. The time-mean SST distribution from MM is close to the observed SST climatology while the mean SST field from DD is about 1.5 C cooler. To identify the responsible processes, we examined the mean heat budgets and the heat balance during the first year (when the difference developed) in the two experiments. The analysis reveals that this is contributed by two factors. One is the difference in latent heat flux. The other is the difference in mixing processes. To further evaluate the responsible processes, we repeated the DD experiment by reducing the based vertical diffusion from 1e-4 to 0.5e-5. The resultant SST field becomes quite closer to the observed SST field. SST variability from the two experiments is generally similar, but the equatorial SST differences between the two experiments show interannual variations. We are investigating the possible mechanisms responsible for the different responses.
Sui, Chung-Hsiung↗
Technical Report Series on Global Modeling and Data Assimilation: A Coupled Ocean-Atmosphere Radiative Model for Global Ocean Biogeochemical Models - Volume 22
An ocean-atmosphere radiative model (OARM) evaluates irradiance availability and quality in the water column to support phytoplankton growth and drive ocean thermodynamics. An atmospheric component incorporates spectral and directional effects of clear and cloudy skies as a function of atmospheric optical constituents, and spectral reflectance across the air-sea interface. An oceanic component evaluates the propagation of spectral and directional irradiance through the water column as a function of water, five phytoplankton groups, and chromophoric dissolved organic matter. It tracks the direct and diffuse streams from the atmospheric component, and a third stream, upwelling diffuse irradiance. The atmospheric component of OARM was compared to data sources at the ocean surface with a coefficient of determination (r2) of 0.97 and a root mean square of 12.1%.
Gregg, Watson W.↗
Combining TOPEX and SeaWinds Data to Refine Models for Ocean Surface Emissivity: Wind Vector Signatures at 18, 21 and 37 GHz
TOPEX Microwave Radiometer (TMR) brightness temperature observations at 18, 21, and 37 GHz have been collocated with near-simultaneous SeaWinds wind vector data as well as with a monthly SST climatological product. The combined data set allows us to study the dependence of ocean surface emissivity (at each frequency) upon both wind speed. Results show clear two-branch wind speed dependence; weak and linear below 6 m/s with an abrupt increase in sensitivity above that point. The analysis also shows that the zenith-directed ocean surface emissivity is polarization dependent above wind speeds of 5-6 m/s with an azimuthal variation related to the wind direction. This last result accords with recent polarimetric radiometer observations collected from aircraft. Implications of these observations to water vapor retrieval algorithms, nadir-viewing polarimetry, and ocean emission modeling will be discussed.
Vandemark, D.↗
Inconsistent Strategies to Spin up Models in CMIP5: Implications for Ocean Biogeochemical Model Performance Assessment
During the fifth phase of the Coupled Model Intercomparison Project (CMIP5) substantial efforts were made to systematically assess the skill of Earth system models. One goal was to check how realistically representative marine biogeochemical tracer distributions could be reproduced by models. In routine assessments model historical hindcasts were compared with available modern biogeochemical observations. However, these assessments considered neither how close modeled biogeochemical reservoirs were to equilibrium nor the sensitivity of model performance to initial conditions or to the spin-up protocols. Here, we explore how the large diversity in spin-up protocols used for marine biogeochemistry in CMIP5 Earth system models (ESMs) contributes to model-to-model differences in the simulated fields. We take advantage of a 500-year spin-up simulation of IPSL-CM5A-LR to quantify the influence of the spin-up protocol on model ability to reproduce relevant data fields. Amplification of biases in selected biogeochemical fields (O2, NO3, Alk-DIC) is assessed as a function of spin-up duration. We demonstrate that a relationship between spin-up duration and assessment metrics emerges from our model results and holds when confronted with a larger ensemble of CMIP5 models. This shows that drift has implications for performance assessment in addition to possibly aliasing estimates of climate change impact. Our study suggests that differences in spin-up protocols could explain a substantial part of model disparities, constituting a source of model-to- model uncertainty. This requires more attention in future model intercomparison exercises in order to provide quantitatively more correct ESM results on marine biogeochemistry and carbon cycle feedbacks.
histories↗
The effect of ocean tides on the earth's rotation as predicted by the results of an ocean tide model
The published ocean tidal angular momentum results of Seiler (1991) are used to predict the effects of the most important semidiurnal, diurnal, and long period ocean tides on the earth's rotation. The separate, as well as combined, effects of ocean tidal currents and sea level height changes on the length-of-day, UT1, and polar motion are computed. The predicted polar motion results reported here account for the presence of the free core nutation and are given in terms of the motion of the celestial ephemeris pole so that they can be compared directly to the results of observations. Outside the retrograde diurnal tidal band, the summed effect of the semidiurnal and diurnal ocean tides studied here predict peak-to-peak polar motion amplitudes as large as 2 mas. Within the retrograde diurnal tidal band, the resonant enhancement caused by the free core nutation leads to predicted polar motion amplitudes as large as 9 mas.
Gross, Richard S.↗
Modeling Mesoscale Eddies in the North Atlantic Ocean
Ocean modeling plays an important role in understanding the current climatic conditions and predicting the future climate change. Modeling the ocean at eddy-permitting and/or eddy resolving resolutions (1/3 degree or higher) has a two-fold objective. One part is to represent the ocean as realistically as possible, because mesoscale eddies have an impact on the large-scale circulation. The second objective is to learn how to represent effects of mesoscale eddies without explicitly resolving them. This is particularly important for climate models which cannot be run at eddy-resolving resolutions because of the computational constraints. At JPL, a 1/6 degree latitude by 1/6 degree longitude with 37 vertical levels Atlantic Ocean model has been developed. The model is based on the Parallel Ocean Program (POP) developed at Los Alamos National Laboratory (LANL). Using the 256-processor Cray T3D, we have conducted a 40-year integration of this Atlantic eddy-resolving ocean model. A regional analysis demonstrate that many observed features associated with the Caribbean Sea eddies can be realistically simulated by this model. Analysis of this Atlantic eddy-resolving ocean model further suggests that these Caribbean Sea eddies are connected with eddies formed outside the Caribbean Sea at the confluence of the North Brazil Current (NBC) and the North Equatorial Countercurrent. The diagram of the model simulated surface current shows that the Caribbean eddies ultimately originate in the NBC retroflection region, traveling more than a year from the North Brazil coast through the Lesser Antilles into the Caribbean Sea and eventually into the Gulf of Mexico. Additional information is contained in the original.
Chao, Yi↗
The use of stellite scatterometer winds to drive a primitive equation model of the Indian Ocean: The impact of bandlike sampling
The aim of this study is to evaluate the impact of the bandlike sampling of spaceborne scatterometers on the ability of scatterometer winds to successfully force the mean flow and seasonal cycle of an ocean model in the context of equatorial and tropical dynamics. The equatorial ocean is simulated with a four-layer, primitive equation, reduced gravity model of the Indian Ocean. The variable wind stress used in this study is derived from one year (1988) of 6-hour analyses of the 10-m wind vector over the Indian Ocean performed at the European Centre for Medium-Range Weather Forecasts (ECMWF). It is applied as a forcing at every grid point of the model to drive a reference circulation. Scatterometer winds are simulated from ECMWF winds, using the nominal configurations and orbital parameters of the European Remote Sensing 1 (ERS-1) and NASA Scatterometer (NSCAT) missions. The model is forced in real time under swaths with the raw scatterometer winds of ERS-1 and NSCAT, with a persistence condition (i.e., the wind is kept constsnt until the next passage of the satellite provides a new value). The circulation obtained for each of the scatterometer experiments is compared with the reference circulation. The seasonal circulation of the Indian Ocean with NSCAT winds is very similar to the reference. The perturbations introduced by the bandlike sampling and the persistance condition have an impact similar to that of a small uncorrelated noise added to the reference forcing. The persistence condition for ERS-1 does not give results which are as good as those obtained for NSCAT.
Barnier, Bernard↗
Vacillations in a coupled ocean-atmosphere model
Results are presented from a 35-year integration of a coupled ocean-atmosphere model. Both ocean and atmosphere are two-level, nonlinear primitive equation models. The global atmospheric model is forced by a steady, zonally symmetric Newtonian heating. The ocean model is solved in a rectangular tropical basin. Heat fluxes between ocean and atmosphere are linear in air-sea temperature differences, and the interfacial stress is proportional to lower-level atmospheric winds. The coupled models produce ENSO-like variability on time scales of 3 to 5 years. Since there is no external time-dependent forcing, these are self-sustained vacillations of the nonlinear system. It is argued that the energetics of the vacillations is that of unstable coupled modes and that the time scale is crucially dependent on the effects of ocean waves propagating in a closed basin.
Schopf, Paul S.↗
An ice-ocean coupled model for the Northern Hemisphere
The Hibler ice model has been modified and adapted to a domain that includes most of the sea ice-covered areas in the Northern Hemisphere. This model, joined with the Cox ocean model, is developed as an enhancement to the U.S. Navy's sea ice forecasting, PIPS, and is termed PIPS2.0. Generally, the modeled ice edge is consistent with the Navy-NOAA Joint Ice Center weekly analysis, and the modeled ice thickness distribution agrees with submarine sonar data in the central Arctic basin.
Cheng, Abe↗
Consequences of Different Air-Sea Feedbacks on Ocean Using MITgcm and MERRA-2 Forcing: Implications for Coupled Data Assimilation Systems
Ocean surface flux estimates from atmospheric and oceanic reanalyses contain errors that compensate for inaccuracies in the respective atmosphere and ocean models used to generate these reanalyses. A conundrum for climate studies is the discrepancy between surface fluxes that minimize model-data differences for an atmosphere-only model vs surface fluxes that minimize model-data differences for an ocean model. As a first step towards a consistent coupled ocean-atmosphere data-assimilation (DA) system, we compare surface net heat flux from a state-of-the-art atmospheric reanalysis, the Modern-Era Retrospective analysis for Research and Applications, Version 2 (MERRA-2), to net heat flux from a state-of-the-art ocean state estimate, the Estimating the Circulation and Climate of the Ocean Version 4 (ECCO-v4). The possible impacts of the MERRA-2 and ECCO-v4 air-sea net heat flux difference in a coupled DA system were assessed using a set of experiments designed to imitate different “flavors” of a coupled DA system in an ocean-only setup. This was done by forcing the ECCO-v4 underlying ocean model - the Massachusetts Institute of Technology general circulation model (MITgcm) - with different sets of MERRA-2 fields and utilizing different forcing methods. By doing so we were able to turn off different air-sea feedbacks which, in a coupled DA setup, are partially muted by the constraining observations. The set of experiments, therefore, represents a range of active feedbacks in different “flavors” of coupled data-assimilation systems. For the period 1992–2011, MERRA-2 net heat flux has a global mean difference of -4.9 Wm(exp -2) relative to ECCO-v4. When MERRA-2 surface fields are used to force MITgcm, imbalances in the energy and the hydrological cycles of MERRA-2, which are directly related to the fact that MERRA-2 was created without an interactive ocean, propagate to the ocean. The experiment in which MITgcm is forced with MERRA-2 fluxes (MERRA-2-flux experiment) results in a 2.5°C global mean Sea Surface Temperature (SST) cooling, a 1m reduction in global mean sea level, and other drastic changes in the large scale ocean circulation relative to those resulting when the MITgcm is forced with the optimized ECCO-v4 net heat flux (the ECCO-v4 experiment itself). When MITgcm is forced with MERRA-2 state variables (MERRA-2-state experiment), the SST is somewhat restored to the observed SST, but the errors are shifted to the water cycle, resulting in a global mean sea level increase of 2.7 m. To further explore the pros and cons of these two approaches, we introduce a new intermediate forcing method in which the ocean is forced with turbulent fluxes but has a long wave feedback. This method, unlike MERRA-2 state, preserves the MERRA-2 water and salinity cycles, and it reduces the SST error compared to the MERRA-2-flux experiment, but the SST is not as good as that in the MERRA-2-state experiment. Our results have implications for ocean-model forcing recipes and clearly reveal the undesirable consequences of limiting the feedbacks in either these types of experiments or in coupled DA.
Ehud Strobach↗
Meridionally propagating interannual-to-interdecadal variability in a linear ocean-atmosphere model
Meridional oscillation modes in a global, primitive-equation coupled ocean-atmosphere model have been analyzed in order to determine whether they contain such meridionally propagating modes as surface-pressure perturbations with years-to-decades oscillation periods. A two-layer global ocean model and a two-level global atmosphere model were then formulated. For realistic parameter values and basic states, meridional modes oscillating at periods of several years to several decades are noted to be present in the coupled ocean-atmosphere model; the oscillation periods, travel times, and meridional structures of surface pressure perturbations in one of the modes are found to be comparable to the corresponding characteristics of observed sea-level pressure perturbations.
Mehta, Vikram M.↗
Comparison of MERRA-2 and ECCO-V4 Ocean Surface Heat Fluxes: Consequences of Different Forcing Feedbacks on Ocean Circulation and Implications for Climate Data Assimilation
Forcing ocean models with reanalysis data is a common practice in ocean modeling. As part of this practice, prescribed atmospheric state variables and interactive ocean SST (Sea Surface Temperature) are used to calculate fluxes between the ocean and the atmosphere. When forcing an ocean model with reanalysis fields, errors in the reanalysis data, errors in the ocean model and errors in the forcing formulation will generate a different solution compared to other ocean reanalysis solutions (which also have their own errors). As a first step towards a consistent coupled ocean-atmosphere reanalysis, we compare surface heat fluxes from a state-of-the-art atmospheric reanalysis, the Modern-Era Retrospective analysis for Research and Applications, Version 2 (MERRA-2), to heat fluxes from a state-of-the-art oceanic reanalysis, the Estimating the Circulation and Climate of the Ocean Version 4, Release 2 (ECCO-v4). Then, we investigate the errors associated with the MITgcm (Massachusetts Institute of Technology general circulation model) ocean model in its ECCO-v4 ocean reanalysis configuration (1992-2011) when it is forced with MERRA- 2 atmospheric reanalysis fields instead of with the ECCO-v4 adjoint optimized ERA-interim state variables. This is done by forcing ECCO-v4 ocean with and without feedbacks from MERRA-2 related to turbulent fluxes of heat and moisture and the outgoing long wave radiation. In addition, we introduce an intermediate forcing method that includes only the feedback from the interactive outgoing long wave radiation. The resulting ocean circulation is compared with ECCO-v4 reanalysis and in-situ observations. We show that, without feedbacks, imbalances in the energy and the hydrological cycles of MERRA-2 (which are directly related to the fact it was created without interactive ocean) result in considerable SST drifts and a large reduction in sea level. The bulk formulae and interactive outgoing long wave radiation, although providing air-sea feedbacks and reducing model-data misfit, strongly relax the ocean to observed SST and may result in unwanted features such as large change in the water budget. These features have implications in a desired forcing recipe to be used. The results strongly and unambiguously argue for next generation data assimilation climate studies to involve fully coupled systems.
MERRA-2↗
A Global Ocean Tide Model From TOPEX/POSEIDON Altimetry: GOT99.2
Goddard Ocean Tide model GOT99.2 is a new solution for the amplitudes and phases of the global oceanic tides, based on over six years of sea-surface height measurements by the TOPEX/POSEIDON satellite altimeter. Comparison with deep-ocean tide-gauge measurements show that this new tidal solution is an improvement over previous global models, with accuracies for the main semidiurnal lunar constituent M2 now below 1.5 cm (deep water only). The new solution benefits from use of prior hydrodynamic models, several in shallow and inland seas as well as the global finite-element model FES94.1. This report describes some of the data processing details involved in handling the altimetry, and it provides a comprehensive set of global cotidal charts of the resulting solutions. Various derived tidal charts are also provided, including tidal loading deformation charts, tidal gravimetric charts, and tidal current velocity (or transport) charts. Finally, low-degree spherical harmonic coefficients are computed by numerical quadrature and are tabulated for the major short-period tides; these are useful for a variety of geodetic and geophysical purposes, especially in combination with similar estimates from satellite laser ranging.
Ray, Richard D.↗
Global ocean tides through assimilation of oceanographic and altimeter satellite data in a hydrodynamic model
Ocean tides must be considered in many scientific disciplines: astronomy, oceanography, geodesy, geophysics, meteorology, and space technologies. Progress in each of these disciplines leads to the need for greater knowledge and more precise predictions of the ocean tide contribution. This is particularly true of satellite altimetry. On one side, the present and future satellite altimetry missions provide and will supply new data that will contribute to the improvement of the present ocean tide solutions. On the other side, tidal corrections included in the Geophysical Data Records must be determined with the maximum possible accuracy. The valuable results obtained with satellite altimeter data thus far have not been penalized by the insufficiencies of the present ocean tide predictions included in the geophysical data records (GDR's) because the oceanic processes investigated have shorter wavelengths than the error field of the tidal predictions, so that the residual errors of the tidal corrections are absorbed in the empirical tilt and bias corrections of the satellite orbit. For future applications to large-scale oceanic phenomena, however, it will no longer be possible to ignore these insufficiencies.
Leprovost, Christian↗
Intraseasonal and interannual oscillations in coupled ocean-atmosphere models
An investigation is presented of coupled ocean-atmosphere models' behavior in an environment where atmospheric wave speeds are substantially reduced from dry atmospheric values by such processes as condensation-moisture convergence. Modes are calculated for zonally periodic, unbounded ocean-atmosphere systems, emphasizing the importance of an inclusion of prognostic atmosphere equations in simple coupled ocean-atmosphere models with a view to simulations of intraseasonal variability and its possible interaction with interannual variability. The dynamics of low and high frequency modes are compared; both classes are sensitive to the degree to which surface wind anomalies are able to affect the evaporation rate.
Hirst, Anthony C.↗
A model of ocean basin crustal magnetization appropriate for satellite elevation anomalies
A model of ocean basin crustal magnetization measured at satellite altitudes is developed which will serve both as background to which anomalous magnetizations can be contrasted and as a beginning point for studies of tectonic modification of normal ocean crust. The model is based on published data concerned with the petrology and magnetization of the ocean crust and consists of viscous magnetization and induced magnetization estimated for individual crustal layers. Thermal remanent magnetization and chemical remanent magnetization are excluded from the model because seafloor spreading anomalies are too short in wavelength to be resolved at satellite altitudes. The exception to this generalization is found at the oceanic magnetic quiet zones where thermal remanent magnetization and chemical remanent magnetization must be considered along with viscous magnetization and induced magnetization.
Thomas, Herman H.↗