Downscaling Ocean Conditions: Initial Results using a Quasigeostrophic and Realistic Ocean Model

NASA Astrophysics Data System (ADS)

Katavouta, Anna; Thompson, Keith

2014-05-01

Previous theoretical work (Henshaw et al, 2003) has shown that the small-scale modes of variability of solutions of the unforced, incompressible Navier-Stokes equation, and Burgers' equation, can be reconstructed with surprisingly high accuracy from the time history of a few of the large-scale modes. Motivated by this theoretical work we first describe a straightforward method for assimilating information on the large scales in order to recover the small scale oceanic variability. The method is based on nudging in specific wavebands and frequencies and is similar to the so-called spectral nudging method that has been used successfully for atmospheric downscaling with limited area models (e.g. von Storch et al., 2000). The validity of the method is tested using a quasigestrophic model configured to simulate a double ocean gyre separated by an unstable mid-ocean jet. It is shown that important features of the ocean circulation including the position of the meandering mid-ocean jet and associated pinch-off eddies can indeed be recovered from the time history of a small number of large-scales modes. The benefit of assimilating additional time series of observations from a limited number of locations, that alone are too sparse to significantly improve the recovery of the small scales using traditional assimilation techniques, is also demonstrated using several twin experiments. The final part of the study outlines the application of the approach using a realistic high resolution (1/36 degree) model, based on the NEMO (Nucleus for European Modelling of the Ocean) modeling framework, configured for the Scotian Shelf of the east coast of Canada. The large scale conditions used in this application are obtained from the HYCOM (HYbrid Coordinate Ocean Model) + NCODA (Navy Coupled Ocean Data Assimilation) global 1/12 degree analysis product. Henshaw, W., Kreiss, H.-O., Ystrom, J., 2003. Numerical experiments on the interaction between the larger- and the small-scale motion of the Navier-Stokes equations. Multiscale Modeling and Simulation 1, 119-149. von Storch, H., Langenberg, H., Feser, F., 2000. A spectral nudging technique for dynamical downscaling purposes. Monthly Weather Review 128, 3664-3673.

The self-consistent dynamic pole tide in global oceans

NASA Technical Reports Server (NTRS)

Dickman, S. R.

1985-01-01

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

A higher-order split-step Fourier parabolic-equation sound propagation solution scheme.

PubMed

Lin, Ying-Tsong; Duda, Timothy F

2012-08-01

A three-dimensional Cartesian parabolic-equation model with a higher-order approximation to the square-root Helmholtz operator is presented for simulating underwater sound propagation in ocean waveguides. The higher-order approximation includes cross terms with the free-space square-root Helmholtz operator and the medium phase speed anomaly. It can be implemented with a split-step Fourier algorithm to solve for sound pressure in the model. Two idealized ocean waveguide examples are presented to demonstrate the performance of this numerical technique.

Elastic parabolic equation solutions for oceanic T-wave generation and propagation from deep seismic sources.

PubMed

Frank, Scott D; Collis, Jon M; Odom, Robert I

2015-06-01

Oceanic T-waves are earthquake signals that originate when elastic waves interact with the fluid-elastic interface at the ocean bottom and are converted to acoustic waves in the ocean. These waves propagate long distances in the Sound Fixing and Ranging (SOFAR) channel and tend to be the largest observed arrivals from seismic events. Thus, an understanding of their generation is important for event detection, localization, and source-type discrimination. Recently benchmarked seismic self-starting fields are used to generate elastic parabolic equation solutions that demonstrate generation and propagation of oceanic T-waves in range-dependent underwater acoustic environments. Both downward sloping and abyssal ocean range-dependent environments are considered, and results demonstrate conversion of elastic waves into water-borne oceanic T-waves. Examples demonstrating long-range broadband T-wave propagation in range-dependent environments are shown. These results confirm that elastic parabolic equation solutions are valuable for characterization of the relationships between T-wave propagation and variations in range-dependent bathymetry or elastic material parameters, as well as for modeling T-wave receptions at hydrophone arrays or coastal receiving stations.

Modeling and simulation of ocean wave propagation using lattice Boltzmann method

NASA Astrophysics Data System (ADS)

Nuraiman, Dian

2017-10-01

In this paper, we present on modeling and simulation of ocean wave propagation from the deep sea to the shoreline. This requires high computational cost for simulation with large domain. We propose to couple a 1D shallow water equations (SWE) model with a 2D incompressible Navier-Stokes equations (NSE) model in order to reduce the computational cost. The coupled model is solved using the lattice Boltzmann method (LBM) with the lattice Bhatnagar-Gross-Krook (BGK) scheme. Additionally, a special method is implemented to treat the complex behavior of free surface close to the shoreline. The result shows the coupled model can reduce computational cost significantly compared to the full NSE model.

DOE Office of Scientific and Technical Information (OSTI.GOV)

Saenz, Juan A.; Chen, Qingshan; Ringler, Todd

Recent work has shown that taking the thickness-weighted average (TWA) of the Boussinesq equations in buoyancy coordinates results in exact equations governing the prognostic residual mean flow where eddy–mean flow interactions appear in the horizontal momentum equations as the divergence of the Eliassen–Palm flux tensor (EPFT). It has been proposed that, given the mathematical tractability of the TWA equations, the physical interpretation of the EPFT, and its relation to potential vorticity fluxes, the TWA is an appropriate framework for modeling ocean circulation with parameterized eddies. The authors test the feasibility of this proposition and investigate the connections between the TWAmore » framework and the conventional framework used in models, where Eulerian mean flow prognostic variables are solved for. Using the TWA framework as a starting point, this study explores the well-known connections between vertical transfer of horizontal momentum by eddy form drag and eddy overturning by the bolus velocity, used by Greatbatch and Lamb and Gent and McWilliams to parameterize eddies. After implementing the TWA framework in an ocean general circulation model, we verify our analysis by comparing the flows in an idealized Southern Ocean configuration simulated using the TWA and conventional frameworks with the same mesoscale eddy parameterization.« less

The response of an ocean general circulation model to surface wind stress produced by an atmospheric general circulation model

DOE Office of Scientific and Technical Information (OSTI.GOV)

Huang, B.; Schneider, E.K.

1995-10-01

Two surface wind stress datasets for 1979-91, one based on observations and the other from an investigation of the COLA atmospheric general circulation model (AGCM) with prescribed SST, are used to drive the GFDL ocean general circulation model. These two runs are referred to as the control and COLA experiments, respectively. Simulated SST and upper-ocean heat contents (HC) in the tropical Pacific Ocean are compared with observations and between experiments. Both simulation reproduced the observed mean SST and HC fields as well as their annual cycles realistically. Major errors common to both runs are colder than observed SST in themore » eastern equatorial ocean and HC in the western Pacific south of the equator, with errors generally larger in the COLA experiment. New errors arising from the AGCM wind forcing include higher SST near the South American coast throughout the year and weaker HC gradients along the equator in boreal spring. The former is associated with suppressed coastal upwelling by weak along shore AGCM winds, and the latter is caused by weaker equatorial easterlies in boreal spring. The low-frequency ENSO fluctuations are also realistic for both runs. Correlations between the observed and simulated SST anomalies from the COLA simulation are as high as those from the control run in the central equatorial Pacific. A major problem in the COLA simulation is the appearance of unrealistic tropical cold anomalies during the boreal spring of mature El Nino years. These anomalies propagate along the equator from the western Pacific to the eastern coast in about three months, and temporarily eliminate the warm SST and HC anomalies in the eastern Pacific. This erroneous oceanic response in the COLA simulation is caused by a reversal of the westerly wind anomalies on the equator, associated with an unrealistic southward shift of the ITCZ in boreal spring during El Nino events. 66 refs., 16 figs.« less

The Influence of Ice-Ocean Interactions on Europa's Overturning Circulation

NASA Astrophysics Data System (ADS)

Zhu, P.; Manucharyan, G. E.; Thompson, A. F.; Goodman, J. C.; Vance, S.

2016-12-01

Jupiter's moon Europa appears to have a global liquid ocean, which is located beneath an ice shell that covers the moon's entire surface. Linking ocean dynamics and ice-ocean interactions is crucial to understanding observed surface features on Europa as well as other satellite measurements. Ocean properties and circulation may also provide clues as to whether the moon has the potential to support extraterrestrial life through chemical transport governed by ice-ocean interactions. Previous studies have identified a Hadley cell-like overturning circulation extending from the equator to mid latitudes. However, these model simulations do not consider ice-ocean interactions. In this study, our goal is to investigate how the ocean circulation may be affected by ice. We study two ice-related processes by building idealized models. One process is horizontal convection driven by an equator-to-pole buoyancy difference due to latitudinal ice transport at the ocean surface, which is found to be much weaker than the convective overturning circulation. The second process we consider is the freshwater layer formed by ice melting at the equator. A strong buoyancy contrast between the freshwater layer and the underlying water suppresses convection and turbulent mixing, which may modify the surface heat flux from the ocean to the bottom of the ice. We find that the salinity of the ocean below the freshwater layer tends to be homogeneous both vertically and horizontally with the presence of an overturning circulation. Critical values of circulation strength constrain the freshwater layer depth, and this relationship is sensitive to the average salinity of the ocean. Further coupling of temperature and salinity of the ice and the ocean that includes mutual influences between the surface heat flux and the freshwater layer may provide additional insights into the ice-ocean feedback, and its influence on the latitudinal difference of heat transport.

Bayesian Inference of High-Dimensional Dynamical Ocean Models

NASA Astrophysics Data System (ADS)

Lin, J.; Lermusiaux, P. F. J.; Lolla, S. V. T.; Gupta, A.; Haley, P. J., Jr.

2015-12-01

This presentation addresses a holistic set of challenges in high-dimension ocean Bayesian nonlinear estimation: i) predict the probability distribution functions (pdfs) of large nonlinear dynamical systems using stochastic partial differential equations (PDEs); ii) assimilate data using Bayes' law with these pdfs; iii) predict the future data that optimally reduce uncertainties; and (iv) rank the known and learn the new model formulations themselves. Overall, we allow the joint inference of the state, equations, geometry, boundary conditions and initial conditions of dynamical models. Examples are provided for time-dependent fluid and ocean flows, including cavity, double-gyre and Strait flows with jets and eddies. The Bayesian model inference, based on limited observations, is illustrated first by the estimation of obstacle shapes and positions in fluid flows. Next, the Bayesian inference of biogeochemical reaction equations and of their states and parameters is presented, illustrating how PDE-based machine learning can rigorously guide the selection and discovery of complex ecosystem models. Finally, the inference of multiscale bottom gravity current dynamics is illustrated, motivated in part by classic overflows and dense water formation sites and their relevance to climate monitoring and dynamics. This is joint work with our MSEAS group at MIT.

Dynamics of coupled ice-ocean system in the marginal ice zone: Study of the mesoscale processes and of constitutive equations for sea ice

NASA Technical Reports Server (NTRS)

Hakkinen, S.

1984-01-01

This study is aimed at the modelling of mesoscale processed such as up/downwelling and ice edge eddies in the marginal ice zones. A 2-dimensional coupled ice-ocean model is used for the study. The ice model is coupled to the reduced gravity ocean model (f-plane) through interfacial stresses. The constitutive equations of the sea ice are formulated on the basis of the Reiner-Rivlin theory. The internal ice stresses are important only at high ice concentrations (90-100%), otherwise the ice motion is essentially free drift, where the air-ice stress is balanced by the ice-water stress. The model was tested by studying the upwelling dynamics. Winds parallel to the ice edge with the ice on the right produce upwilling because the air-ice momentum flux is much greater that air-ocean momentum flux, and thus the Ekman transport is bigger under the ice than in the open water. The upwelling simulation was extended to include temporally varying forcing, which was chosen to vary sinusoidally with a 4 day period. This forcing resembles successive cyclone passings. In the model with a thin oceanic upper layer, ice bands were formed.

The pole tide in deep oceans

NASA Technical Reports Server (NTRS)

Dickman, S. R.

1990-01-01

The fluid-dynamical theory of the pole tide is examined by describing the oceanic response to the Chandler wobble and assessing its implications for mantle anelasticity and low-frequency ocean dynamics. The Laplace tide equations accounting for bottom friction are given, and a spherical harmonic approach is delineated in which the time-independent portion of the tide height is expanded. Pole-tide height and related inertia products are linearly proportional to wobble amplitude, and the final equations are modified to account for mantle elasticity and oceanic loading. Results for pole tide effects are given for various earth models with attention to the role of boundary constraints. A dynamic effect is identified which lengthens the Chandler period by about 1 day more than static lengthening, a contribution that suggests a vigorous low-frequency response. The values derived are shown to agree with previous models that do not incorporate the effects of the pole tide.

Sensitivity of equatorial Pacific and Indian Ocean watermasses to the position of the Indonesian Throughflow

NASA Astrophysics Data System (ADS)

Rodgers, Keith B.; Latif, Mojib; Legutke, Stephanie

2000-09-01

The sensitivity of the thermal structure of the equatorial Pacific and Indian Ocean pycnoclines to a model's representation of the Indonesian Straits connecting the two basins is investigated. Two integrations are performed using the global HOPE ocean model. The initial conditions and surface forcing for both cases are identical; the only difference between the runs is that one has an opening for the Indonesian Straits which spans the equator on the Pacific side, and the other has an opening which lies fully north of the equator. The resulting sensitivity throughout much of the upper ocean is greater than 0.5°C for both the equatorial Indian and Pacific. A realistic simulation of net Indonesian Throughflow (ITF) transport (measured in Sverdrups) is not sufficient for an adequate simulation of equatorial watermasses. The ITF must also contain a realistic admixture of northern and southern Pacific source water.

Intraseasonal and interannual oscillations in coupled ocean-atmosphere models

NASA Technical Reports Server (NTRS)

Hirst, Anthony C.; Lau, K.-M.

1990-01-01

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.

Ocean modelling on the CYBER 205 at GFDL

NASA Technical Reports Server (NTRS)

Cox, M.

1984-01-01

At the Geophysical Fluid Dynamics Laboratory, research is carried out for the purpose of understanding various aspects of climate, such as its variability, predictability, stability and sensitivity. The atmosphere and oceans are modelled mathematically and their phenomenology studied by computer simulation methods. The present state-of-the-art in the computer simulation of large scale oceans on the CYBER 205 is discussed. While atmospheric modelling differs in some aspects, the basic approach used is similar. The equations of the ocean model are presented along with a short description of the numerical techniques used to find their solution. Computational considerations and a typical solution are presented in section 4.

Eddy-resolving 1/10° model of the World Ocean

NASA Astrophysics Data System (ADS)

Ibrayev, R. A.; Khabeev, R. N.; Ushakov, K. V.

2012-02-01

The first results on simulating the intra-annual variability of the World Ocean circulation by use of the eddy-resolving model are considered. For this purpose, a model of the World Ocean with a 1/10° horizontal resolution and 49 vertical levels was developed (a 1/10 × 1/10 × 49 model of the World Ocean). This model is based on the traditional system of three-dimensional equations of the large-scale dynamics of the ocean and boundary conditions with an explicit allowance for water fluxes on the free surface of the ocean. The equations are written in the tripolar coordinate system. The numerical method is based on the separation of the barotropic and baroclinic components of the solution. Discretization in time is implemented using explicit schemes allowing effective parallelization for a large number of processors. The model uses the sub-models of the boundary layer of the atmosphere and the submodel of sea-ice thermodynamics. The model of the World Ocean was developed at the Institute of Numerical Mathematics of the Russian Academy of Sciences (INM RAS) and the P.P. Shirshov Institute of Oceanogy (IO RAS). The formulation of the problem of simulating the intra-annual variability of thermohydrodynamic processes of the World Ocean and the parameterizations that were used are considered. In the numerical experiment, the temporal evolution of the atmospheric effect is determined by the normal annual cycle according to the conditions of the international Coordinated Ocean-Ice Reference Experiment (CORE-I). The calculation was carried out on a multiprocessor computer with distributed memory; 1601 computational cores were used. The presented analysis demonstrates that the obtained results are quite satisfactory when compared to the results that were obtained by other eddy-resolving models of the global ocean. The analysis of the model solution is, to a larger extent, of a descriptive character. A detailed analysis of the results is to be presented in following works. This experiment is a significant first step in developing the eddy-resolving model of the World Ocean.

Using Laboratory Experiments to Improve Ice-Ocean Parameterizations

NASA Astrophysics Data System (ADS)

McConnochie, C. D.; Kerr, R. C.

2017-12-01

Numerical models of ice-ocean interactions are typically unable to resolve the transport of heat and salt to the ice face. Instead, models rely upon parameterizations that have not been sufficiently validated by observations. Recent laboratory experiments of ice-saltwater interactions allow us to test the standard parameterization of heat and salt transport to ice faces - the three-equation model. The three-equation model predicts that the melt rate is proportional to the fluid velocity while the experimental results typically show that the melt rate is independent of the fluid velocity. By considering an analysis of the boundary layer that forms next to a melting ice face, we suggest a resolution to this disagreement. We show that the three-equation model makes the implicit assumption that the thickness of the diffusive sublayer next to the ice is set by a shear instability. However, at low flow velocities, the sublayer is instead set by a convective instability. This distinction leads to a threshold velocity of approximately 4 cm/s at geophysically relevant conditions, above which the form of the parameterization should be valid. In contrast, at flow speeds below 4 cm/s, the three-equation model will underestimate the melt rate. By incorporating such a minimum velocity into the three-equation model, predictions made by numerical simulations could be easily improved.

A laboratory examination of the three-equation model of ice-ocean interactions

NASA Astrophysics Data System (ADS)

McConnochie, Craig; Kerr, Ross

2017-11-01

Numerical models of ice-ocean interactions are typically unable to resolve the transport of heat and salt to the ice face. As such, models rely upon parameterizations that have not been properly validated by data. Recent laboratory experiments of ice-saltwater interactions allow us to test the standard parameterization of heat and salt transport to ice faces - the `three equation model'. We find a significant disagreement in the dependence of the melt rate on the fluid velocity. The three-equation model predicts that the melt rate is proportional to the fluid velocity while the experimental results typically show that the melt rate is independent of the fluid velocity. By considering a theoretical analysis of the boundary layer next to a melting ice face we suggest a resolution to this disagreement. We show that the three-equation model assumes that the thickness of the diffusive sublayer is set by a shear instability. However, at low flow velocities, the sublayer is instead set by a convective instability. This distinction leads to a threshold velocity of approximately 4 cm/s at geophysically relevant conditions, above which the form of the parameterization should be valid. In contrast, at flow speeds below 4 cm/s, the three-equation model will underestimate the melt rate. ARC DP120102772.

Meridionally propagating interannual-to-interdecadal variability in a linear ocean-atmosphere model

NASA Technical Reports Server (NTRS)

Mehta, Vikram M.

1992-01-01

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.

The ocean quasi-homogeneous layer model and global cycle of carbon dioxide in system of atmosphere-ocean

NASA Astrophysics Data System (ADS)

Glushkov, Alexander; Glushkov, Alexander; Loboda, Nataliya; Khokhlov, Valery; Serbov, Nikoly; Svinarenko, Andrey

The purpose of this paper is carrying out the detailed model of the CO2 global turnover in system of "atmosphere-ocean" with using the ocean quasi-homogeneous layer model. Practically all carried out models are functioning in the average annual regime and accounting for the carbon distribution in bio-sphere in most general form (Glushkov et al, 2003). We construct a modified model for cycle of the carbon dioxide, which allows to reproduce a season dynamics of carbon turnover in ocean with account of zone ocean structure (up quasi-homogeneous layer, thermocline and deepest layer). It is taken into account dependence of the CO2 transfer through the bounder between atmosphere and ocean upon temperature of water and air, wind velocity, buffer mechanism of the CO2 dissolution. The same program is realized for atmosphere part of whole system. It is obtained a tempo-ral and space distribution for concentration of non-organic carbon in ocean, partial press of dissolute CO2 and value of exchange on the border between atmosphere and ocean. It is estimated a role of the wind intermixing of the up ocean layer. The increasing of this effect leads to increasing the plankton mass and further particles, which are transferred by wind, contribute to more quick immersion of microscopic shells and organic material. It is fulfilled investigation of sen-sibility of the master differential equations system solutions from the model parameters. The master differential equa-tions system, describing a dynamics of the CO2 cycle, is numerically integrated by the four order Runge-Cutt method under given initial values of valuables till output of solution on periodic regime. At first it is indicated on possible real-zation of the chaos scenario in system. On our data, the difference of the average annual values for the non-organic car-bon concentration in the up quasi-homogeneous layer between equator and extreme southern zone is 0.15 mol/m3, be-tween the equator and extreme northern zone is 0.12 mol/m3. the maximum amplitude of season oscillations (40° -50° n.l.) is 0.07 mol/m3. A link between global cycle of carbon dioxide and global climate change is investigated. Refrences: Glushkov A.V., Khokhlov V.N., Prepelitsa G.P., Tsenenko I.A., Optics of atmosphere and ocean.-2004.-Vol.14,N7.-p.219-223; Glushkov A.V., Loboda N.S., Khokhlov V.N., Atmospheric Research (Elseiver).-2005.-Vol.77.-P.100-113;Glushkov A.V., Loboda N.S., Khokhlov V.N., Lovett L. Journal of Hydrology (Elseiver).-2006.-Vol. 322. N1-4.-P.14-24; Glushkov A.V., Khokhlov V.N., Loboda N.S., Quart.J.Royal Meteorol. Soc.-2006.-Vol.132.- pp.447-465; Glushkov A.V., Khokhlov V.N., Loboda N.S., Ponomarenko E.L., Environm. Inf. Arch.-2003.-Vol.1.-P.125-130.

Indian Ocean Dipolelike Variability in the CSIRO Mark 3 Coupled Climate Model.

NASA Astrophysics Data System (ADS)

Cai, Wenju; Hendon, Harry H.; Meyers, Gary

2005-05-01

Coupled ocean-atmosphere variability in the tropical Indian Ocean is explored with a multicentury integration of the Commonwealth Scientific and Industrial Research Organisation (CSIRO) Mark 3 climate model, which runs without flux adjustment. Despite the presence of some common deficiencies in this type of coupled model, zonal dipolelike variability is produced. During July through November, the dominant mode of variability of sea surface temperature resembles the observed zonal dipole and has out-of-phase rainfall variations across the Indian Ocean basin, which are as large as those associated with the model El Niño-Southern Oscillation (ENSO). In the positive dipole phase, cold SST anomaly and suppressed rainfall south of the equator on the Sumatra-Java coast drives an anticyclonic circulation anomaly that is consistent with the steady response (Gill model) to a heat sink displaced south of the equator. The northwest-southeast tilting Sumatra-Java coast results in cold sea surface temperature (SST) centered south of the equator, which forces anticylonic winds that are southeasterly along the coast, which thus produces local upwelling, cool SSTs, and promotes more anticylonic winds; on the equator, the easterlies raise the thermocline to the east via upwelling Kelvin waves and deepen the off-equatorial thermocline to the west via off-equatorial downwelling Rossby waves. The model dipole mode exhibits little contemporaneous relationship with the model ENSO; however, this does not imply that it is independent of ENSO. The model dipole often (but not always) develops in the year following El Niño. It is triggered by an unrealistic transmission of the model's ENSO discharge phase through the Indonesian passages. In the model, the ENSO discharge Rossby waves arrive at the Sumatra-Java coast some 6 to 9 months after an El Niño peaks, causing the majority of model dipole events to peak in the year after an ENSO warm event. In the observed ENSO discharge, Rossby waves arrive at the Australian northwest coast. Thus the model Indian Ocean dipolelike variability is triggered by an unrealistic mechanism. The result highlights the importance of properly representing the transmission of Pacific Rossby waves and Indonesian throughflow in the complex topography of the Indonesian region in coupled climate models.

Simulations of the 2.5D inviscid primitive equations in a limited domain

NASA Astrophysics Data System (ADS)

Chen, Qingshan; Temam, Roger; Tribbia, Joseph J.

2008-12-01

The primitive equations (PEs) of the atmosphere and the oceans without viscosity are considered. These equations are not well-posed for any set of local boundary conditions. In space dimension 2.5 a set of nonlocal boundary conditions has been proposed in Chen et al. [Q. Chen, J. Laminie, A. Rousseau, R. Temam, J. Tribbia, A 2.5D Model for the equations of the ocean and the atmosphere, Anal. Appl. 5(3) (2007) 199-229]. The present article is aimed at testing the validity of these boundary conditions with physically relevant data. The issues tested are the well-posedness in the nonlinear case and the computational efficiency of the boundary conditions for limited area models [T.T. Warner, R.A. Peterson, R.E. Treadon, A tutorial on lateral boundary conditions as a basic and potentially serious limitation to regional numerical weather prediction, Bull. Amer. Meteor. Soc. 78(11) (1997) 2599-2617].

A low-order model of the equatorial ocean-atmosphere system

NASA Astrophysics Data System (ADS)

Legnani, Roberto

A low order model of the equatorial ocean-atmosphere coupled system is presented. The model atmosphere includes a hydrological cycle with cloud-radiation interaction. The model ocean is based on mixed layer dynamics with a parameterization of entrainment processes. The coupling takes place via transfer to momentum, sensible heat, latent heat and short wave and long wave radiation through the ocean surface. The dynamical formulation is that of the primitive equations of an equatorial beta-plane, with zonally periodic and meridionally infinite geometry. The system is expanded into the set of normal modes pertinent to the linear problem and severly truncated to a few modes; 54 degrees of freedom are retained. Some nonlinear terms of the equations are evaluated in physical space and then projected onto the functional space; other terms are evaluated directly in the functional space. Sensitivity tests to variations of the parameters are performed, and some results from 10-year initial value simulations are presented. The model is capable of supporting oscillations of different time scales, ranging from a few days to a few years; it prefers a particular zonally asymmetric state, but temporarily switches to a different (opposite) zonally asymmetric state in an event-like fashion.

a Low-Order Model of the Equatorial Ocean-Atmosphere System.

NASA Astrophysics Data System (ADS)

Legnani, Roberto

A low order model of the equatorial ocean-atmosphere coupled system is presented. The model atmosphere includes a hydrological cycle with cloud-radiation interaction. The model ocean is based on mixed layer dynamics with a parameterization of entrainment processes. The coupling takes place via transfer to momentum, sensible heat, latent heat and short -wave and long-wave radiation through the ocean surface. The dynamical formulation is that of the primitive equations of an equatorial beta-plane, with zonally periodic and meridionally infinite geometry. The system is expanded into the set of normal modes pertinent to the linear problem and severely truncated to a few modes; 54 degrees of freedom are retained. Some nonlinear terms of the equations are evaluated in physical space and then projected onto the functional space; other terms are evaluated directly in the functional space. Sensitivity tests to variations of the parameters are performed, and some results from 10-year initial value simulations are presented. The model is capable of supporting oscillations of different time scales, ranging from a few days to a few years; it prefers a particular zonally asymmetric state, but temporarily switches to a different (opposite) zonally asymmetric state in an event-like fashion.

Operational seasonal and interannual predictions of ocean conditions

NASA Technical Reports Server (NTRS)

Leetmaa, Ants

1992-01-01

Dr. Leetmaa described current work at the U.S. National Meteorological Center (NMC) on coupled systems leading to a seasonal prediction system. He described the way in which ocean thermal data is quality controlled and used in a four dimensional data assimilation system. This consists of a statistical interpolation scheme, a primitive equation ocean general circulation model, and the atmospheric fluxes that are required to force this. This whole process generated dynamically consist thermohaline and velocity fields for the ocean. Currently routine weekly analyses are performed for the Atlantic and Pacific oceans. These analyses are used for ocean climate diagnostics and as initial conditions for coupled forecast models. Specific examples of output products were shown both in the Pacific and the Atlantic Ocean.

Accurate polynomial expressions for the density and specific volume of seawater using the TEOS-10 standard

NASA Astrophysics Data System (ADS)

Roquet, F.; Madec, G.; McDougall, Trevor J.; Barker, Paul M.

2015-06-01

A new set of approximations to the standard TEOS-10 equation of state are presented. These follow a polynomial form, making it computationally efficient for use in numerical ocean models. Two versions are provided, the first being a fit of density for Boussinesq ocean models, and the second fitting specific volume which is more suitable for compressible models. Both versions are given as the sum of a vertical reference profile (6th-order polynomial) and an anomaly (52-term polynomial, cubic in pressure), with relative errors of ∼0.1% on the thermal expansion coefficients. A 75-term polynomial expression is also presented for computing specific volume, with a better accuracy than the existing TEOS-10 48-term rational approximation, especially regarding the sound speed, and it is suggested that this expression represents a valuable approximation of the TEOS-10 equation of state for hydrographic data analysis. In the last section, practical aspects about the implementation of TEOS-10 in ocean models are discussed.

The numerics of hydrostatic structured-grid coastal ocean models: State of the art and future perspectives

NASA Astrophysics Data System (ADS)

Klingbeil, Knut; Lemarié, Florian; Debreu, Laurent; Burchard, Hans

2018-05-01

The state of the art of the numerics of hydrostatic structured-grid coastal ocean models is reviewed here. First, some fundamental differences in the hydrodynamics of the coastal ocean, such as the large surface elevation variation compared to the mean water depth, are contrasted against large scale ocean dynamics. Then the hydrodynamic equations as they are used in coastal ocean models as well as in large scale ocean models are presented, including parameterisations for turbulent transports. As steps towards discretisation, coordinate transformations and spatial discretisations based on a finite-volume approach are discussed with focus on the specific requirements for coastal ocean models. As in large scale ocean models, splitting of internal and external modes is essential also for coastal ocean models, but specific care is needed when drying & flooding of intertidal flats is included. As one obvious characteristic of coastal ocean models, open boundaries occur and need to be treated in a way that correct model forcing from outside is transmitted to the model domain without reflecting waves from the inside. Here, also new developments in two-way nesting are presented. Single processes such as internal inertia-gravity waves, advection and turbulence closure models are discussed with focus on the coastal scales. Some overview on existing hydrostatic structured-grid coastal ocean models is given, including their extensions towards non-hydrostatic models. Finally, an outlook on future perspectives is made.

Parallel Computation of the Regional Ocean Modeling System (ROMS)

DOE Office of Scientific and Technical Information (OSTI.GOV)

Wang, P; Song, Y T; Chao, Y

2005-04-05

The Regional Ocean Modeling System (ROMS) is a regional ocean general circulation modeling system solving the free surface, hydrostatic, primitive equations over varying topography. It is free software distributed world-wide for studying both complex coastal ocean problems and the basin-to-global scale ocean circulation. The original ROMS code could only be run on shared-memory systems. With the increasing need to simulate larger model domains with finer resolutions and on a variety of computer platforms, there is a need in the ocean-modeling community to have a ROMS code that can be run on any parallel computer ranging from 10 to hundreds ofmore » processors. Recently, we have explored parallelization for ROMS using the MPI programming model. In this paper, an efficient parallelization strategy for such a large-scale scientific software package, based on an existing shared-memory computing model, is presented. In addition, scientific applications and data-performance issues on a couple of SGI systems, including Columbia, the world's third-fastest supercomputer, are discussed.« less

Response of an eddy-permitting ocean model to the assimilation of sparse in situ data

NASA Astrophysics Data System (ADS)

Li, Jian-Guo; Killworth, Peter D.; Smeed, David A.

2003-04-01

The response of an eddy-permitting ocean model to changes introduced by data assimilation is studied when the available in situ data are sparse in both space and time (typical for the majority of the ocean). Temperature and salinity (T&S) profiles from the WOCE upper ocean thermal data set were assimilated into a primitive equation ocean model over the North Atlantic, using a simple nudging scheme with a time window of about 2 days and a horizontal spatial radius of about 1°. When data are sparse the model returns to its unassimilated behavior, locally "forgetting" or rejecting the assimilation, on timescales determined by the local advection and diffusion. Increasing the spatial weighting radius effectively reduces both processes and hence lengthens the model restoring time (and with it, the impact of assimilation). Increasing the nudging factor enhances the assimilation effect but has little effect on the model restoring time.

GNuMe: A Galerkin-based Numerical Modeliing Environment for modeling geophysical fluid dynamics applications ranging from the Atmosphere to the Ocean

NASA Astrophysics Data System (ADS)

Giraldo, Francis; Abdi, Daniel; Kopera, Michal

2017-04-01

We have built a Galerkin-based Numerical Modeling Environment (GNuMe) for non hydrostatic atmospheric and ocean processes. GNuMe uses continuous Galerkin and Discontinuous Galerkin (CG/DG) discetizations as well as non-conforming adaptive mesh refinement (AMR), along with advanced time-integration methods that exploits both CG/DG and AMR capabilities. GNuMe currently solves the compressible and incompressible Navier-Stokes equations, the shallow water equations (with wetting and drying), and work is underway for inclusion of other types of equations. Moreover, GNuMe can run in both 2D and 3D modes on any type of accelerator hardware such as Nvidia GPUs and Intel KNL, and on standard X86 cores. In this talk, we shall present representative solutions obtained with GNuMe and will discuss where we think such a modeling framework could fit within standard Earth Systems Models. For further information on GNuMe please visit: http://frankgiraldo.wixsite.com/mysite/gnume.

Ocean tide models for satellite geodesy and Earth rotation

NASA Technical Reports Server (NTRS)

Dickman, Steven R.

1991-01-01

A theory is presented which predicts tides in turbulent, self-gravitating, and loading oceans possessing linearized bottom friction, realistic bathymetry, and continents (at coastal boundaries no-flow conditions are imposed). The theory is phrased in terms of spherical harmonics, which allows the tide equations to be reduced to linear matrix equations. This approach also allows an ocean-wide mass conservation constraint to be applied. Solutions were obtained for 32 long and short period luni-solar tidal constituents (and the pole tide), including the tidal velocities in addition to the tide height. Calibrating the intensity of bottom friction produces reasonable phase lags for all constituents; however, tidal amplitudes compare well with those from observation and other theories only for long-period constituents. In the most recent stage of grant research, traditional theory (Liouville equations) for determining the effects of angular momentum exchange on Earth's rotation were extended to encompass high-frequency excitations (such as short-period tides).

Ray Modeling Methods for Range Dependent Ocean Environments

DTIC Science & Technology

1983-12-01

the eikonal equation, gives rise to equations for ray paths which are perpendicular to the wave fronts. Equation II.4, the transport equation, leads... databases for use by MEDUSA. The author has assisted in the installation of MEDUSA at computer facilities which possess databases containing archives of...sound velocity profiles, bathymetry, and bottom loss data. At each computer site, programs convert the archival data retrieved by the database system

On the coupled evolution of oceanic internal waves and quasi-geostrophic flow

NASA Astrophysics Data System (ADS)

Wagner, Gregory LeClaire

Oceanic motion outside thin boundary layers is primarily a mixture of quasi-geostrophic flow and internal waves with either near-inertial frequencies or the frequency of the semidiurnal lunar tide. This dissertation seeks a deeper understanding of waves and flow through reduced models that isolate their nonlinear and coupled evolution from the Boussinesq equations. Three physical-space models are developed: an equation that describes quasi-geostrophic evolution in an arbitrary and prescribed field of hydrostatic internal waves; a three-component model that couples quasi-geostrophic flow to both near-inertial waves and the near-inertial second harmonic; and a model for the slow evolution of hydrostatic internal tides in quasi-geostrophic flow of near-arbitrary scale. This slow internal tide equation opens the path to a coupled model for the energetic interaction of quasi-geostrophic flow and oceanic internal tides. Four results emerge. First, the wave-averaged quasi-geostrophic equation reveals that finite-amplitude waves give rise to a mean flow that advects quasi-geostrophic potential vorticity. Second is the definition of a new material invariant: Available Potential Vorticity, or APV. APV isolates the part of Ertel potential vorticity available for balanced-flow evolution in Eulerian frames and proves necessary in the separating waves and quasi-geostrophic flow. The third result, hashed out for near-inertial waves and quasi-geostrophic flow, is that wave-flow interaction leads to energy exchange even under conditions of weak nonlinearity. For storm-forced oceanic near-inertial waves the interaction often energizes waves at the expense of flow. We call this extraction of balanced quasi-geostrophic energy 'stimulated generation' since it requires externally-forced rather than spontaneously-generated waves. The fourth result is that quasi-geostrophic flow can encourage or 'catalyze' a nonlinear interaction between a near-inertial wave field and its second harmonic that transfers energy to the small near-inertial vertical scales of wave breaking and mixing.

Horizontal Residual Mean Circulation: Evaluation of Spatial Correlations in Coarse Resolution Ocean Models

NASA Astrophysics Data System (ADS)

Li, Y.; McDougall, T. J.

2016-02-01

Coarse resolution ocean models lack knowledge of spatial correlations between variables on scales smaller than the grid scale. Some researchers have shown that these spatial correlations play a role in the poleward heat flux. In order to evaluate the poleward transport induced by the spatial correlations at a fixed horizontal position, an equation is obtained to calculate the approximate transport from velocity gradients. The equation involves two terms that can be added to the quasi-Stokes streamfunction (based on temporal correlations) to incorporate the contribution of spatial correlations. Moreover, these new terms do not need to be parameterized and is ready to be evaluated by using model data directly. In this study, data from a high resolution ocean model have been used to estimate the accuracy of this HRM approach for improving the horizontal property fluxes in coarse-resolution ocean models. A coarse grid is formed by sub-sampling and box-car averaging the fine grid scale. The transport calculated on the coarse grid is then compared to the transport on original high resolution grid scale accumulated over a corresponding number of grid boxes. The preliminary results have shown that the estimate on coarse resolution grids roughly match the corresponding transports on high resolution grids.

Numerical simulation of the world ocean circulation

NASA Technical Reports Server (NTRS)

Takano, K.; Mintz, Y.; Han, Y. J.

1973-01-01

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.

Open ocean tide modelling

NASA Technical Reports Server (NTRS)

Parke, M. E.

1978-01-01

Two trends evident in global tidal modelling since the first GEOP conference in 1972 are described. The first centers on the incorporation of terms for ocean loading and gravitational self attraction into Laplace's tidal equations. The second centers on a better understanding of the problem of near resonant modelling and the need for realistic maps of tidal elevation for use by geodesists and geophysicists. Although new models still show significant differences, especially in the South Atlantic, there are significant similarities in many of the world's oceans. This allows suggestions to be made for future locations for bottom pressure gauge measurements. Where available, estimates of M2 tidal dissipation from the new models are significantly lower than estimates from previous models.

The Predictability of Near-Coastal Currents Using a Baroclinic Unstructured Grid Model

DTIC Science & Technology

2011-12-28

clinic simulations. ADCIRC solves the time-dependent scalar transport equation for salinity and temperature. Through the equation of state...described by McDougall ct al. (2003), ADCIRC uses the temperature, salinity , and pressure in determining the density field. In order to avoid spurious...model. 2.3 Initialization and boundary forcing Temperature, salinity , elevation, and velocity fields from a regional ocean model are needed both to

A numerical world ocean general circulation model Part I. Basic design and barotropic experiment

NASA Astrophysics Data System (ADS)

Han, Young-June

1984-08-01

A new six-layer world ocean general circulation model based on the primitive system of equations is described in detail and its performance in the case of a homogeneous ocean is described. These test integrations show that the model is capable of reproducing the observed mean barotropic or vertically-integrated transport, as well as the seasonal variability of the major ocean gyres. The surface currents, however, are dominated by the Ekman transport, and such non-linear features as the western boundary currents and the equatorial countercurrents are poorly represented. The abyssal boundary countercurrents are also absent due to the lack of thermohaline forcing. The most conspicuous effect of the bottom topography on a homogeneous ocean is seen in the Southern ocean where the calculated Antarctic circumpolar transport through the Drake passage ( ≈ 10 Sv, with bathymetry included) greatly underestimates the observed transport (≈ 100 Sv).

Crustal control of dissipative ocean tides in Enceladus and other icy moons

NASA Astrophysics Data System (ADS)

Beuthe, Mikael

2016-12-01

Could tidal dissipation within Enceladus' subsurface ocean account for the observed heat flow? Earthlike models of dynamical tides give no definitive answer because they neglect the influence of the crust. I propose here the first model of dissipative tides in a subsurface ocean, by combining the Laplace Tidal Equations with the membrane approach. For the first time, it is possible to compute tidal dissipation rates within the crust, ocean, and mantle in one go. I show that oceanic dissipation is strongly reduced by the crustal constraint, and thus contributes little to Enceladus' present heat budget. Tidal resonances could have played a role in a forming or freezing ocean less than 100 m deep. The model is general: it applies to all icy satellites with a thin crust and a shallow ocean. Scaling rules relate the resonances and dissipation rate of a subsurface ocean to the ones of a surface ocean. If the ocean has low viscosity, the westward obliquity tide does not move the crust. Therefore, crustal dissipation due to dynamical obliquity tides can differ from the static prediction by up to a factor of two.

Subsurface Ocean Tides in Enceladus and Other Icy Moons

NASA Astrophysics Data System (ADS)

Beuthe, M.

2016-12-01

Could tidal dissipation within Enceladus' subsurface ocean account for the observed heat flow? Earthlike models of dynamical tides give no definitive answer because they neglect the influence of the crust. I propose here the first model of dissipative tides in a subsurface ocean, by combining the Laplace Tidal Equations with the membrane approach. For the first time, it is possible to compute tidal dissipation rates within the crust, ocean, and mantle in one go. I show that oceanic dissipation is strongly reduced by the crustal constraint, and thus contributes little to Enceladus' present heat budget. Tidal resonances could have played a role in a forming or freezing ocean less than 100 meters deep. The model is general: it applies to all icy satellites with a thin crust and a shallow or stratified ocean. Scaling rules relate the resonances and dissipation rate of a subsurface ocean to the ones of a surface ocean. If the ocean has low viscosity, the westward obliquity tide does not move the crust. Therefore, crustal dissipation due to dynamical obliquity tides can differ from the static prediction by up to a factor of two.

M2, S2, K1 models of the global ocean tide

NASA Technical Reports Server (NTRS)

Parke, M. E.; Hendershott, M. C.

1979-01-01

Ocean tidal signals appear in many geophysical measurements. Geophysicists need realistic tidal models to aid in interpretation of their data. Because of the closeness to resonance of dissipationless ocean tides, it is difficult for numerical models to correctly represent the actual open ocean tide. As an approximate solution to this problem, test functions derived by solving Laplace's Tidal Equations with ocean loading and self gravitation are used as a basis for least squares dynamic interpolation of coastal and island tidal data for the constituents M2, S2, and Kl. The resulting representations of the global tide are stable over at least a ?5% variation in the mean depth of the model basin, and they conserve mass. Maps of the geocentric tide, the induced free space potential, the induced vertical component of the solid earth tide, and the induced vertical component of the gravitational field for each contituent are presented.

Deriving Equations of State for Specific Lakes and Inland Seas from Laboratory Measurements

NASA Astrophysics Data System (ADS)

Andrulionis, Natalia; Zavialov, Ivan; Zavialov, Peter; Osadchiev, Alexander; Kolokolova, Alexandra; Alukaeva, Alevtina; Izhitskiy, Alexander; Izhitskaya, Elena

2017-04-01

The equation of state is the dependence of water density on temperature, salinity, and pressure. It is important in many respects, in particular, for numerical modeling of marine systems. The widely used UNESCO equation of state, as well as the more recent and general TEOS-10 equation, are intended for the ocean waters. Hence, they are confined to salinities below 40 ‰ and, even more restrictively, valid only for ionic salt composition characteristic for the ocean. Both conditions do not hold for many lakes. Moreover, significant deviations of the ionic composition from the oceanic one have been documented for coastal zones, especially those exposed to river discharges. Therefore, the objective of this study was to find equations of state for areas or water bodies with non-oceanic ionic salt composition. In order to obtain the required equations, we analyzed water samples obtained in expeditions of 2014-2016 from the Black Sea, the Aral Sea, Lake Issyk-Kul and Caspian Sea. The filtered samples were submitted to high accuracy (up to 0.00001 g/cm3) density measurements in laboratory using the Anton Paar DMA 5000M in the temperature range from 1 to 29°C. The absolute salinity values of the initial samples were obtained through the dry residue method. Further, we diluted the samples by purified deionized water to produce different salinities. To control the accuracy of the dilution process, we used a reference sample of standard IAPSO-certified seawater at 35‰. The density versus salinity and temperature data obtained thereby were then approximated by a best fitting 2-order polynomial surface using the least squares method. This procedure yielded the approximate empirical equations of state for the selected marine areas (the Russian Black Sea shelf) and inland water bodies (the Aral Sea, the Lake Issyk-Kul, the Caspian Sea). The newly derived equations - even the one for the Black Sea shelf - are different from the oceanic equation significantly within the confidence intervals. We also analyzed the salt content in all samples using the ionic chromotography method and the potentiometric titration method and discussed the relations between the ionic composition on the one hand and density on the other.

Evidence for self-refraction in a convergence zone: NPE (Nonlinear progressive wave equation) model results

NASA Technical Reports Server (NTRS)

Mcdonald, B. Edward; Plante, Daniel R.

1989-01-01

The nonlinear progressive wave equation (NPE) model was developed by the Naval Ocean Research and Development Activity during 1982 to 1987 to study nonlinear effects in long range oceanic propagation of finite amplitude acoustic waves, including weak shocks. The NPE model was applied to propagation of a generic shock wave (initial condition provided by Sandia Division 1533) in a few illustrative environments. The following consequences of nonlinearity are seen by comparing linear and nonlinear NPE results: (1) a decrease in shock strength versus range (a well-known result of entropy increases at the shock front); (2) an increase in the convergence zone range; and (3) a vertical meandering of the energy path about the corresponding linear ray path. Items (2) and (3) are manifestations of self-refraction.

Nonlinear and Dissipation Characteristics of Ocean Surface Waves in Estuarine Environments

DTIC Science & Technology

2013-09-30

developed models while using the general framework of operational wave models. We will conduct robustness tests of the system to determine the...and Guza (1984) model is weakly dispersive, in line with the assumptions behind the Boussinesq equations from which it was derived. The Kaihatu and...interactions across both frequency and directions. This system of equations is solved over a 2D frequency (f) and shore parallel wave number (κ) space. The

Global Paleobathymetry Reconstruction with Realistic Shelf-Slope and Sediment Wedge

NASA Astrophysics Data System (ADS)

Goswami, A.; Hinnov, L. A.; Gnanadesikan, A.; Olson, P.

2013-12-01

We present paleo-ocean bathymetry reconstructions in a 0.1°x0.1° resolution, using simple geophysical models (Plate Model Equation for oceanic lithosphere), published ages of the ocean floor (Müller et al. 2008), and modern world sediment thickness data (Divins 2003). The motivation is to create realistic paleobathymetry to understand the effect of ocean floor roughness on tides and heat transport in paleoclimate simulations. The values for the parameters in the Plate Model Equation are deduced from Crosby et al. (2006) and are used together with ocean floor age to model Depth to Basement. On top of the Depth to Basement, we added an isostatically adjusted multilayer sediment layer, as indicated from sediment thickness data of the modern oceans and marginal seas (Divins 2003). We also created another version of the sediment layer from the Müller et al. dataset. The Depth to Basement with the appropriate sediment layer together represent a realistic paleobathymetry. A Sediment Wedge was modeled to complement the reconstructed paleobathymetry by extending it to the coastlines. In this process we added a modeled Continental Shelf and Continental Slope to match the extent of the reconstructed paleobathymetry. The Sediment Wedge was prepared by studying the modern ocean where a complete history of seafloor spreading is preserved (north, south and central Atlantic Ocean, Southern Ocean between Australia-Antarctica, and the Pacific Ocean off the west coast of South America). The model takes into account the modern continental shelf-slope structure (as evident from ETOPO1/ETOPO5), tectonic margin type (active vs. passive margin) and age of the latest tectonic activity (USGS & CGMW). Once the complete ocean bathymetry is modeled, we combine it with PALEOMAP (Scotese, 2011) continental reconstructions to produce global paleoworld elevation-bathymetry maps. Modern time (00 Ma) was assumed as a test case. Using the above-described methodology we reconstructed modern ocean bathymetry, starting with age of the oceanic crust. We then reconstructed paleobathymetry for PETM (55 Ma) and Cenomanian-Turonian (90 Ma) times. For each case, the final products are: a) a global depth to basement measurement map based on plate model and EarthByte published age of the ocean crust for modern world; b) global oceanic crust bathymetry maps with a multilayer sediment layer (two versions with two types of sediment layers based on: i) observed total sediment thickness of the modern oceans and marginal seas, and ii) EarthByte-estimated global sediment data for 00 Ma); c) global oceanic bathymetry maps (two versions with two types of sediment layers) with reconstructed shelf and slope; and d) global elevation-bathymetry maps (two versions with two types of sediment layers) with continental elevations (PALEOMAP) and ocean bathymetry. Similar maps for other geological times can be produced using this method provided that ocean crustal age is known.

The Specific Features of Pollution Transport in the Northwest Pacific Ocean

NASA Astrophysics Data System (ADS)

Diansky, Nikolay; Fomin, Vladimir; Gusev, Anatoly

2013-04-01

Two calculations of pollutant dispersal in the Northwest Pacific Ocean are presented: (1) during possible shipwrecks in the process of spent nuclear fuel transportation from Petropavlovsk-Kamchatsky and (2) pollutant spread from the Japanese coast after the Fukushima 1 nuclear disaster on March 11, 2011. The circulation was simulated using a σ - coordinate ocean model INMOM (Institute of Numerical Mathematics Ocean Model) developed at the INM RAS. The INMOM is based on the primitive equations using the spherical σ - coordinate system with a free ocean surface. The INMOM was realized for the Pacific Ocean basin from the equator to the Bering Strait with a high 1/8° spatial resolution for reproducing the mesoscale ocean variability. The pollutant dispersal in the case of possible shipwrecks was estimated for currents for a statistically average year with atmospheric forcing from Common Ocean-ice Reference Experiments (CORE) for normal year data. The pollution spread from the Fukushima 1 nuclear power plant (NPP) was estimated for currents calculated with the real atmospheric forcing in accordance with the NCEP GFS (0.5 degree grid). The simulation period of pollutant dispersal from Fukushima 1 was 17 days: from March 11 to 28, 2011. The results of numerical simulation show that pollutant dispersal from the Fukushima 1 spread eastward according to the Kuroshio. Moreover, exceeding of natural background radiation level was simulated in the narrow region of the Japanese coast with width of less than 50 km.

``Supemodeling" by Coupling Multiple Atmospheres to A Single Ocean Simulates Single-ITCZ Climatology

NASA Astrophysics Data System (ADS)

Duane, G. S.; Shen, M. L.; Keenlyside, N. S.

2017-12-01

If the members of an ensemble of different models are allowed to interact with one another in run time, predictive skill can be improved as compared to that of any individual model or any average of individual model outputs. Inter-model connections in such an interactive ensemble can be trained, using historical data, so that the resulting ``supermodel" synchronizes with reality when observations are continuously assimilated, as in weather prediction. In climate-projection mode, the supermodel, after training, reproduces the attractor of the real system. We consider a variant of full supermodeling in which the models are only connected via the fluxes at the ocean-atmosphere interface. Two ECHAM atmospheres that differ in their convection parametrization schemes are thus connected to a single, shared ocean. The atmospheres partially synchronize at lower levels in the tropics, giving more realistic SST patterns than the constituent models: Although the constituent models both exhibit double ITCZ's, with cold tongues that extend too far west, the supermodel has the desired single ITCZ [Shen et al., Geophys. Res. Lett. 2016]. Here we explain the physical mechanism through which the supermodel removes even defects that are shared. One model (Nordeng) produces an unrealistically large zonal wind stress that results in upwelling of cold water and westward extension of the cold tongue. The other model (Tiedtke) produces an unrealistically small zonal wind stress that also implies a reduced wind stress curl off the equator because of Hadley-Walker coupling. The reduced wind stress curl leads to downwelling off the equator, and resultant upwelling of cold water at the equator through the tropical ocean cell. Thus the two constituent models give erroneous patterns of the same type, while the supermodel, which combines the models dynamically, avoids the error. If the models were linear, the errors of the two models would average; the success of the supermodel depends on nonlinearities in the east-west and north-south ocean-atmosphere feedbacks. It is argued that such behavior is widespread: supermodeling near-critical behavior in the coupled Earth System can give results that depend non-monotonically on the weights attached to the constituent models, thus surpassing those models, even when they err in the same way.

Ocean Modeling in an Eddying Regime

NASA Astrophysics Data System (ADS)

Hecht, Matthew W.; Hasumi, Hiroyasu

This monograph is the first to survey progress in realistic simulation in a strongly eddying regime made possible by recent increases in computational capability. Its contributors comprise the leading researchers in this important and constantly evolving field. Divided into three parts, • Oceanographic Processes and Regimes: Fundamental Questions • Ocean Dynamics and State: From Regional to Global Scale, and • Modeling at the Mesoscale: State of the Art and Future Directions the volume details important advances in physical oceanography based on eddy resolving ocean modeling. It captures the state of the art and discusses issues that ocean modelers must consider in order to effectively contribute to advancing current knowledge, from subtleties of the underlying fluid dynamical equations to meaningful comparison with oceanographic observations and leading-edge model development. It summarizes many of the important results which have emerged from ocean modeling in an eddying regime, for those interested broadly in the physical science. More technical topics are intended to address the concerns of those actively working in the field.

A four-dimensional primitive equation model for coupled coastal-deep ocean studies

NASA Technical Reports Server (NTRS)

Haidvogel, D. B.

1981-01-01

A prototype four dimensional continental shelf/deep ocean model is described. In its present form, the model incorporates the effects of finite amplitude topography, advective nonlinearities, and variable stratification and rotation. The model can be forced either directly by imposed atmospheric windstress and surface pressure distributions, and energetic mean currents imposed by the exterior oceanic circulation; or indirectly by initial distributions of shoreward propagation mesoscale waves and eddies. To avoid concerns over the appropriate specification of 'open' boundary conditions on the cross-shelf and seaward model boundaries, a periodic channel geometry (oriented along-coast) is used. The model employs a traditional finite difference expansion in the cross-shelf direction, and a Fourier (periodic) representation in the long-shelf coordinate.

Seasonal simulations using a coupled ocean-atmosphere model with data assimilation

NASA Astrophysics Data System (ADS)

Larow, Timothy Edward

1997-10-01

A coupled ocean-atmosphere initialization scheme using Newtonian relaxation has been developed for the Florida State University coupled ocean-atmosphere global general circulation model. The coupled model is used for seasonal predictions of the boreal summers of 1987 and 1988. The atmosphere model is a modified version of the Florida State University global spectral model, resolution triangular truncation 42 waves. The ocean general circulation model consists of a slightly modified version developed by Latif (1987). Coupling is synchronous with exchange of information every two model hours. Using daily analysis from ECMWF and observed monthly mean SSTs from NCEP, two - one year, time dependent, Newtonian relaxation were conducted using the coupled model prior to the seasonal forecasts. Relaxation was selectively applied to the atmospheric vorticity, divergence, temperature, and dew point depression equations, and to the ocean's surface temperature equation. The ocean's initial conditions are from a six year ocean-only simulation which used observed wind stresses and a relaxation towards observed SSTs for forcings. Coupled initialization was conducted from 1 June 1986 to 1 June 1987 for the 1987 boreal forecast and from 1 June 1987 to 1 June 1988 for the 1988 boreal forecast. Examination of annual means of net heat flux, freshwater flux and wind stress obtained by from the initialization show close agreement with Oberhuber (1988) climatology and the Florida State University pseudo wind stress analysis. Sensitivity of the initialization/assimilation scheme was tested by conducting two - ten member ensemble integrations. Each member was integrated for 90 days (June-August) of the respective year. Initial conditions for the ensembles consisted of the same ocean state as used by the initialize forecasts, while the atmospheric initial conditions were from ECMWF analysis centered on 1 June of the respective year. Root mean square error and anomaly correlations between observed and forecasted SSTs in the Nino 3 and Nino 4 regions show greater skill between the initialized forecasts than the ensemble forecasts. It is hypothesized that differences in the specific humidity within the planetary boundary layer are responsible for the large SST errors noted with the ensembles.

Mixing parametrizations for ocean climate modelling

NASA Astrophysics Data System (ADS)

Gusev, Anatoly; Moshonkin, Sergey; Diansky, Nikolay; Zalesny, Vladimir

2016-04-01

The algorithm is presented of splitting the total evolutionary equations for the turbulence kinetic energy (TKE) and turbulence dissipation frequency (TDF), which is used to parameterize the viscosity and diffusion coefficients in ocean circulation models. The turbulence model equations are split into the stages of transport-diffusion and generation-dissipation. For the generation-dissipation stage, the following schemes are implemented: the explicit-implicit numerical scheme, analytical solution and the asymptotic behavior of the analytical solutions. The experiments were performed with different mixing parameterizations for the modelling of Arctic and the Atlantic climate decadal variability with the eddy-permitting circulation model INMOM (Institute of Numerical Mathematics Ocean Model) using vertical grid refinement in the zone of fully developed turbulence. The proposed model with the split equations for turbulence characteristics is similar to the contemporary differential turbulence models, concerning the physical formulations. At the same time, its algorithm has high enough computational efficiency. Parameterizations with using the split turbulence model make it possible to obtain more adequate structure of temperature and salinity at decadal timescales, compared to the simpler Pacanowski-Philander (PP) turbulence parameterization. Parameterizations with using analytical solution or numerical scheme at the generation-dissipation step of the turbulence model leads to better representation of ocean climate than the faster parameterization using the asymptotic behavior of the analytical solution. At the same time, the computational efficiency left almost unchanged relative to the simple PP parameterization. Usage of PP parametrization in the circulation model leads to realistic simulation of density and circulation with violation of T,S-relationships. This error is majorly avoided with using the proposed parameterizations containing the split turbulence model. The high sensitivity of the eddy-permitting circulation model to the definition of mixing is revealed, which is associated with significant changes of density fields in the upper baroclinic ocean layer over the total considered area. For instance, usage of the turbulence parameterization instead of PP algorithm leads to increasing circulation velocity in the Gulf Stream and North Atlantic Current, as well as the subpolar cyclonic gyre in the North Atlantic and Beaufort Gyre in the Arctic basin are reproduced more realistically. Consideration of the Prandtl number as a function of the Richardson number significantly increases the modelling quality. The research was supported by the Russian Foundation for Basic Research (grant № 16-05-00534) and the Council on the Russian Federation President Grants (grant № MK-3241.2015.5)

Importance of ocean salinity for climate and habitability

PubMed Central

Cullum, Jodie; Stevens, David P.; Joshi, Manoj M.

2016-01-01

Modeling studies of terrestrial extrasolar planetary climates are now including the effects of ocean circulation due to a recognition of the importance of oceans for climate; indeed, the peak equator-pole ocean heat transport on Earth peaks at almost half that of the atmosphere. However, such studies have made the assumption that fundamental oceanic properties, such as salinity, temperature, and depth, are similar to Earth. This assumption results in Earth-like circulations: a meridional overturning with warm water moving poleward at the surface, being cooled, sinking at high latitudes, and traveling equatorward at depth. Here it is shown that an exoplanetary ocean with a different salinity can circulate in the opposite direction: an equatorward flow of polar water at the surface, sinking in the tropics, and filling the deep ocean with warm water. This alternative flow regime results in a dramatic warming in the polar regions, demonstrated here using both a conceptual model and an ocean general circulation model. These results highlight the importance of ocean salinity for exoplanetary climate and consequent habitability and the need for its consideration in future studies. PMID:27044090

Importance of ocean salinity for climate and habitability.

PubMed

Cullum, Jodie; Stevens, David P; Joshi, Manoj M

2016-04-19

Modeling studies of terrestrial extrasolar planetary climates are now including the effects of ocean circulation due to a recognition of the importance of oceans for climate; indeed, the peak equator-pole ocean heat transport on Earth peaks at almost half that of the atmosphere. However, such studies have made the assumption that fundamental oceanic properties, such as salinity, temperature, and depth, are similar to Earth. This assumption results in Earth-like circulations: a meridional overturning with warm water moving poleward at the surface, being cooled, sinking at high latitudes, and traveling equatorward at depth. Here it is shown that an exoplanetary ocean with a different salinity can circulate in the opposite direction: an equatorward flow of polar water at the surface, sinking in the tropics, and filling the deep ocean with warm water. This alternative flow regime results in a dramatic warming in the polar regions, demonstrated here using both a conceptual model and an ocean general circulation model. These results highlight the importance of ocean salinity for exoplanetary climate and consequent habitability and the need for its consideration in future studies.

Elastic parabolic equation solutions for underwater acoustic problems using seismic sources.

PubMed

Frank, Scott D; Odom, Robert I; Collis, Jon M

2013-03-01

Several problems of current interest involve elastic bottom range-dependent ocean environments with buried or earthquake-type sources, specifically oceanic T-wave propagation studies and interface wave related analyses. Additionally, observed deep shadow-zone arrivals are not predicted by ray theoretic methods, and attempts to model them with fluid-bottom parabolic equation solutions suggest that it may be necessary to account for elastic bottom interactions. In order to study energy conversion between elastic and acoustic waves, current elastic parabolic equation solutions must be modified to allow for seismic starting fields for underwater acoustic propagation environments. Two types of elastic self-starter are presented. An explosive-type source is implemented using a compressional self-starter and the resulting acoustic field is consistent with benchmark solutions. A shear wave self-starter is implemented and shown to generate transmission loss levels consistent with the explosive source. Source fields can be combined to generate starting fields for source types such as explosions, earthquakes, or pile driving. Examples demonstrate the use of source fields for shallow sources or deep ocean-bottom earthquake sources, where down slope conversion, a known T-wave generation mechanism, is modeled. Self-starters are interpreted in the context of the seismic moment tensor.

The use of stellite scatterometer winds to drive a primitive equation model of the Indian Ocean: The impact of bandlike sampling

NASA Technical Reports Server (NTRS)

Barnier, Bernard; Capella, Jorge; O'Brien, James J.

1994-01-01

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.

Dynamic ocean-tide effects on Earth's rotation

NASA Technical Reports Server (NTRS)

Dickman, S. R.

1993-01-01

This article develops 'broad-band' Liouville equations which are capable of determining the effects on the rotation of the Earth of a periodic excitation even at frequencies as high as semi-diurnal; these equations are then used to predict the rotational effects of altimetric, numerical and 32-constituent spherical harmonic ocean-tide models. The rotational model includes a frequency-dependent decoupled core, the effects of which are especially marked near retrograde diurnal frequencies; and a fully dynamic oceanic response, whose effects appear to be minor despite significant frequency dependence. The model also includes solid-earth effects which are frequency dependent as the result of both anelasticity at long periods and the fluid-core resonance at nearly diurnal periods. The effects of both tidal inertia and relative angular momentum on Earth rotation (polar motion, length of day, 'nutation' and Universal Time) are presented for 32 long- and short-period ocean tides determined as solutions to the author's spherical harmonic tide theory. The lengthening of the Chandler wobble period by the pole tide is also re-computed using the author's full theory. Additionally, using the spherical harmonic theory, tidal currents and their effects on rotation are determined for available numerical and altimetric tide height models. For all models, we find that the effects of tidal currents are at least as important as those of tide height for diurnal and semi-diurnal constituents.

Polar Motion Constraints on Models of the Fortnightly Tide

NASA Technical Reports Server (NTRS)

Ray, Richard D.; Egbert, G. D.; Smith, David E. (Technical Monitor)

2002-01-01

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.

Should tsunami models use a nonzero initial condition for horizontal velocity?

NASA Astrophysics Data System (ADS)

Nava, G.; Lotto, G. C.; Dunham, E. M.

2017-12-01

Tsunami propagation in the open ocean is most commonly modeled by solving the shallow water wave equations. These equations require two initial conditions: one on sea surface height and another on depth-averaged horizontal particle velocity or, equivalently, horizontal momentum. While most modelers assume that initial velocity is zero, Y.T. Song and collaborators have argued for nonzero initial velocity, claiming that horizontal displacement of a sloping seafloor imparts significant horizontal momentum to the ocean. They show examples in which this effect increases the resulting tsunami height by a factor of two or more relative to models in which initial velocity is zero. We test this claim with a "full-physics" integrated dynamic rupture and tsunami model that couples the elastic response of the Earth to the linearized acoustic-gravitational response of a compressible ocean with gravity; the model self-consistently accounts for seismic waves in the solid Earth, acoustic waves in the ocean, and tsunamis (with dispersion at short wavelengths). We run several full-physics simulations of subduction zone megathrust ruptures and tsunamis in geometries with a sloping seafloor, using both idealized structures and a more realistic Tohoku structure. Substantial horizontal momentum is imparted to the ocean, but almost all momentum is carried away in the form of ocean acoustic waves. We compare tsunami propagation in each full-physics simulation to that predicted by an equivalent shallow water wave simulation with varying assumptions regarding initial conditions. We find that the initial horizontal velocity conditions proposed by Song and collaborators consistently overestimate the tsunami amplitude and predict an inconsistent wave profile. Finally, we determine tsunami initial conditions that are rigorously consistent with our full-physics simulations by isolating the tsunami waves (from ocean acoustic and seismic waves) at some final time, and backpropagating the tsunami waves to their initial state by solving the adjoint problem. The resulting initial conditions have negligible horizontal velocity.

Effect of surface wave propagation in a four-layered oceanic crust model

NASA Astrophysics Data System (ADS)

Paul, Pasupati; Kundu, Santimoy; Mandal, Dinbandhu

2017-12-01

Dispersion of Rayleigh type surface wave propagation has been discussed in four-layered oceanic crust. It includes a sandy layer over a crystalline elastic half-space and over it there are two more layers—on the top inhomogeneous liquid layer and under it a liquid-saturated porous layer. Frequency equation is obtained in the form of determinant. The effects of the width of different layers as well as the inhomogeneity of liquid layer, sandiness of sandy layer on surface waves are depicted and shown graphically by considering all possible case of the particular model. Some special cases have been deduced, few special cases give the dispersion equation of Scholte wave and Stoneley wave, some of which have already been discussed elsewhere.

Simulation of Wave-Current Interaction Using a Three-Dimensional Hydrodynamic Model Coupled With a Phase Averaged Wave Model

NASA Astrophysics Data System (ADS)

Marsooli, R.; Orton, P. M.; Georgas, N.; Blumberg, A. F.

2016-02-01

The Stevens Institute of Technology Estuarine and Coastal Ocean Model (sECOM) has been coupled with a more advanced surface wave model to simulate wave‒current interaction, and results have been validated in estuarine and nearshore waters. sECOM is a three‒dimensional, hydrostatic, free surface, primitive equation model. It solves the Navier‒Stokes equations and the conservation equations for temperature and salinity using a finite‒difference method on an Arakawa C‒grid with a terrain‒following (sigma) vertical coordinate and orthogonal curvilinear horizontal coordinate system. The model is coupled with the surface wave model developed by Mellor et al. (2008), which solves the spectral equation and takes into account depth and current refraction, and deep and shallow water. The wave model parameterizes the energy distribution in frequency space and the wave‒wave interaction process by using a specified spectrum shape. The coupled wave‒hydrodynamic model considers the wave‒current interaction through wave‒induced bottom stress, depth‒dependent radiation stress, and wave effects on wind‒induced surface stress. The model is validated using the data collected at a natural sandy beach at Duck, North Carolina, during the DUCK94 experiment. This test case reveals the capability of the model to simulate the wave‒current interaction in nearshore coastal systems. The model is further validated using the data collected in Jamaica Bay, a semi‒enclosed body of water located in New York City region. This test reveals the applicability of the model to estuarine systems. These validations of the model and comparisons to its prior wave model, the Great Lakes Environmental Research Laboratory (GLERL) wave model (Donelan 1977), are presented and discussed. ReferencesG.L. Mellor, M.A. Donelan, and L‒Y. Oey, 2008, A Surface Wave Model for Coupling with Numerical Ocean Circulation Models. J. Atmos. Oceanic Technol., 25, 1785‒1807.Donelan, M. A 1977. A simple numerical model for wave and wind stress application. Report, National Water Research Institute, Burlington, Ontario, Canada, 28 pp.

Modeling the Gulf Stream System: How Far from Reality?

NASA Technical Reports Server (NTRS)

Choa, Yi; Gangopadhyay, Avijit; Bryan, Frank O.; Holland, William R.

1996-01-01

Analyses of a primitive equation ocean model simulation of the Atlantic Ocean circulation at 1/6 deg horizontal resolution are presented with a focus on the Gulf Stream region. Among many successful features of this simulation, this letter describes the Gulf Stream separation from the coast of North America near Cape Hatteras, meandering of the Gulf Stream between Cape Hatteras and the Grand Banks, and the vertical structure of temperature and velocity associated with the Gulf Stream. These results demonstrate significant improvement in modeling the Gulf Stream system using basin- to global scale ocean general circulation models. Possible reasons responsible for the realistic Gulf Stream simulation are discussed, contrasting the major differences between the present model configuration and those of previous eddy resolving studies.

Global Paleobathymetry for the Cenomanian-Turonian (90 Ma)

NASA Astrophysics Data System (ADS)

Goswami, A.; Olson, P.; Hinnov, L. A.; Gnanadesikan, A.

2014-12-01

We present a paleo-ocean bathymetry reconstruction for Cenomanian-Turonian (90 Ma) time in a 0.1°x0.1° resolution for use in paleo-climate studies. Age of the ocean floor for the Cenomanian-Turonian (90 Ma) is from Müller et al. (2008 a,b); coastlines are from the PALEOMAP Project (Scotese, 2011). To reconstruct paleo-ocean bathymetry, we use a plate model equation to model depth to basement (Turcotte and Schubert, 2002). We estimate plate model equation parameter values from measurements of modern oceans (Crosby et al., 2006). On top of the depth to basement, we isostatically add a multilayer sediment model derived from area-corrected sediment thickness data (Divins, 2003; Whittaker et al., 2013). Lastly, we parameterize the modern continental shelf, slope, and rise in a "sediment wedge model" to connect the coastline with the closest ocean crust as defined by Müller et al. (2008 a, b). These parameters are defined using empirical relationships obtained from study of modern ocean transects where a complete rifting history is preserved (Atlantic and Southern oceans), and the closest approach of the respective oceanic crust (Müller et al., 2008a,b) to the coastline. We use the modern ocean as a test, comparing maps and cross sections of modern ocean bathymetry modeled using our reconstruction method with that of ETOPO1 (Amante and Eakins, 2009). Adding sea plateaus and seamounts minimize the difference between our modeled bathymetry and ETOPO1. Finally, we also present a comparison of our reconstructed paleo-bathymetry to that of Müller et al. (2008 a,b) for the Cenomanian-Turonian (90 Ma). References: Amante, C., Eakins, B.W., 2009, NOAA Tech. Memo. NESDIS NGDC-24, 19 p. Crosby, A., McKenzie, D., Sclater, J.G., 2006, Geophysical Journal Int. 166.2, 553-573. Divins, D., 2003, NOAA NGDC, Boulder, CO. Müller, R., Sdrolias, M., Gaina, C., Roest, W., 2008b, Geochemistry, Geophysics, Geosystems, 9, Q04006, doi:10.1029/2007GC001743 Müller, R., Sdrolias, M., Gaina, C., Steinberger, B., Heine, C., 2008a, Science, 319, 1357-1362. Scotese, C., 2011, PALEOMAP Project, Arlington, Texas. Turcotte, D., Schubert, G., 2002, Cambridge University Press, Cambridge, 456 p. Whittaker, J., Goncharov, A., Williams, S., Müller, R., Leitchenkov, G., 2013, Geochemistry, Geophysics, Geosystems. DOI:10.1002/ggge.20181

Spectrum Analysis of Inertial and Subinertial Motions Based on Analyzed Winds and Wind-Driven Currents from a Primitive Equation General Ocean Circulation Model.

DTIC Science & Technology

1982-12-01

1Muter.Te Motions Based on Ana lyzed Winds and wind-driven December 1982 Currents from. a Primitive Squat ion General a.OW -love"*..* Oean Circulation...mew se"$ (comeS.... do oISN..u am ae~ 00do OWaor NUN Fourier and Rotary Spc , Analysis Modeled Inertial and Subinrtial Motion 4 Primitive Equation

Environmental Acoustics and Intensity Vector Acoustics with Emphasis on Shallow Water Effects and the Sea Surface

DTIC Science & Technology

2015-09-30

work of Ballard [2] in which horizontal (bathymetric) refraction is solved using a parabolic equation (PE) in Cartesian coordinates. In particular for...wedge-like ocean with penetrable ocean,” J. Acoust. Soc. Am., 82 198-210, 1987 [2] M. S. Ballard , “Modeling three-dimensional propagation in a

Swimming Behavior and Flow Geometry: A Fluid Mechanical Study of the Feeding Currents in Calanoid Copepods

NASA Astrophysics Data System (ADS)

Jiang, Houshuo; Meneveau, Charles; Osborn, Thomas R.

2003-11-01

Copepods are small crustaceans living in oceans and fresh waters and play an important role in the marine and freshwater food webs. As they are the biggest biomass in the oceans some call them "the insects of the sea". Previous laboratory observations have shown that the fluid mechanical phenomena occurring at copepod body scale are crucial for the survival of copepods. One of the interesting phenomena is that many calanoid copepods display various behaviors to create the feeding currents for the purpose of capturing food particles. We have developed a fluid mechanical model to study the feeding currents. The model is a self-propelled body model in that the Navier-Stokes equations are properly coupled with the dynamic equations for the copepod's body. The model has been solved both analytically using the Stokes approximation with a spherical body shape and numerically using CFD with a realistic body shape.

A numerical study of the string function using a primitive equation ocean model

NASA Astrophysics Data System (ADS)

Tyler, R. H.; Käse, R.

We use results from a primitive-equation ocean numerical model (SCRUM) to test a theoretical 'string function' formulation put forward by Tyler and Käse in another article in this issue. The string function acts as a stream function for the large-scale potential energy flow under the combined beta and topographic effects. The model results verify that large-scale anomalies propagate along the string function contours with a speed correctly given by the cross-string gradient. For anomalies having a scale similar to the Rossby radius, material rates of change in the layer mass following the string velocity are balanced by material rates of change in relative vorticity following the flow velocity. It is shown that large-amplitude anomalies can be generated when wind stress is resonant with the string function configuration.

Southwest Pacific Ocean Circulation and Climate Experiment (SPICE) scientific advances and future west pacific coordination

NASA Astrophysics Data System (ADS)

Ganachaud, A. S.; Sprintall, J.; Lin, X.; Ando, K.

2016-02-01

The Southwest Pacific Ocean Circulation and Climate Experiment (SPICE) is an international research program under the auspices of CLIVAR (Climate Variability and Predictability). The key objectives are to understand the Southwest Pacific Ocean circulation and Convergence Zone (SPCZ) dynamics, as well as their influence on regional and basin-scale climate patterns. It was designed to measure and monitor the ocean circulation, and to validate and improve numerical models. South Pacific oceanic waters are carried from the subtropical gyre centre in the westward flowing South Equatorial Current (SEC), towards the southwest Pacific-a major circulation pathway that redistributes water from the subtropics to the equator and Southern Ocean. Water transit through the Coral and Solomon Seas is potentially of great importance to tropical climate prediction because changes in either the temperature or the amount of water arriving at the equator have the capability to modulate ENSO and produce basin-scale climate feedbacks. On average, the oceanic circulation is driven by the Trade Winds, and subject to substantial variability, related with the SPCZ position and intensity. The circulation is complex, with the SEC splitting into zonal jets upon encountering island archipelagos, before joining either the East Australian Current or the New Guinea Costal UnderCurrent towards the equator. SPICE included large, coordinated in situ measurement programs and high resolution numerical simulations of the area. After 8 years of substantial in situ oceanic observational and modeling efforts, our understanding of the region has much improved. We have a refined description of the SPCZ behavior, boundary currents, pathways, and water mass transformation, including the previously undocumented Solomon Sea. The transports are large and vary substantially in a counter-intuitive way, with asymmetries and gating effects that depend on time scales. We will review the recent advancements and discuss our current knowledge gaps and important emerging research directions. In particular we will discuss how SPICE, along with the Northwestern Pacific Ocean Circulation and Climate Experiment (NPOCE) and Indonesian ThroughFlow (ITF) programs could evolve toward an integrative observing system under CLIVAR coordination.

A theoretical study of radar return and radiometric emission from the sea

NASA Technical Reports Server (NTRS)

Peake, W. H.

1972-01-01

The applicability of the various electromagnetic models of scattering from the ocean are reviewed. These models include the small perturbation method, the geometric optics solution, the composite model, and the exact integral equation solution. The restrictions on the electromagnetic models are discussed.

Optimisation of an idealised primitive equation ocean model using stochastic parameterization

NASA Astrophysics Data System (ADS)

Cooper, Fenwick C.

2017-05-01

Using a simple parameterization, an idealised low resolution (biharmonic viscosity coefficient of 5 × 1012 m4s-1 , 128 × 128 grid) primitive equation baroclinic ocean gyre model is optimised to have a much more accurate climatological mean, variance and response to forcing, in all model variables, with respect to a high resolution (biharmonic viscosity coefficient of 8 × 1010 m4s-1 , 512 × 512 grid) equivalent. For example, the change in the climatological mean due to a small change in the boundary conditions is more accurate in the model with parameterization. Both the low resolution and high resolution models are strongly chaotic. We also find that long timescales in the model temperature auto-correlation at depth are controlled by the vertical temperature diffusion parameter and time mean vertical advection and are caused by short timescale random forcing near the surface. This paper extends earlier work that considered a shallow water barotropic gyre. Here the analysis is extended to a more turbulent multi-layer primitive equation model that includes temperature as a prognostic variable. The parameterization consists of a constant forcing, applied to the velocity and temperature equations at each grid point, which is optimised to obtain a model with an accurate climatological mean, and a linear stochastic forcing, that is optimised to also obtain an accurate climatological variance and 5 day lag auto-covariance. A linear relaxation (nudging) is not used. Conservation of energy and momentum is discussed in an appendix.

Should tsunami simulations include a nonzero initial horizontal velocity?

NASA Astrophysics Data System (ADS)

Lotto, Gabriel C.; Nava, Gabriel; Dunham, Eric M.

2017-08-01

Tsunami propagation in the open ocean is most commonly modeled by solving the shallow water wave equations. These equations require initial conditions on sea surface height and depth-averaged horizontal particle velocity or, equivalently, horizontal momentum. While most modelers assume that initial velocity is zero, Y.T. Song and collaborators have argued for nonzero initial velocity, claiming that horizontal displacement of a sloping seafloor imparts significant horizontal momentum to the ocean. They show examples in which this effect increases the resulting tsunami height by a factor of two or more relative to models in which initial velocity is zero. We test this claim with a "full-physics" integrated dynamic rupture and tsunami model that couples the elastic response of the Earth to the linearized acoustic-gravitational response of a compressible ocean with gravity; the model self-consistently accounts for seismic waves in the solid Earth, acoustic waves in the ocean, and tsunamis (with dispersion at short wavelengths). Full-physics simulations of subduction zone megathrust ruptures and tsunamis in geometries with a sloping seafloor confirm that substantial horizontal momentum is imparted to the ocean. However, almost all of that initial momentum is carried away by ocean acoustic waves, with negligible momentum imparted to the tsunami. We also compare tsunami propagation in each simulation to that predicted by an equivalent shallow water wave simulation with varying assumptions regarding initial velocity. We find that the initial horizontal velocity conditions proposed by Song and collaborators consistently overestimate the tsunami amplitude and predict an inconsistent wave profile. Finally, we determine tsunami initial conditions that are rigorously consistent with our full-physics simulations by isolating the tsunami waves from ocean acoustic and seismic waves at some final time, and backpropagating the tsunami waves to their initial state by solving the adjoint problem. The resulting initial conditions have negligible horizontal velocity.[Figure not available: see fulltext.

Parameterization of typhoon-induced ocean cooling using temperature equation and machine learning algorithms: an example of typhoon Soulik (2013)

NASA Astrophysics Data System (ADS)

Wei, Jun; Jiang, Guo-Qing; Liu, Xin

2017-09-01

This study proposed three algorithms that can potentially be used to provide sea surface temperature (SST) conditions for typhoon prediction models. Different from traditional data assimilation approaches, which provide prescribed initial/boundary conditions, our proposed algorithms aim to resolve a flow-dependent SST feedback between growing typhoons and oceans in the future time. Two of these algorithms are based on linear temperature equations (TE-based), and the other is based on an innovative technique involving machine learning (ML-based). The algorithms are then implemented into a Weather Research and Forecasting model for the simulation of typhoon to assess their effectiveness, and the results show significant improvement in simulated storm intensities by including ocean cooling feedback. The TE-based algorithm I considers wind-induced ocean vertical mixing and upwelling processes only, and thus obtained a synoptic and relatively smooth sea surface temperature cooling. The TE-based algorithm II incorporates not only typhoon winds but also ocean information, and thus resolves more cooling features. The ML-based algorithm is based on a neural network, consisting of multiple layers of input variables and neurons, and produces the best estimate of the cooling structure, in terms of its amplitude and position. Sensitivity analysis indicated that the typhoon-induced ocean cooling is a nonlinear process involving interactions of multiple atmospheric and oceanic variables. Therefore, with an appropriate selection of input variables and neuron sizes, the ML-based algorithm appears to be more efficient in prognosing the typhoon-induced ocean cooling and in predicting typhoon intensity than those algorithms based on linear regression methods.

On the time-splitting scheme used in the Princeton Ocean Model

NASA Astrophysics Data System (ADS)

Kamenkovich, V. M.; Nechaev, D. A.

2009-05-01

The analysis of the time-splitting procedure implemented in the Princeton Ocean Model (POM) is presented. The time-splitting procedure uses different time steps to describe the evolution of interacting fast and slow propagating modes. In the general case the exact separation of the fast and slow modes is not possible. The main idea of the analyzed procedure is to split the system of primitive equations into two systems of equations for interacting external and internal modes. By definition, the internal mode varies slowly and the crux of the problem is to determine the proper filter, which excludes the fast component of the external mode variables in the relevant equations. The objective of this paper is to examine properties of the POM time-splitting procedure applied to equations governing the simplest linear non-rotating two-layer model of constant depth. The simplicity of the model makes it possible to study these properties analytically. First, the time-split system of differential equations is examined for two types of the determination of the slow component based on an asymptotic approach or time-averaging. Second, the differential-difference scheme is developed and some criteria of its stability are discussed for centered, forward, or backward time-averaging of the external mode variables. Finally, the stability of the POM time-splitting schemes with centered and forward time-averaging is analyzed. The effect of the Asselin filter on solutions of the considered schemes is studied. It is assumed that questions arising in the analysis of the simplest model are inherent in the general model as well.

Energy-optimal path planning in the coastal ocean

NASA Astrophysics Data System (ADS)

Subramani, Deepak N.; Haley, Patrick J.; Lermusiaux, Pierre F. J.

2017-05-01

We integrate data-driven ocean modeling with the stochastic Dynamically Orthogonal (DO) level-set optimization methodology to compute and study energy-optimal paths, speeds, and headings for ocean vehicles in the Middle-Atlantic Bight (MAB) region. We hindcast the energy-optimal paths from among exact time-optimal paths for the period 28 August 2006 to 9 September 2006. To do so, we first obtain a data-assimilative multiscale reanalysis, combining ocean observations with implicit two-way nested multiresolution primitive-equation simulations of the tidal-to-mesoscale dynamics in the region. Second, we solve the reduced-order stochastic DO level-set partial differential equations (PDEs) to compute the joint probability of minimum arrival time, vehicle-speed time series, and total energy utilized. Third, for each arrival time, we select the vehicle-speed time series that minimize the total energy utilization from the marginal probability of vehicle-speed and total energy. The corresponding energy-optimal path and headings are obtained through the exact particle-backtracking equation. Theoretically, the present methodology is PDE-based and provides fundamental energy-optimal predictions without heuristics. Computationally, it is 3-4 orders of magnitude faster than direct Monte Carlo methods. For the missions considered, we analyze the effects of the regional tidal currents, strong wind events, coastal jets, shelfbreak front, and other local circulations on the energy-optimal paths. Results showcase the opportunities for vehicles that intelligently utilize the ocean environment to minimize energy usage, rigorously integrating ocean forecasting with optimal control of autonomous vehicles.

How does ice sheet loading affect ocean flow around Antarctica?

NASA Astrophysics Data System (ADS)

Dijkstra, H. A.; Rugenstein, M. A.; Stocchi, P.; von der Heydt, A. S.

2012-12-01

Interactions and dynamical feedbacks between ocean circulation, heat and atmospheric moisture transport, ice sheet evolution, and Glacial Isostatic Adjustment (GIA) are overlooked issues in paleoclimatology. Here we will present first results on how ocean flows were possibly affected by the glaciation of Antarctica across the Eocene-Oligocene Transition (~ 34 Ma) through GIA and bathymetry variations. GIA-induced gravitationally self-consistent bathymetry variations are determined by solving the Sea Level Equation (SLE), which describes the time dependent shape of (i) the solid Earth and (ii) the equipotential surface of gravity. Since the ocean circulation equations are defined relative to the equipotential surface of gravity, only bathymetry variations can influence ocean flows, although the sea surface slope will also change through time due to gravitational attraction. We use the Hallberg Isopycnal Model under late Eocene conditions to calculate equilibrium ocean flows in a domain in which the bathymetry evolves under ice loading according to the SLE. The bathymetric effects of the glaciation of Antarctica lead to substantial spatial changes in ocean flows, and close to the coast, the flow even reverses direction. Volume transports through the Drake Passage and Tasman Seaway adjust to the new bathymetry. The results indicate that GIA-induced ocean flow variations alone may have had an impact on sedimentation and erosion patterns, the repositioning of fronts, ocean heat transport and grounding line and ice sheet stability.

Dynamics, Heat Transport, Spectral Composition and Acoustic Signatures of Mesoscale Variability in the Ocean

DTIC Science & Technology

2013-12-01

Eastward background flow EOS Equation of state GDEM Generalized Digital Environmental Model GRB Growth Rate Balance model HPCMP High Performance...the Naval Research Lab (NRL) Generalized Digital Environmental Model ( GDEM ). This provides a realistic and detailed profile for a known turbulent

Impact of errors in short wave radiation and its attenuation on modeled upper ocean heat content

DTIC Science & Technology

Photosynthetically available radiation (PAR) and its attenuation with the depth represent a forcing (source) term in the governing equation for the...and vertical attenuation of PAR have on the upper ocean model heat content. In the Monterey Bay area, we show that with a decrease in water clarity...attenuation coefficient. For Jerlov’s type IA water (attenuation coefficient is 0.049 m1), the relative error in surface PAR introduces an error

A delay differential model of ENSO variability: parametric instability and the distribution of extremes

NASA Astrophysics Data System (ADS)

Ghil, M.; Zaliapin, I.; Thompson, S.

2008-05-01

We consider a delay differential equation (DDE) model for El-Niño Southern Oscillation (ENSO) variability. The model combines two key mechanisms that participate in ENSO dynamics: delayed negative feedback and seasonal forcing. We perform stability analyses of the model in the three-dimensional space of its physically relevant parameters. Our results illustrate the role of these three parameters: strength of seasonal forcing b, atmosphere-ocean coupling κ, and propagation period τ of oceanic waves across the Tropical Pacific. Two regimes of variability, stable and unstable, are separated by a sharp neutral curve in the (b, τ) plane at constant κ. The detailed structure of the neutral curve becomes very irregular and possibly fractal, while individual trajectories within the unstable region become highly complex and possibly chaotic, as the atmosphere-ocean coupling κ increases. In the unstable regime, spontaneous transitions occur in the mean "temperature" (i.e., thermocline depth), period, and extreme annual values, for purely periodic, seasonal forcing. The model reproduces the Devil's bleachers characterizing other ENSO models, such as nonlinear, coupled systems of partial differential equations; some of the features of this behavior have been documented in general circulation models, as well as in observations. We expect, therefore, similar behavior in much more detailed and realistic models, where it is harder to describe its causes as completely.

Modeling Vertical Structure and Heat Transport within the Oceans of Ice-covered Worlds (Invited)

NASA Astrophysics Data System (ADS)

Goodman, J. C.

2010-12-01

Indirect observational evidence provides a strong case for liquid oceans beneath the icy crust of Europa and several other frozen moons in the outer solar system. However, little is known about the fluid circulation within these exotic oceans. As a first step toward understanding circulations driven by buoyancy (rather than mechanical forcing from tides), one must understand the typical vertical structure of temperature, salinity, and thus density within the ocean. Following a common approach from terrestrial oceanography, I have built a "single column convection model" for icy world oceans, which describes the density structure of the ocean as a function of depth only: horizontal variations are ignored. On Earth, this approach is of limited utility, because of the strong influence of horizontal wind-driven currents and sea-surface temperature gradients set in concert with the overlying atmosphere. Neither of these confounding issues is present in an icy world's ocean. In the model, mixing of fluid properties via overturning convection is modeled as a strong diffusive process which only acts when the ocean is vertically unstable. "Double diffusive" processes (salt fingering and diffusive layering) are included: these are mixing processes resulting from the unequal molecular diffusivities of heat and salt. Other important processes, such as heating on adiabatic compression, and freshwater fluxes from melting overlying ice, are also included. As a simple test case, I considered an ocean of Europa-like depth (~100 km) and gravity, heated from the seafloor. To simplify matters, I specified an equation of state appropriate to terrestrial seawater, and a simple isothermal ocean as an initial condition. As expected, convection gradually penetrates upward, warming the ocean to an adiabatic, unstratified equilibrium density profile on a timescale of 50 kyr if 4.5 TW of heat are emitted by the silicate interior; the same result is achieved in proportionally more/less time for weaker/stronger internal heating. Unlike Earth's oceans, I predict that since icy worlds' oceans are heated from below, they will generally be unstratified, with constant potential density from top to bottom. There will be no pycnocline as on Earth, so global ocean currents supported by large-scale density gradients seem unlikely. However, icy world oceans may be "weird" in ways which are unheard-of in terrestrial oceanography The density of sulfate brine has a very different equation of state than chloride brines: does this affect the vertical structure? If the ocean water is very pure, cold water can be less dense than warm. Can this lead to periodic catastrophic overturning, as proposed by other authors? These and other questions are currently being investigated using the single-column convection model as a primary tool.

Twin Convergence Zones

NASA Technical Reports Server (NTRS)

2002-01-01

NASA's QuikSCAT satellite has confirmed a 30-year old largely unproven theory that there are two areas near the equator where the winds converge year after year and drive ocean circulation south of the equator. By analyzing winds, QuikSCAT has found a year-round southern and northern Intertropical Convergence Zone. This find is important to climate modelers and weather forecasters because it provides more detail on how the oceans and atmosphere interact near the equator. The Intertropical Convergence Zone (ITCZ) is the region that circles the Earth near the equator, where the trade winds of both the Northern and Southern Hemispheres come together. North of the equator, strong sun and warm water of the equator heats the air in the ITCZ, drawing air in from north and south and causing the air to rise. As the air rises it cools, releasing the accumulated moisture in an almost perpetual series of thunderstorms. Satellite data, however, has confirmed that there is an ITCZ north of the equator and a parallel ITCZ south of the equator. Variation in the location of the ITCZ is important to people around the world because it affects the north-south atmospheric circulation, which redistributes energy. It drastically affects rainfall in many equatorial nations, resulting in the wet and dry seasons of the tropics rather than the cold and warm seasons of higher latitudes. Longer term changes in the ITCZ can result in severe droughts or flooding in nearby areas. 'The double ITCZ is usually only identified in the Pacific and Atlantic Oceans on a limited and seasonal basis,' said Timothy Liu, of NASA's Jet Propulsion Laboratory and California Institute of Technology, Pasadena, Calif., and lead researcher on the project. In the eastern Pacific Ocean, the southern ITCZ is usually seen springtime. In the western Atlantic Ocean, the southern ITCZ was recently clearly identified only in the summertime. However, QuikSCAT's wind data has seen the southern ITCZ in all seasons across the entire Atlantic Ocean and the eastern Pacific. 'QuikSCAT's wind data confirms there is a double ITCZ, and that they exist all year long,' Liu said. This is a major find for the science community, as the existence, location, and seasonality of the double ITCZ had remained controversial since 1969. full text: Satellite Sees Double Zones of Converging Tropical Winds around The World For more about convergence zones, read: The Intertropical Convergence Zone and Convergence Zones: Where the Action Is Image courtesy Liu and Xie, NASA JPL

Response of the tropical Pacific Ocean to El Niño versus global warming

DOE Office of Scientific and Technical Information (OSTI.GOV)

Liu, Fukai; Luo, Yiyong; Lu, Jian

Climate models project an El Niño-like SST response in the tropical Pacific Ocean to global warming (GW). By employing the Community Earth System Model (CESM) and applying an overriding technique to its ocean component, Parallel Ocean Program version 2 (POP2), this study investigates the similarity and difference of formation mechanism for the changes in the tropical Pacific Ocean under El Niño and GW. Results show that, despite sharing some similarities between the two scenarios, there are many significant distinctions between GW and El Niño: 1) the phase locking of the seasonal cycle reduction is more notable under GW compared withmore » El Niño, implying more extreme El Niño events in the future; 2) in contrast to the penetration of the equatorial subsurface temperature anomaly that appears to propagate in the form of an oceanic equatorial upwelling Kelvin wave during El Niño, the GW-induced subsurface temperature anomaly manifest in the form of off-equatorial upwelling Rossby waves; 3) while significant across-equator northward heat transport (NHT) is induced by the wind stress anomalies associated with El Niño, little NHT is found at the equator due to a symmetric change in the shallow meridional overturning circulation that appears to be weakened in both North and South Pacific under GW; and 4) the maintaining mechanisms for the eastern equatorial Pacific warming are also substantially different.« less

Stochastic Forcing for High-Resolution Regional and Global Ocean and Atmosphere-Ocean Coupled Ensemble Forecast System

NASA Astrophysics Data System (ADS)

Rowley, C. D.; Hogan, P. J.; Martin, P.; Thoppil, P.; Wei, M.

2017-12-01

An extended range ensemble forecast system is being developed in the US Navy Earth System Prediction Capability (ESPC), and a global ocean ensemble generation capability to represent uncertainty in the ocean initial conditions has been developed. At extended forecast times, the uncertainty due to the model error overtakes the initial condition as the primary source of forecast uncertainty. Recently, stochastic parameterization or stochastic forcing techniques have been applied to represent the model error in research and operational atmospheric, ocean, and coupled ensemble forecasts. A simple stochastic forcing technique has been developed for application to US Navy high resolution regional and global ocean models, for use in ocean-only and coupled atmosphere-ocean-ice-wave ensemble forecast systems. Perturbation forcing is added to the tendency equations for state variables, with the forcing defined by random 3- or 4-dimensional fields with horizontal, vertical, and temporal correlations specified to characterize different possible kinds of error. Here, we demonstrate the stochastic forcing in regional and global ensemble forecasts with varying perturbation amplitudes and length and time scales, and assess the change in ensemble skill measured by a range of deterministic and probabilistic metrics.

Impact of Seawater Nonlinearities on Nordic Seas Circulation

NASA Astrophysics Data System (ADS)

Helber, R. W.; Wallcraft, A. J.; Shriver, J. F.

2017-12-01

The Nordic Seas (Greenland, Iceland, and Norwegian Seas) form an ocean basin important for Arctic-mid-latitude climate linkages. Cold fresh water from the Arctic Ocean and warm salty water from the North Atlantic Ocean meet in the Nordic Seas, where a delicate balance between temperature and salinity variability results in deep water formation. Seawater non-linearities are stronger at low temperatures and salinities making high-latitude oceans highly subject to thermbaricity and cabbeling. This presentation highlights and quantifies the impact of seawater non-linearities on the Nordic Seas circulation. We use two layered ocean circulation models, the Hybrid Coordinate Ocean Model (HYOCM) and the Modular Ocean Model version 6 (MOM6), that enable accurate representation of processes along and across density or neutral density surfaces. Different equations-of-state and vertical coordinates are evaluated to clarify the impact of seawater non-linearities. Present Navy systems, however, do not capture some features in the Nrodic Seas vertical structure. For example, observations from the Greenland Sea reveal a subsurface temperature maximum that deepens from approximately 1500 m during 1998 to 1800 m during 2005. We demonstrate that in terms of density, salinity is the largest source of error in Nordic Seas Navy forecasts, regional scale models can represent mesoscale features driven by thermobaricity, vertical coordinates are a critical issue in Nordic Sea circulation modeling.

Preliminary testing of turbulence and radionuclide transport modeling in deep ocean environment

DOE Office of Scientific and Technical Information (OSTI.GOV)

Onishi, Y.; Dummuller, D.C.; Trent, D.S.

Pacific Northwest Laboratory (PNL) performed a study for the US Environmental Protection Agency's Office of Radiation Programs to (1) identify candidate models for regional modeling of low-level waste ocean disposal sites in the mid-Atlantic ocean; (2) evaluate mathematical representation of the model's eddy viscosity/dispersion coefficients; and (3) evaluate the adequacy of the k-{epsilon} turbulence model and the feasibility of one of the candidate models, TEMPEST{copyright}/FLESCOT{copyright}, to deep-ocean applications on a preliminary basis. PNL identified the TEMPEST{copyright}/FLESCOT{copyright}, FLOWER, Blumberg's, and RMA 10 models as appropriate candidates for the regional radionuclide modeling. Among these models, TEMPEST/FLESCOT is currently the only model thatmore » solves distributions of flow, turbulence (with the k-{epsilon} model), salinity, water temperature, sediment, dissolved contaminants, and sediment-sorbed contaminants. Solving the Navier-Stokes equations using higher order correlations is not practical for regional modeling because of the prohibitive computational requirements; therefore, the turbulence modeling is a more practical approach. PNL applied the three-dimensional code, TEMPEST{copyright}/FLESCOT{copyright} with the k-{epsilon} model, to a very simple, hypothetical, two-dimensional, deep-ocean case, producing at least qualitatively appropriate results. However, more detailed testing should be performed for the further testing of the code. 46 refs., 39 figs., 6 tabs.« less

Using Vertically Integrated Ocean Fields to Characterize Greenland Icebergs' Distribution and Lifetime

NASA Astrophysics Data System (ADS)

Marson, Juliana M.; Myers, Paul G.; Hu, Xianmin; Le Sommer, Julien

2018-05-01

Icebergs represent approximately half of Greenland's yearly mass loss, having important implications for biological productivity, freshwater fluxes in the ocean, and navigation. This study applies an iceberg model that uses integrated ocean fields (from surface to iceberg keel) to simulate the drift and decay of Greenland icebergs. This version of iceberg model (VERT) is compared with a more widely adopted version (SURF) which only uses surface ocean fields in its equations. We show that icebergs in VERT tend to drift along the shelf break, while in SURF they concentrate along the coastline. Additionally, we show that Greenland's southeast coast is the source of ˜60% of the icebergs that cross the interior of the Labrador Sea—a region that stages buoyancy-driven convection and is, therefore, sensitive to freshwater input.

A Preliminary Model Study of the Large-Scale Seasonal Cycle in Bottom Pressure Over the Global Ocean

NASA Technical Reports Server (NTRS)

Ponte, Rui M.

1998-01-01

Output from the primitive equation model of Semtner and Chervin is used to examine the seasonal cycle in bottom pressure (Pb) over the global ocean. Effects of the volume-conserving formulation of the model on the calculation Of Pb are considered. The estimated seasonal, large-scale Pb signals have amplitudes ranging from less than 1 cm over most of the deep ocean to several centimeters over shallow, boundary regions. Variability generally increases toward the western sides of the basins, and is also larger in some Southern Ocean regions. An oscillation between subtropical and higher latitudes in the North Pacific is clear. Comparison with barotropic simulations indicates that, on basin scales, seasonal Pb variability is related to barotropic dynamics and the seasonal cycle in Ekman pumping, and results from a small, net residual in mass divergence from the balance between Ekman and Sverdrup flows.

A two-dimensional ocean model for long-term climatic simulations: Stability and coupling to atmospheric and sea ice models

NASA Astrophysics Data System (ADS)

Harvey, L. D. Danny

1992-06-01

A two-dimensional (latitude-depth) deep ocean model is presented which is coupled to a sea ice model and an Energy Balance Climate Model (EBCM), the latter having land-sea and surface-air resolution. The processes which occur in the ocean model are thermohaline overturning driven by the horizontal density gradient, shallow wind-driven overturning cells, convective overturning, and vertical and horizontal diffusion of heat and salt. The density field is determined from the temperature and salinity fields using a nonlinear equation of state. Mixed layer salinity is affected by evaporation, precipitation, runoff from continents, and sea ice freezing and melting, as well as by advective, convective, and diffusive exchanges with the deep ocean. The ocean model is first tested in an uncoupled mode, in which hemispherically symmetric mixed layer temperature and salinity, or salinity flux, are specified as upper boundary conditions. An experiment performed with previous models is repeated in which a mixed layer salinity perturbation is introduced in the polar half of one hemisphere after switching from a fixed salinity to a fixed salinity flux boundary condition. For small values of the vertical diffusion coefficient KV, the model undergoes self-sustained oscillations with a period of about 1500 years. With larger values of KV, the model locks into either an asymmetric mode with a single overturning cell spanning both hemispheres, or a symmetric quiescent state with downwelling near the equator, upwelling at high latitudes, and a warm deep ocean (depending on the value of KV). When the ocean model is forced with observed mixed layer temperature and salinity, no oscillations occur. The model successfully simulates the very weak meridional overturning and strong Antarctic Circumpolar Current at the latitudes of the Drake Passage. The coupled EBCM-deep ocean model displays internal oscillations with a period of 3000 years if the ocean fraction is uniform with latitude and KV and the horizontal diffusion coefficient in the mixed layer are not too large. Globally averaged atmospheric temperature changes of 2 K are driven by oscillations in the heat flux into or out of the deep ocean, with the sudden onset of a heat flux out of the deep ocean associated with the rapid onset of thermohaline overturning after a quiescent period, and the sudden onset of a heat flux into the deep ocean associated with the collapse of thermohaline overturning. When the coupled model is run with prescribed parameters (such as land-sea fraction and precipitation) varying with latitude based on observations, the model does not oscillate and produces a reasonable deep ocean temperature field but a completely unrealistic salinity field. Resetting the mixed layer salinity to observations on each time step (equivalent to the "flux correction" method used in atmosphere-ocean general circulation models) is sufficient to give a realistic salinity field throughout the ocean depth, but dramatically alters the flow field and associated heat transport. Although the model is highly idealized, the finding that the maximum perturbation in globally averaged heat flux from the deep ocean to the surface over a 100-year period is 1.4 W m-2 suggests that effect of continuing greenhouse gas increases, which could result in a heating perturbation of 10 W m-2 by the end of the next century, will swamp possible surface heating perturbations due to changes in oceanic circulation. On the other hand, the extreme sensitivity of the oceanic flow field to variations in precipitation and evaporation suggests that it will not be possible to produce accurate projections of regional climatic change in the near term, if at all.

Application of a Reduced Order Kalman Filter to Initialize a Coupled Atmosphere-Ocean Model: Impact on the Prediction of El Nino

NASA Technical Reports Server (NTRS)

Ballabrera-Poy, J.; Busalacchi, A.; Murtugudde, R.

2000-01-01

A reduced order Kalman Filter, based on a simplification of the Singular Evolutive Extended Kalman (SEEK) filter equations, is used to assimilate observed fields of the surface wind stress, sea surface temperature and sea level into the nonlinear coupled ocean-atmosphere model of Zebiak and Cane. The SEEK filter projects the Kalman Filter equations onto a subspace defined by the eigenvalue decomposition of the error forecast matrix, allowing its application to high dimensional systems. The Zebiak and Cane model couples a linear reduced gravity ocean model with a single vertical mode atmospheric model of Zebiak. The compatibility between the simplified physics of the model and each observed variable is studied separately and together. The results show the ability of the model to represent the simultaneous value of the wind stress, SST and sea level, when the fields are limited to the latitude band 10 deg S - 10 deg N In this first application of the Kalman Filter to a coupled ocean-atmosphere prediction model, the sea level fields are assimilated in terms of the Kelvin and Rossby modes of the thermocline depth anomaly. An estimation of the error of these modes is derived from the projection of an estimation of the sea level error over such modes. This method gives a value of 12 for the error of the Kelvin amplitude, and 6 m of error for the Rossby component of the thermocline depth. The ability of the method to reconstruct the state of the equatorial Pacific and predict its time evolution is demonstrated. The method is shown to be quite robust for predictions up to six months, and able to predict the onset of the 1997 warm event fifteen months before its occurrence.

Application of a Reduced Order Kalman Filter to Initialize a Coupled Atmosphere-Ocean Model: Impact on the Prediction of El Nino

NASA Technical Reports Server (NTRS)

Ballabrera-Poy, Joaquim; Busalacchi, Antonio J.; Murtugudde, Ragu

2000-01-01

A reduced order Kalman Filter, based on a simplification of the Singular Evolutive Extended Kalman (SEEK) filter equations, is used to assimilate observed fields of the surface wind stress, sea surface temperature and sea level into the nonlinear coupled ocean-atmosphere model. The SEEK filter projects the Kalman Filter equations onto a subspace defined by the eigenvalue decomposition of the error forecast matrix, allowing its application to high dimensional systems. The Zebiak and Cane model couples a linear reduced gravity ocean model with a single vertical mode atmospheric model of Zebiak. The compatibility between the simplified physics of the model and each observed variable is studied separately and together. The results show the ability of the model to represent the simultaneous value of the wind stress, SST and sea level, when the fields are limited to the latitude band 10 deg S - 10 deg N. In this first application of the Kalman Filter to a coupled ocean-atmosphere prediction model, the sea level fields are assimilated in terms of the Kelvin and Rossby modes of the thermocline depth anomaly. An estimation of the error of these modes is derived from the projection of an estimation of the sea level error over such modes. This method gives a value of 12 for the error of the Kelvin amplitude, and 6 m of error for the Rossby component of the thermocline depth. The ability of the method to reconstruct the state of the equatorial Pacific and predict its time evolution is demonstrated. The method is shown to be quite robust for predictions I up to six months, and able to predict the onset of the 1997 warm event fifteen months before its occurrence.

Parallel Computation of Ocean-Atmosphere-Wave Coupled Storm Surge Model

NASA Astrophysics Data System (ADS)

Kim, K.; Yamashita, T.

2003-12-01

Ocean-atmosphere interactions are very important in the formation and development of tropical storms. These interactions are dominant in exchanging heat, momentum, and moisture fluxes. Heat flux is usually computed using a bulk equation. In this equation air-sea interface supplies heat energy to the atmosphere and to the storm. Dynamical interaction is most often one way in which it is the atmosphere that drives the ocean. The winds transfer momentum to both ocean surface waves and ocean current. The wind wave makes an important role in the exchange of the quantities of motion, heat and a substance between the atmosphere and the ocean. Storm surges can be considered as the phenomena of mean sea-level changes, which are the result of the frictional stresses of strong winds blowing toward the land and causing the set level and the low atmospheric pressure at the centre of the cyclone can additionally raise the sea level. In addition to the rise in water level itself, another wave factor must be considered. A rise of mean sea level due to white-cap wave dissipation should be considered. In bounded bodies of water, such as small seas, wind driven sea level set up is much serious than inverted barometer effects, in which the effects of wind waves on wind-driven current play an important role. It is necessary to develop the coupled system of the full spectral third-generation wind-wave model (WAM or WAVEWATCH III), the meso-scale atmosphere model (MM5) and the coastal ocean model (POM) for simulating these physical interactions. As the component of coupled system is so heavy for personal usage, the parallel computing system should be developed. In this study, first, we developed the coupling system of the atmosphere model, ocean wave model and the coastal ocean model, in the Beowulf System, for the simulation of the storm surge. It was applied to the storm surge simulation caused by Typhoon Bart (T9918) in the Yatsushiro Sea. The atmosphere model and the ocean model have been made the parallel codes by SPMD methods. The wave-current interface model was developed by defining the wave breaking stresses. And we developed the coupling program to collect and distribute the exchanging data with the parallel system. Every models and coupler are executed at same time, and they calculate own jobs and pass data with organic system. MPMD method programming was performed to couple the models. The coupler and each models united by the separated group, and they calculated by the group unit. Also they passed message when exchanging data by global unit. The data are exchanged every 60-second model time that is the least common multiple time of the atmosphere model, the wave model and the ocean model. The model was applied to the storm surge simulation in the Yatsushiro Sea, in which we could not simulated the observed maximum surge height with the numerical model that did not include the wave breaking stress. It is confirmed that the simulation which includes the wave breaking stress effects can produce the observed maximum height, 450 cm, at Matsuai.

Calculation of wind-driven surface currents in the North Atlantic Ocean

NASA Technical Reports Server (NTRS)

Rees, T. H.; Turner, R. E.

1976-01-01

Calculations to simulate the wind driven near surface currents of the North Atlantic Ocean are described. The primitive equations were integrated on a finite difference grid with a horizontal resolution of 2.5 deg in longitude and latitude. The model ocean was homogeneous with a uniform depth of 100 m and with five levels in the vertical direction. A form of the rigid-lid approximation was applied. Generally, the computed surface current patterns agreed with observed currents. The development of a subsurface equatorial countercurrent was observed.

The importance of planetary rotation period for ocean heat transport.

PubMed

Cullum, J; Stevens, D; Joshi, M

2014-08-01

The climate and, hence, potential habitability of a planet crucially depends on how its atmospheric and ocean circulation transports heat from warmer to cooler regions. However, previous studies of planetary climate have concentrated on modeling the dynamics of atmospheres, while dramatically simplifying the treatment of oceans, which neglects or misrepresents the effect of the ocean in the total heat transport. Even the majority of studies with a dynamic ocean have used a simple so-called aquaplanet that has no continental barriers, which is a configuration that dramatically changes the ocean dynamics. Here, the significance of the response of poleward ocean heat transport to planetary rotation period is shown with a simple meridional barrier--the simplest representation of any continental configuration. The poleward ocean heat transport increases significantly as the planetary rotation period is increased. The peak heat transport more than doubles when the rotation period is increased by a factor of ten. There are also significant changes to ocean temperature at depth, with implications for the carbon cycle. There is strong agreement between the model results and a scale analysis of the governing equations. This result highlights the importance of both planetary rotation period and the ocean circulation when considering planetary habitability.

Adaptive scaling model of the main pycnocline and the associated overturning circulation

NASA Astrophysics Data System (ADS)

Fuckar, Neven-Stjepan

This thesis examines a number of crucial factors and processes that control the structure of the main pycnocline and the associated overturning circulation that maintains the ocean stratification. We construct an adaptive scaling model: a semi-empirical low-order theory based on the total transformation balance that linearly superimposes parameterized transformation rate terms of various mechanisms that participate in the water-mass conversion between the warm water sphere and the cold water sphere. The depth of the main pycnocline separates the light-water domain from the dense-water domain beneath the surface, hence we introduce a new definition in an integral form that is dynamically based on the large-scale potential vorticity (i.e., vertical density gradient is selected for the kernel function of the normalized vertical integral). We exclude the abyssal pycnocline from our consideration and limit our domain of interest to the top 2 km of water column. The goal is to understand the controlling mechanisms, and analytically predict and describe a wide spectrum of ocean steady states in terms of key large-scale indices relevant for understanding the ocean's role in climate. A devised polynomial equation uses the average depth of the main pycnocline as a single unknown (the key vertical scale of the upper ocean stratification) and gives us an estimate for the northern hemisphere deep water production and export across the equator from the parts of this equation. The adaptive scaling model aims to elucidate the roles of a limited number of dominant processes that determine some key upper ocean circulation and stratification properties. Additionally, we use a general circulation model in a series of simplified single-basin ocean configurations and surface forcing fields to confirm the usefulness of our analytical model and further clarify several aspects of the upper ocean structure. An idealized numerical setup, containing all the relevant physical and dynamical properties, is key to obtaining a clear understanding, uncomplicated by the effect of the real world geometry or intricacy of realistic surface radiative and turbulent fluxes. We show that wind-driven transformation processes can be decomposed into two terms separately driven by the mid-latitude westerlies and the low-latitude easterlies. Our analytical model smoothly connects all the classical limits describing different ocean regimes in a single-basin single-hemisphere geometry. The adjective "adaptive" refers to a simple and quantitatively successful adjustment to the description of a single-basin two-hemisphere ocean, with and without a circumpolar channel under the hemispherically symmetric surface buoyancy. For example, our water-mass conversion framework, unifying wind-driven and thermohaline processes, provides us with further insight into the "Drake Passage effect without Drake Passage". The modification of different transformation pathways in the Southern Hemisphere results in the equivalent net conversion changes. The introduction of hemispheric asymmetry in the surface density can lead to significant hemispheric differences in the main pycnocline structure. This demonstrates the limitations of our analytical model based on only one key vertical scale. Also, we show a strong influence of the northern hemisphere surface density change in high latitudes on the southern hemisphere stratification and circumpolar transport.

Longitudinal biases in the Seychelles Dome simulated by 35 ocean-atmosphere coupled general circulation models

NASA Astrophysics Data System (ADS)

Nagura, Motoki; Sasaki, Wataru; Tozuka, Tomoki; Luo, Jing-Jia; Behera, Swadhin K.; Yamagata, Toshio

2013-02-01

Seychelles Dome refers to the shallow climatological thermocline in the southwestern Indian Ocean, where ocean wave dynamics efficiently affect sea surface temperature, allowing sea surface temperature anomalies to be predicted up to 1-2 years in advance. Accurate reproduction of the dome by ocean-atmosphere coupled general circulation models (CGCMs) is essential for successful seasonal predictions in the Indian Ocean. This study examines the Seychelles Dome as simulated by 35 CGCMs, including models used in phase five of the Coupled Model Intercomparison Project (CMIP5). Among the 35 CGCMs, 14 models erroneously produce an upwelling dome in the eastern half of the basin whereas the observed Seychelles Dome is located in the southwestern tropical Indian Ocean. The annual mean Ekman pumping velocity in these models is found to be almost zero in the southern off-equatorial region. This result is inconsistent with observations, in which Ekman upwelling acts as the main cause of the Seychelles Dome. In the models reproducing an eastward-displaced dome, easterly biases are prominent along the equator in boreal summer and fall, which result in shallow thermocline biases along the Java and Sumatra coasts via Kelvin wave dynamics and a spurious upwelling dome in the region. Compared to the CMIP3 models, the CMIP5 models are even worse in simulating the dome longitudes.

Global ocean tide mapping using TOPEX/Poseidon altimetry

NASA Technical Reports Server (NTRS)

Sanchez, Braulio V.; Cartwright, D. E.; Estes, R. H.; Williamson, R. G.; Colombo, O. L.

1991-01-01

The investigation's main goals are to produce accurate tidal maps of the main diurnal, semidiurnal, and long-period tidal components in the world's deep oceans. This will be done by the application of statistical estimation techniques to long time series of altimeter data provided by the TOPEX/POSEIDON mission, with additional information provided by satellite tracking data. In the prelaunch phase, we will use in our simulations and preliminary work data supplied by previous oceanographic missions, such as Seasat and Geosat. These results will be of scientific interest in themselves. The investigation will also be concerned with the estimation of new values, and their uncertainties, for tidal currents and for the physical parameters appearing in the Laplace tidal equations, such as bottom friction coefficients and eddy viscosity coefficients. This will be done by incorporating the altimetry-derived charts of vertical tides as boundary conditions in the integration of those equations. The methodology of the tidal representation will include the use of appropriate series expansions such as ocean-basin normal modes and spherical harmonics. The results of the investigation will be space-determined tidal models of coverage and accuracy superior to that of the present numerical models of the ocean tides, with the concomitant benefits to oceanography and associated disciplinary fields.

Constraints on Oceanic Meridional Transport of Heat and Carbon from Combined Oceanic and Atmospheric Measurements.

NASA Astrophysics Data System (ADS)

Resplandy, L.; Keeling, R. F.; Stephens, B. B.; Bent, J. D.; Jacobson, A. R.; Rödenbeck, C.; Khatiwala, S.

2016-02-01

The global ocean transports heat northward. The magnitude of this asymmetry between the two hemispheres is a key factor of the climate system through the displacement of tropical precipitation north of the equator and its influence on Arctic temperature and sea-ice extent. These asymmetric influences on heat are however not well constrained by observations or models. We identify a robust link between the ocean heat asymmetry and the large-scale distribution in atmospheric oxygen, using both atmospheric and oceanic observations and a suite of models (oceanic, climate and inverse). Novel aircraft observations from the pole-to-pole HIPPO campaign reveal that the ocean northward heat transport necessary to explain the atmospheric oxygen distribution is in the upper range of previous estimates from hydrographic sections and atmospheric reanalyses. Finally, we evidence a strong link between the oceanic transports of heat and natural carbon. This supports the existence of a strong southward transport of natural carbon at the global scale, a feature present at pre-industrial times and still underlying the anthropogenic signal today. We find that current climate models systematically underestimate these natural large-scale ocean meridional transports of heat and carbon, which bears on future climate projections, in particular concerning Arctic climate, possible shifts in rainfall and carbon sinks partition between the land and the ocean.

Assimilation of TOPEX/Poseidon altimeter data into a global ocean circulation model: How good are the results?

NASA Astrophysics Data System (ADS)

Fukumori, Ichiro; Raghunath, Ramanujam; Fu, Lee-Lueng; Chao, Yi

1999-11-01

The feasibility of assimilating satellite altimetry data into a global ocean general circulation model is studied. Three years of TOPEX/Poseidon data are analyzed using a global, three-dimensional, nonlinear primitive equation model. The assimilation's success is examined by analyzing its consistency and reliability measured by formal error estimates with respect to independent measurements. Improvements in model solution are demonstrated, in particular, properties not directly measured. Comparisons are performed with sea level measured by tide gauges, subsurface temperatures and currents from moorings, and bottom pressure measurements. Model representation errors dictate what can and cannot be resolved by assimilation, and its identification is emphasized.

Crystal settling and crystal growth caused by Ostwald Ripening in a terrestrial magma ocean under rotation

NASA Astrophysics Data System (ADS)

Maas, C.; Moeller, A.; Hansen, U.

2013-12-01

About 4.5 billion years ago the earth was covered by a heavily convecting and rotating global magma ocean which was caused by an impact of a mars-sized impactor in a later stage of the earth's accretion. After the separation of metal and silicate (see A. Möller, U. Hansen (2013)) and the formation of the earth's core it began to crystallize. Small silicate crystals emerge and grow by Ostwald Ripening when the fluid is supersaturated. This process results in shrinking of small crystals and growing of large crystals on behalf of the smaller ones. This leads to an altering of the crystal settling time. One question which is still under great debate is whether fractional or equilibrium crystallization occurred in the magma ocean. Fractional crystallization means that different mineral fractions settle one after the other which would lead to a strongly differentiated mantle after solidification of the magma ocean. In contrast to that equilibrium crystallization would result in a well mixed mantle. Whether fractional or equilibrium crystallization occurred is for example important for the starting model of plate tectonics or the understanding of the mantle development until today. To study the change of crystal radius in a convecting and rotating magma ocean we employed a 3D numerical model. Due to the low viscosity and strong rotation the influence of rotation on the early magma Ocean cannot be neglected. In the model the crystals are able to influence each other and the fluid flow. They are able to grow, shrink, vanish and form and gravitational, Coriolis and drag forces due to the fluid act on them. In our present work we study the crystal settling depending on different rotation rates and rotation axes with two configurations. For the polar setting the rotation axis is parallel, at the equator it is perpendicular to gravity. Low rotation at the pole leads to a large fraction of suspended crystals. With increasing rotation the crystals settle and form a thick layer at the bottom of the magma ocean. At the equator we find three regimes (see A. Möller, U. Hansen (2013)) depending on the rotation strength. At low rotation a high fraction of silicate crystals settle at the bottom. At higher rotation the crystals form a thick layer in the bottom 1/3 of box. At high rotation all crystals are suspended and we observe a ribbon structure in the middle of the box. With a second model we investigate growing and shrinking of crystals by Ostwald Ripening and include formation and melting. In general we observe the same behaviour and regimes as described above, however due to Ostwald Ripening the evolution of crystal radius with time depends on the strength of rotation and on the orientation of the rotation axis. Very first results show that at the pole the growth of the silicate crystals is limited. The resulting small radius leads to a slow crystal settling. At the equator the crystals are able to grow larger than at the pole and therefore settle faster. This could lead to an asymmetrical crystallization of the magma ocean. In an extreme case due to the different settling times this could lead to a well mixed mantle at the pole whereas at the equator the mantle could be strongly differentiated after the solidification of the magma ocean.

The effect of numerical methods on the simulation of mid-ocean ridge hydrothermal models

NASA Astrophysics Data System (ADS)

Carpio, J.; Braack, M.

2012-01-01

This work considers the effect of the numerical method on the simulation of a 2D model of hydrothermal systems located in the high-permeability axial plane of mid-ocean ridges. The behavior of hot plumes, formed in a porous medium between volcanic lava and the ocean floor, is very irregular due to convective instabilities. Therefore, we discuss and compare two different numerical methods for solving the mathematical model of this system. In concrete, we consider two ways to treat the temperature equation of the model: a semi-Lagrangian formulation of the advective terms in combination with a Galerkin finite element method for the parabolic part of the equations and a stabilized finite element scheme. Both methods are very robust and accurate. However, due to physical instabilities in the system at high Rayleigh number, the effect of the numerical method is significant with regard to the temperature distribution at a certain time instant. The good news is that relevant statistical quantities remain relatively stable and coincide for the two numerical schemes. The agreement is larger in the case of a mathematical model with constant water properties. In the case of a model with nonlinear dependence of the water properties on the temperature and pressure, the agreement in the statistics is clearly less pronounced. Hence, the presented work accentuates the need for a strengthened validation of the compatibility between numerical scheme (accuracy/resolution) and complex (realistic/nonlinear) models.

Dynamics of Andaman Sea circulation and its role in connecting the equatorial Indian Ocean to the Bay of Bengal

NASA Astrophysics Data System (ADS)

Chatterjee, Abhisek; Shankar, D.; McCreary, J. P.; Vinayachandran, P. N.; Mukherjee, A.

2017-04-01

Circulation in the Bay of Bengal (BoB) is driven not only by local winds, but are also strongly forced by the reflection of equatorial Kelvin waves (EKWs) from the eastern boundary of the Indian Ocean. The equatorial influence attains its peak during the monsoon-transition period when strong eastward currents force the strong EKWs along the equator. The Andaman Sea, lying between the Andaman and Nicobar island chains to its west and Indonesia, Thailand, and Myanmar to the south, east, and north, is connected to the equatorial ocean and the BoB by three primary passages, the southern (6°N), middle (10°N), and northern (15°N) channels. We use ocean circulation models, together with satellite altimeter data, to study the pathways by which equatorial signals pass through the Andaman Sea to the BoB and associated dynamical interactions in the process. The mean coastal circulation within the Andaman Sea and around the islands is primarily driven by equatorial forcing, with the local winds forcing a weak sea-level signal. On the other hand, the current forced by local winds is comparable to that forced remotely from the equator. Our results suggest that the Andaman and Nicobar Islands not only influence the circulation within the Andaman Sea, but also significantly alter the circulation in the interior bay and along the east coast of India, implying that they need to be represented accurately in numerical models of the Indian Ocean.

Motionally-induced electromagnetic fields generated by idealized ocean currents

NASA Astrophysics Data System (ADS)

Tyler, R. H.; Mysak, L. A.

Using the induction equation, we investigate the generation of electromagnetic fields by the motional electromagnetic induction due to ocean currents. In this paper, solutions are presented for a linear induction equation for the magnetic flux density vector which contains prescribed time-independent ocean current and conductivity fields. Once the magnetic flux density is known, the electric field and electric current density are easily obtained by differentiation. Solutions are given for several examples of idealized flow which include: 1) Vertically and horizontally sheared plane-parallel flow with depth-dependent conductivity; 2) A simple Stommel circulation gyre; and 3) Symmetric gyres. The results indicate that typical ocean current features induce magnetic fields with magnitudes reaching 100's of nT within the water and about 1-10 outside of the water. For the case of a field of gyres, the ocean-induced magnetic fields decay away from the ocean on spatial scales set by the horizontal scale of the ocean feature. At the altitudes of magnetic field satellite surveys, ocean-induced magnetic fields may retain values of a few nT, which are strong enough to be detected. Thus it is concluded that satellite observations of the earth's main magnetic field and, in particular, the observed temporal variations, could be affected by the ocean circulation. Summary and discussion In Section 3, we found exact solutions to the induction equation for idealized flows. The results gave magnitudes of about tens to hundreds of nT for the magnetic fields bH, about 10-5 V/m for the electric fields E, and about 10-5 A/m2 for the electric current density J induced by the ocean currents. These figures are in general agreement with the calculations of Lilley et al. (1993). In Section 4.2 we obtained solutions for the magnetic field above the ocean surface for the case of a Stommel gyre and a field of symmetric gyres. It was found in the last case that ocean gyres with a total transport of 100 Sv would have field magnitudes above the ocean of up to 23 nT and more typically between about 1 and 10 nT. The results also indicate that the decay scales for the magnetic field away from the ocean are of the same order as the horizontal scale of the flow. We calculated that flow features with scales of 100 km or more may retain magnitudes that are strong enough (a few nT) to be detected at satellite altitudes (300-500 km). Flows of smaller scale, however, would probably not be detected by current satellites. We have not explicitly solved for the magnetic fields that would be observed at magnetic observations on land since we have only treated cases of horizontally homogeneous conductivity. However, we conjecture that inland magnetic observations will also be affected by the ocean induction. The spatial decay inland from the ocean would again be set by the horizontal scale of the flow except when the electric currents involved in the induction are close to the coast. In this case there is a coastal effect, with the stronger fields decaying over a shorter scale. We have presented arguments which indicate that for uniform ocean conductivity s over a sediment layer of low conductivity, barotropic currents are efficient generators of electric fields but poor generators of electrical current and magnetic fields, while baroclinic currents are efficient generators of electrical current and magnetic fields, and (in our simple examples) poor generators of electric fields. When, however, s varies in the vertical, it appears that virtually all realistic forms of ocean circulation will be reasonably efficient generators of electrical current and magnetic fields. It is conceivable that the geomagnetic record from land and satellite observatories has captured oceanic signals. Another task before ocean and geomagnetic records can be linked, however, is to isolate the oceanic signal from the records which are known to be swarmed with magnetic signals from many other sources. How can we know that measured magnetic variations are due to variations in the ocean induction and not due to sources in the ionosphere or earth's core? This is a difficult problem, but there may be some ways to resolve it. Variations in ocean circulation or conductivity are rather slow compared to the rapid magnetic storms and most other variations due to external sources. Also the external effects at the earth's surface tend to have large spatial scales which allows removal using techniques such as 'remote referencing' as done by Lilley et al. (1993). With regard to sources in the earth's core, the geometric and electromagnetic filtering by the mantle are thought to prevent all but the lowest frequencies from reaching the earth's surface. Hence, it is conceivable that at the earth's surface, magnetic fields due to the earth's core can only appear as relatively smooth, slowly-varying fields with periods of decadal scale and longer. Hence, there is probably a fortuitous 'spectral window' through which we can view interannual variations in the ocean-induced fields. It is also important that the accuracy in measuring the time rate of change of the magnetic field on these time scales is greater than the accuracy of the field values (at least at the land observatories). This is because when differencing the magnetic series, errors in the baseline drift are reduced. Hence, it is probable that fluctuations in the ocean-induced magnetic fields would be easier to detect than the steady-state fields. The results presented here should also be helpful in designing future strategies for numerically modelling the ocean-induced electromagnetic fields. As we mentioned, (73) is similar in principle to the two-dimensional equation solved by Stephenson and Bryan (1992) but is more general, allowing for a more realistic description of conductivity and allowing for horizontal divergence in the conductivity transport. Other results in this paper may be helpful in finding new approaches to the problem of numerically modelling the oceanic induction. The flux form of the induction equation (18) could be a more convenient form to be used in established numerical algorithms using flux conservation. Also, (11) or (18) could be used with the Stokes or Divergence theorems to create a box-model description of global oceanic induction. When the scalar fluid property Λ described in Section 4.2.3 is taken to be a function of the conductivity, equation (76) is greatly simplified, suggesting that a three-dimensional numerical model using constant conductivity layers [similar in essence to the constant-density layer approach used in the OPYC ocean-circulation model written by Joseph Oberhuber (1993)] may be profitably exploited.

Tsunami Simulation using CIP Method with Characteristic Curve Equations and TVD-MacCormack Method

NASA Astrophysics Data System (ADS)

Fukazawa, Souki; Tosaka, Hiroyuki

2015-04-01

After entering 21st century, we already had two big tsunami disasters associated with Mw9 earthquakes in Sumatra and Japan. To mitigate the damages of tsunami, the numerical simulation technology combined with information technologies could provide reliable predictions in planning countermeasures to prevent the damage to the social system, making safety maps, and submitting early evacuation information to the residents. Shallow water equations are still solved not only for global scale simulation of the ocean tsunami propagation but also for local scale simulation of overland inundation in many tsunami simulators though three-dimensional model starts to be used due to improvement of CPU. One-dimensional shallow water equations are below: partial bm{Q}/partial t+partial bm{E}/partial x=bm{S} in which bm{Q}=( D M )), bm{E}=( M M^2/D+gD^2/2 )), bm{S}=( 0 -gDpartial z/partial x-gn2 M|M| /D7/3 )). where D[m] is total water depth; M[m^2/s] is water flux; z[m] is topography; g[m/s^2] is the gravitational acceleration; n[s/m1/3] is Manning's roughness coefficient. To solve these, the staggered leapfrog scheme is used in a lot of wide-scale tsunami simulator. But this scheme has a problem that lagging phase error occurs when courant number is small. In some practical simulation, a kind of diffusion term is added. In this study, we developed two wide-scale tsunami simulators with different schemes and compared usual scheme and other schemes in practicability and validity. One is a total variation diminishing modification of the MacCormack method (TVD-MacCormack method) which is famous for the simulation of compressible fluids. The other is the Cubic Interpolated Profile (CIP) method with characteristic curve equations transformed from shallow water equations. Characteristic curve equations derived from shallow water equations are below: partial R_x±/partial t+C_x±partial R_x±/partial x=∓ g/2partial z/partial x in which R_x±=√{gD}± u/2, C_x±=u± √{gD}. where u[m/s] is water velocity. It is difficult to solve the inundation on the land with these methods though These two methods are applicable to the ocean tsunami propagation. We studied how to apply these methods to overland inundation and how to couple the ocean global model with the land local model. Simple case studies of ocean tsunami propagation and overland tsunami inundation were performed to validate three methods comparing the results with theoretical solution. Finally, we performed case studies of the Great East Japan Earthquake in 2011 and confirmed the applicability to the actual tsunami.

Simulation of heat storages and associated heat budgets in the Pacific Ocean: 2. Interdecadal timescale

NASA Astrophysics Data System (ADS)

Auad, Guillermo; Miller, Arthur J.; White, Warren B.

1998-11-01

We use a primitive equation isopycnal model of the Pacific Ocean to simulate and diagnose the anomalous heat balance on interdecadal timescales associated with heat storage changes observed from 1970-1988 in the expendable bathythermograph (XBT) data set. Given the smallness of the interdecadal signals compared to the El Niño-Southern Oscillation (ENSO) signal, the agreement between model and observations is remarkably good. The total anomalous heat balance is made up of two parts, the diabatic part (from the model temperature equation) and the adiabatic part (from the model mass conservation equation) due to thermocline heave. We therefore describe our analysis of both the total and diabatic anomalous heat balances in four areas of the tropical and subtropical North Pacific Ocean in the upper 400 m. The interdecadal total (diabatic plus adiabatic) heat balance in the North Pacific Ocean is characterized by a complicated interplay of different physical processes, especially revealed in basin-scale averages of the heat budget components that have comparable amounts of variance. In smaller subregions, simpler balances hold. For example, in the western equatorial Pacific (area 1) the total heat content tendency term is nearly zero, so that a simple balance exists between surface heat flux, vertical heat transport, and horizontal mixing. In the western subtropical Pacific the total heat content tendency balances the three-dimensional divergence of the heat flux. We speculate that this complexity is indicative of multiple physical mechanisms involved in the generation of North Pacific interdecadal variability. The diabatic heat balance north of 24°N, a region of special interest to The World Ocean Circulation Experiment (WOCE), can be simplified to a balance between the tendency term, surface heat flux, and meridional advection, the last term dominated by anomalous advection of mean temperature gradients. For the western equatorial region the diabatic heat content tendency is nearly zero and the steady balance involves simply horizontal advection and the surface heat flux, which at these latitudes has a damping role in the model. An important finding of this study is the identification of two interdecadal timescales, roughly 10 and 20 years, both similar to those reported by other investigators in recent years. [Tourre et al., 1998; Latif and Barnett, 1994; Robertson, 1995; White et al, 1997; Gu and Philander, 1997; Jacobs et al., 1994]. The 20-year timescale is only present in diabatic heat budget components, while the 10-year timescale is present in both diabatic and adiabatic components. The 10-year timescale can also be seen in the surface heat flux time series, but it occurs in the ocean adiabatic components which demonstrates the importance of oceanic adjustment through Rossby wave dynamics on decadal timescales.

Quantification of A Tropical Missing Source From Ocean For The Carbonyl Sulfide Global Budget

NASA Astrophysics Data System (ADS)

Kuai, Le; Worden, John; Campbell, Elliott; Kulawik, Susan; Lee, Meemong; Montzka, Stephen; Berry, Joe; Baker, Ian; Denning, Scott; Kawa, Randy; Bian, Huisheng; Yung, Yuk

2015-04-01

Quantifying the carbonyl sulfide (OCS) surface fluxes contributes to the understanding of both sulfur cycle and carbon cycle. Although the major sources and sinks of OCS are well recognized, the uncertainties of individual types of the fluxes remain large. With the understanding of a large underestimate of ecosystem uptake, it suggests a large missing ocean source over tropical region to compensate the increased sink. However before AURA Tropospheric Emissions Spectrometer (TES) OCS data is released, no direct measurements have been taken to test this hypothesis. In this study, we performed a flux inversion to update the fluxes from TES OCS. Then we compared three experimental GEOS-Chem forward model runs driven by different fluxes based on TES inversion to HIPPO aircraft estimates in free troposphere and also to NOAA near surface observations. The TES data supports the hypothesis that a large source from tropical ocean is missing in the current OCS global budget and suggests that the source is even larger than that proposed in Berry et al., (2013). Consequently, it leads to a larger land uptake and increase the estimates of GPP. TES data also suggests the missing oceanic source is not symmetric about equator. It is strong and distributed further north of the equator (to 40°N) but weak and narrow south of the equator (to 20°S).

Are Hydrostatic Models Still Capable of Simulating Oceanic Fronts

DTIC Science & Technology

2016-11-10

Coriolis effect is added to the model momentum equations...nonhydrostatic (NH) models to address the relevance of NH effects on the evolution of density fronts and the development of meso- and submeso-scale vertical...nonhydrostatic (NH) models to address the relevance of NH effects on the evolution of density fronts and the development of meso- and submeso-scale vertical

Dynamics of the global meridional ice flow of Europa's icy shell

NASA Astrophysics Data System (ADS)

Ashkenazy, Yosef; Sayag, Roiy; Tziperman, Eli

2018-01-01

Europa is one of the most probable places in the solar system to find extra-terrestrial life1,2, motivating the study of its deep ( 100 km) ocean3-6 and thick icy shell3,7-11. The chaotic terrain patterns on Europa's surface12-15 have been associated with vertical convective motions within the ice8,10. Horizontal gradients of ice thickness16,17 are expected due to the large equator-to-pole gradient of surface temperature and can drive a global horizontal ice flow, yet such a flow and its observable implications have not been studied. We present a global ice flow model for Europa composed of warm, soft ice flowing beneath a cold brittle rigid ice crust3. The model is coupled to an underlying (diffusive) ocean and includes the effect of tidal heating and convection within the ice. We show that Europa's ice can flow meridionally due to pressure gradients associated with equator-to-pole ice thickness differences, which can be up to a few km and can be reduced both by ice flow and due to ocean heat transport. The ice thickness and meridional flow direction depend on whether the ice convects or not; multiple (convecting and non-convecting) equilibria are found. Measurements of the ice thickness and surface temperature from future Europa missions18,19 can be used with our model to deduce whether Europa's icy shell convects and to constrain the effectiveness of ocean heat transport.

Role of Compressibility on Tsunami Propagation

NASA Astrophysics Data System (ADS)

Abdolali, Ali; Kirby, James T.

2017-12-01

In the present paper, we aim to reduce the discrepancies between tsunami arrival times evaluated from tsunami models and real measurements considering the role of ocean compressibility. We perform qualitative studies to reveal the phase speed reduction rate via a modified version of the Mild Slope Equation for Weakly Compressible fluid (MSEWC) proposed by Sammarco et al. (2013). The model is validated against a 3-D computational model. Physical properties of surface gravity waves are studied and compared with those for waves evaluated from an incompressible flow solver over realistic geometry for 2011 Tohoku-oki event, revealing reduction in phase speed.