A Three-Dimensional Variational Data Assimilation Scheme for the Regional Ocean Modeling System: Implementation and Basic Experiments

NASA Technical Reports Server (NTRS)

Li, Zhijin; Chao, Yi; McWilliams, James C.; Ide, Kayo

2008-01-01

A three-dimensional variational data assimilation scheme for the Regional Ocean Modeling System (ROMS), named ROMS3DVAR, has been described in the work of Li et al. (2008). In this paper, ROMS3DVAR is applied to the central California coastal region, an area characterized by inhomogeneity and anisotropy, as well as by dynamically unbalanced flows. A method for estimating the model error variances from limited observations is presented, and the construction of the inhomogeneous and anisotropic error correlations based on the Kronecker product is demonstrated. A set of single observation experiments illustrates the inhomogeneous and anisotropic error correlations and weak dynamic constraints used. Results are presented from the assimilation of data gathered during the Autonomous Ocean Sampling Network (AOSN) experiment during August 2003. The results show that ROMS3DVAR is capable of reproducing complex flows associated with upwelling and relaxation, as well as the rapid transitions between them. Some difficulties encountered during the experiment are also discussed.

Local stability of a five dimensional food chain model in the ocean

NASA Astrophysics Data System (ADS)

Kusumawinahyu, W. M.; Hidayatulloh, M. R.

2014-02-01

This paper discuss a food chain model on a microbiology ecosystem in the ocean, where predation process occurs. Four population growth rates are discussed, namely bacteria, phytoplankton, zooplankton, and protozoa growth rate. When the growth of nutrient density is also considered, the model is governed by a five dimensional dynamical system. The system considered in this paper is a modification of a model proposed by Hadley and Forbes [1], by taking Holling Type I as the functional response. For sake of simplicity, the model needs to be scaled. Dynamical behavior, such as existence condition of equilibrium points and their local stability are addressed. There are eight equilibrium points, where two of them exist under certain conditions. Three equilibrium points are unstable, while two points stable under certain conditions and the other three points are stable if the Ruth-Hurwitz criteria are satisfied. Numerical simulations are carried out to illustrate analytical findings.

Nutrient Distributions, Transports, and Budgets on the Inner Margin of a River-Dominated Continental Shelf

DTIC Science & Technology

2013-10-02

and budgets on the inner margin of a river-dominated continental shelf, J. Geophys. Res. Oceans , 118, 4822–4838, doi:10.1002/jgrc.20362. 1...13/10.1002/jgrc.20362 4822 JOURNAL OF GEOPHYSICAL RESEARCH: OCEANS , VOL. 118, 4822–4838, doi:10.1002/jgrc.20362, 2013 Report Documentation Page Form...shelf, and current velocities obtained from a three-dimensional (3-D) hydro- dynamic model (the Navy Coastal Ocean Model). The budget terms were used to

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

Moment tensor inversion of the 2016 southeast offshore Mie earthquake in the Tonankai region using a three-dimensional velocity structure model: effects of the accretionary prism and subducting oceanic plate

NASA Astrophysics Data System (ADS)

Takemura, Shunsuke; Kimura, Takeshi; Saito, Tatsuhiko; Kubo, Hisahiko; Shiomi, Katsuhiko

2018-03-01

The southeast offshore Mie earthquake occurred on April 1, 2016 near the rupture area of the 1944 Tonankai earthquake, where seismicity around the interface of the Philippine Sea plate had been very low until this earthquake. Since this earthquake occurred outside of seismic arrays, the focal mechanism and depth were not precisely constrained using a one-dimensional velocity model, as in a conventional approach. We conducted a moment tensor inversion of this earthquake by using a three-dimensional velocity structure model. Before the analysis of observed data, we investigated the effects of offshore heterogeneous structures such as the seawater, accretionary prism, and subducting oceanic plate by using synthetic seismograms in a full three-dimensional model and simpler models. The accretionary prism and subducting oceanic plate play important roles in the moment tensor inversion for offshore earthquakes in the subduction zone. Particularly, the accretionary prism, which controls the excitation and propagation of long-period surface waves around the offshore region, provides better estimations of the centroid depths and focal mechanisms of earthquakes around the Nankai subduction zone. The result of moment tensor inversion for the 2016 southeast offshore Mie earthquake revealed low-angle thrust faulting with a moment magnitude of 5.6. According to geophysical surveys in the Nankai Trough, our results suggest that the rupture of this earthquake occurred on the interface of the Philippine Sea plate, rather than on a mega-splay fault. Detailed comparisons of first-motion polarizations provided additional constraints of the rupture that occurred on the interface of the Philippine Sea plate.

Linear stability analysis of the three-dimensional thermally-driven ocean circulation: application to interdecadal oscillations

NASA Astrophysics Data System (ADS)

Huck, Thierry; Vallis, Geoffrey K.

2001-08-01

What can we learn from performing a linear stability analysis of the large-scale ocean circulation? Can we predict from the basic state the occurrence of interdecadal oscillations, such as might be found in a forward integration of the full equations of motion? If so, do the structure and period of the linearly unstable modes resemble those found in a forward integration? We pursue here a preliminary study of these questions for a case in idealized geometry, in which the full nonlinear behavior can also be explored through forward integrations. Specifically, we perform a three-dimensional linear stability analysis of the thermally-driven circulation of the planetary geostrophic equations. We examine the resulting eigenvalues and eigenfunctions, comparing them with the structure of the interdecadal oscillations found in the fully nonlinear model in various parameter regimes. We obtain a steady state by running the time-dependent, nonlinear model to equilibrium using restoring boundary conditions on surface temperature. If the surface heat fluxes are then diagnosed, and these values applied as constant flux boundary conditions, the nonlinear model switches into a state of perpetual, finite amplitude, interdecadal oscillations. We construct a linearized version of the model by empirically evaluating the tangent linear matrix at the steady state, under both restoring and constant-flux boundary conditions. An eigen-analysis shows there are no unstable eigenmodes of the linearized model with restoring conditions. In contrast, under constant flux conditions, we find a single unstable eigenmode that shows a striking resemblance to the fully-developed oscillations in terms of three-dimensional structure, period and growth rate. The mode may be damped through either surface restoring boundary conditions or sufficiently large horizontal tracer diffusion. The success of this simple numerical method in idealized geometry suggests applications in the study of the stability of the ocean circulation in more realistic configurations, and the possibility of predicting potential oceanic modes, even weakly damped, that might be excited by stochastic atmospheric forcing or mesoscale ocean eddies.

Quasi-Geostrophic Diagnosis of Mixed-Layer Dynamics Embedded in a Mesoscale Turbulent Field

NASA Astrophysics Data System (ADS)

Chavanne, C. P.; Klein, P.

2016-02-01

A new quasi-geostrophic model has been developed to diagnose the three-dimensional circulation, including the vertical velocity, in the upper ocean from high-resolution observations of sea surface height and buoyancy. The formulation for the adiabatic component departs from the classical surface quasi-geostrophic framework considered before since it takes into account the stratification within the surface mixed-layer that is usually much weaker than that in the ocean interior. To achieve this, the model approximates the ocean with two constant-stratification layers : a finite-thickness surface layer (or the mixed-layer) and an infinitely-deep interior layer. It is shown that the leading-order adiabatic circulation is entirely determined if both the surface streamfunction and buoyancy anomalies are considered. The surface layer further includes a diabatic dynamical contribution. Parameterization of diabatic vertical velocities is based on their restoring impacts of the thermal-wind balance that is perturbed by turbulent vertical mixing of momentum and buoyancy. The model skill in reproducing the three-dimensional circulation in the upper ocean from surface data is checked against the output of a high-resolution primitive-equation numerical simulation. Correlation between simulated and diagnosed vertical velocities are significantly improved in the mixed-layer for the new model compared to the classical surface quasi-geostrophic model, reaching 0.9 near the surface.

Multi-Scale Three-Dimensional Variational Data Assimilation System for Coastal Ocean Prediction

NASA Technical Reports Server (NTRS)

Li, Zhijin; Chao, Yi; Li, P. Peggy

2012-01-01

A multi-scale three-dimensional variational data assimilation system (MS-3DVAR) has been formulated and the associated software system has been developed for improving high-resolution coastal ocean prediction. This system helps improve coastal ocean prediction skill, and has been used in support of operational coastal ocean forecasting systems and field experiments. The system has been developed to improve the capability of data assimilation for assimilating, simultaneously and effectively, sparse vertical profiles and high-resolution remote sensing surface measurements into coastal ocean models, as well as constraining model biases. In this system, the cost function is decomposed into two separate units for the large- and small-scale components, respectively. As such, data assimilation is implemented sequentially from large to small scales, the background error covariance is constructed to be scale-dependent, and a scale-dependent dynamic balance is incorporated. This scheme then allows effective constraining large scales and model bias through assimilating sparse vertical profiles, and small scales through assimilating high-resolution surface measurements. This MS-3DVAR enhances the capability of the traditional 3DVAR for assimilating highly heterogeneously distributed observations, such as along-track satellite altimetry data, and particularly maximizing the extraction of information from limited numbers of vertical profile observations.

Partial spline models for the inclusion of tropopause and frontal boundary information in otherwise smooth two- and three-dimensional objective analysis

NASA Technical Reports Server (NTRS)

Shiau, Jyh-Jen; Wahba, Grace; Johnson, Donald R.

1986-01-01

A new method, based on partial spline models, is developed for including specified discontinuities in otherwise smooth two- and three-dimensional objective analyses. The method is appropriate for including tropopause height information in two- and three-dimensinal temperature analyses, using the O'Sullivan-Wahba physical variational method for analysis of satellite radiance data, and may in principle be used in a combined variational analysis of observed, forecast, and climate information. A numerical method for its implementation is described and a prototype two-dimensional analysis based on simulated radiosonde and tropopause height data is shown. The method may also be appropriate for other geophysical problems, such as modeling the ocean thermocline, fronts, discontinuities, etc.

The impacts of oceanic deep temperature perturbations in the North Atlantic on decadal climate variability and predictability

NASA Astrophysics Data System (ADS)

Germe, Agathe; Sévellec, Florian; Mignot, Juliette; Fedorov, Alexey; Nguyen, Sébastien; Swingedouw, Didier

2017-12-01

Decadal climate predictability in the North Atlantic is largely related to ocean low frequency variability, whose sensitivity to initial conditions is not very well understood. Recently, three-dimensional oceanic temperature anomalies optimally perturbing the North Atlantic Mean Temperature (NAMT) have been computed via an optimization procedure using a linear adjoint to a realistic ocean general circulation model. The spatial pattern of the identified perturbations, localized in the North Atlantic, has the largest magnitude between 1000 and 4000 m depth. In the present study, the impacts of these perturbations on NAMT, on the Atlantic meridional overturning circulation (AMOC), and on climate in general are investigated in a global coupled model that uses the same ocean model as was used to compute the three-dimensional optimal perturbations. In the coupled model, these perturbations induce AMOC and NAMT anomalies peaking after 5 and 10 years, respectively, generally consistent with the ocean-only linear predictions. To further understand their impact, their magnitude was varied in a broad range. For initial perturbations with a magnitude comparable to the internal variability of the coupled model, the model response exhibits a strong signature in sea surface temperature and precipitation over North America and the Sahel region. The existence and impacts of these ocean perturbations have important implications for decadal prediction: they can be seen either as a source of predictability or uncertainty, depending on whether the current observing system can detect them or not. In fact, comparing the magnitude of the imposed perturbations with the uncertainty of available ocean observations such as Argo data or ocean state estimates suggests that only the largest perturbations used in this study could be detectable. This highlights the importance for decadal climate prediction of accurate ocean density initialisation in the North Atlantic at intermediate and greater depths.

NUMERICAL MODELS AS ENABLING TOOLS FOR TIDAL-STREAM ENERGY EXTRACTION AND ENVIRONMENTAL IMPACT ASSESSMENT

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

Yang, Zhaoqing; Wang, Taiping

This paper presents a modeling study conducted to evaluate tidal-stream energy extraction and its associated potential environmental impacts using a three-dimensional unstructured-grid coastal ocean model, which was coupled with a water-quality model and a tidal-turbine module.

Assimilation of seasonal chlorophyll and nutrient data into an adjoint three-dimensional ocean carbon cycle model: Sensitivity analysis and ecosystem parameter optimization

NASA Astrophysics Data System (ADS)

Tjiputra, Jerry F.; Polzin, Dierk; Winguth, Arne M. E.

2007-03-01

An adjoint method is applied to a three-dimensional global ocean biogeochemical cycle model to optimize the ecosystem parameters on the basis of SeaWiFS surface chlorophyll observation. We showed with identical twin experiments that the model simulated chlorophyll concentration is sensitive to perturbation of phytoplankton and zooplankton exudation, herbivore egestion as fecal pellets, zooplankton grazing, and the assimilation efficiency parameters. The assimilation of SeaWiFS chlorophyll data significantly improved the prediction of chlorophyll concentration, especially in the high-latitude regions. Experiments that considered regional variations of parameters yielded a high seasonal variance of ecosystem parameters in the high latitudes, but a low variance in the tropical regions. These experiments indicate that the adjoint model is, despite the many uncertainties, generally capable to optimize sensitive parameters and carbon fluxes in the euphotic zone. The best fit regional parameters predict a global net primary production of 36 Pg C yr-1, which lies within the range suggested by Antoine et al. (1996). Additional constraints of nutrient data from the World Ocean Atlas showed further reduction in the model-data misfit and that assimilation with extensive data sets is necessary.

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

Three-Dimensional Shallow Water Adaptive Hydraulics (ADH-SW3) Validation: Galveston Bay Hydrodynamics and Salinity Transport

DTIC Science & Technology

2015-04-01

model mesh with elements (vertical co-ordinate in meters). ....................... 5 Figure 3. Ocean tidal boundary (Hour 0 = 1 Jan 1990, 12:00 a.m...7 Figure 4. Ocean salt boundary (Hour 0 = 1 Jan 1990, 12:00 a.m...channel and the connections of Galveston Bay to the open ocean . Figures 1 and 2 illustrate the distribution of vertical layers and resolution in the

A reanalysis dataset of the South China Sea.

PubMed

Zeng, Xuezhi; Peng, Shiqiu; Li, Zhijin; Qi, Yiquan; Chen, Rongyu

2014-01-01

Ocean reanalysis provides a temporally continuous and spatially gridded four-dimensional estimate of the ocean state for a better understanding of the ocean dynamics and its spatial/temporal variability. Here we present a 19-year (1992-2010) high-resolution ocean reanalysis dataset of the upper ocean in the South China Sea (SCS) produced from an ocean data assimilation system. A wide variety of observations, including in-situ temperature/salinity profiles, ship-measured and satellite-derived sea surface temperatures, and sea surface height anomalies from satellite altimetry, are assimilated into the outputs of an ocean general circulation model using a multi-scale incremental three-dimensional variational data assimilation scheme, yielding a daily high-resolution reanalysis dataset of the SCS. Comparisons between the reanalysis and independent observations support the reliability of the dataset. The presented dataset provides the research community of the SCS an important data source for studying the thermodynamic processes of the ocean circulation and meso-scale features in the SCS, including their spatial and temporal variability.

A reanalysis dataset of the South China Sea

PubMed Central

Zeng, Xuezhi; Peng, Shiqiu; Li, Zhijin; Qi, Yiquan; Chen, Rongyu

2014-01-01

Ocean reanalysis provides a temporally continuous and spatially gridded four-dimensional estimate of the ocean state for a better understanding of the ocean dynamics and its spatial/temporal variability. Here we present a 19-year (1992–2010) high-resolution ocean reanalysis dataset of the upper ocean in the South China Sea (SCS) produced from an ocean data assimilation system. A wide variety of observations, including in-situ temperature/salinity profiles, ship-measured and satellite-derived sea surface temperatures, and sea surface height anomalies from satellite altimetry, are assimilated into the outputs of an ocean general circulation model using a multi-scale incremental three-dimensional variational data assimilation scheme, yielding a daily high-resolution reanalysis dataset of the SCS. Comparisons between the reanalysis and independent observations support the reliability of the dataset. The presented dataset provides the research community of the SCS an important data source for studying the thermodynamic processes of the ocean circulation and meso-scale features in the SCS, including their spatial and temporal variability. PMID:25977803

Three-dimensional S-wave tomography under Axial Seamount

NASA Astrophysics Data System (ADS)

Baillard, C.; Wilcock, W. S. D.; Arnulf, A. F.; Tolstoy, M.; Waldhauser, F.

2017-12-01

Axial Seamount is a submarine volcano located at the intersection of the Juande Fuca Ridge and the Cobb-Eickelberg hotspot 500 km off the coast of thenorthwestern United States. The seamount, which rises 1 km above the seafloor, ischaracterized by a shallow caldera that is elongated in the N-S direction, measure 8km by 3 km and sits on top of a 14 km by 3 km magma reservoir. Two eruptive eventsin 1998 and 2011 motivated the deployment in 2014 of a real time cabled observatorywithin the Axial caldera, as part of the Ocean Observatories Initiative (OOI).Theobservatory includes a network of seven seismometers that span the southern half ofthe caldera. Five months after the observatory came on-line in November 2014, thevolcano erupted on April 24, 2015. Well over 100,000 events were located in thevicinity of the caldera, delineating an outward dipping ring fault that extends fromnear the surface to the magma body at 2 km depth and which accommodatesinflation and deflation of the volcano.The initial earthquake locations have beenobtained with a one-dimensional velocity model but the travel time residuals suggeststrong heterogeneities. A three-dimensional P-wave velocity model, obtained bycombining multichannel and ocean bottom seismometer refraction data, is being usedto refine locations but the three-dimensional S-wave structure is presently unknown.In most mid-ocean ridge settings, the distribution of earthquakes is not conducive forjoint inversions for S-wave velocity and hypocentral parameters because there are fewcrossing ray paths but at Axial the presence of a ring fault that is seismically active atall depths on both the east and west side of the caldera, provides a reasonablegeometry for such efforts. We will present the results of joint inversions that assumethe existing three-dimensional P wave velocity model and solve for VP/VS structure andhypocentral parameters using LOTOS, an algorithm that solves the forward problemusing ray bending.The resulting model of S-wave velocities will provide newconstraints on the volcanic structure of the caldera, the distribution and characteristicsof fractures, and the effects of hydrothermal circulation. The model will also lead toimproved earthquakes locations that are critical for a fine scale interpretation of thefault system.

Three-dimensional vector modeling and restoration of flat finite wave tank radiometric measurements

NASA Technical Reports Server (NTRS)

Truman, W. M.; Balanis, C. A.

1977-01-01

The three-dimensional vector interaction between a microwave radiometer and a wave tank was modeled. Computer programs for predicting the response of the radiometer to the brightness temperature characteristics of the surroundings were developed along with a computer program that can invert (restore) the radiometer measurements. It is shown that the computer programs can be used to simulate the viewing of large bodies of water, and is applicable to radiometer measurements received from satellites monitoring the ocean. The water temperature, salinity, and wind speed can be determined.

Modeling of the Global Water Cycle - Analytical Models

Treesearch

Yongqiang Liu; Roni Avissar

2005-01-01

Both numerical and analytical models of coupled atmosphere and its underlying ground components (land, ocean, ice) are useful tools for modeling the global and regional water cycle. Unlike complex three-dimensional climate models, which need very large computing resources and involve a large number of complicated interactions often difficult to interpret, analytical...

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.

On the modeling of the 2010 Gulf of Mexico Oil Spill

NASA Astrophysics Data System (ADS)

Mariano, A. J.; Kourafalou, V. H.; Srinivasan, A.; Kang, H.; Halliwell, G. R.; Ryan, E. H.; Roffer, M.

2011-09-01

Two oil particle trajectory forecasting systems were developed and applied to the 2010 Deepwater Horizon Oil Spill in the Gulf of Mexico. Both systems use ocean current fields from high-resolution numerical ocean circulation model simulations, Lagrangian stochastic models to represent unresolved sub-grid scale variability to advect oil particles, and Monte Carlo-based schemes for representing uncertain biochemical and physical processes. The first system assumes two-dimensional particle motion at the ocean surface, the oil is in one state, and the particle removal is modeled as a Monte Carlo process parameterized by a one number removal rate. Oil particles are seeded using both initial conditions based on observations and particles released at the location of the Maconda well. The initial conditions (ICs) of oil particle location for the two-dimensional surface oil trajectory forecasts are based on a fusing of all available information including satellite-based analyses. The resulting oil map is digitized into a shape file within which a polygon filling software generates longitude and latitude with variable particle density depending on the amount of oil present in the observations for the IC. The more complex system assumes three (light, medium, heavy) states for the oil, each state has a different removal rate in the Monte Carlo process, three-dimensional particle motion, and a particle size-dependent oil mixing model. Simulations from the two-dimensional forecast system produced results that qualitatively agreed with the uncertain "truth" fields. These simulations validated the use of our Monte Carlo scheme for representing oil removal by evaporation and other weathering processes. Eulerian velocity fields for predicting particle motion from data-assimilative models produced better particle trajectory distributions than a free running model with no data assimilation. Monte Carlo simulations of the three-dimensional oil particle trajectory, whose ensembles were generated by perturbing the size of the oil particles and the fraction in a given size range that are released at depth, the two largest unknowns in this problem. 36 realizations of the model were run with only subsurface oil releases. An average of these results yields that after three months, about 25% of the oil remains in the water column and that most of the oil is below 800 m.

Dyamical Systems Theory and Lagrangian Data Assimilation in 4D Geophysical Fluid Dynamics

DTIC Science & Technology

The long-term goal of our project (known as OCEAN 3D +1) was to better understand and predict ocean circulation features that are fundamentally three...dimensional in space and that vary in time. In particular, we sought to quantify the dynamical processes that govern the formation , evolution, and...predictability of 3D +1 transport pathways in the ocean. Our approach was to develop algorithms to thoroughly analyze a hierarchy of model and

SOMAR-LES: A framework for multi-scale modeling of turbulent stratified oceanic flows

NASA Astrophysics Data System (ADS)

Chalamalla, Vamsi K.; Santilli, Edward; Scotti, Alberto; Jalali, Masoud; Sarkar, Sutanu

2017-12-01

A new multi-scale modeling technique, SOMAR-LES, is presented in this paper. Localized grid refinement gives SOMAR (the Stratified Ocean Model with Adaptive Resolution) access to small scales of the flow which are normally inaccessible to general circulation models (GCMs). SOMAR-LES drives a LES (Large Eddy Simulation) on SOMAR's finest grids, forced with large scale forcing from the coarser grids. Three-dimensional simulations of internal tide generation, propagation and scattering are performed to demonstrate this multi-scale modeling technique. In the case of internal tide generation at a two-dimensional bathymetry, SOMAR-LES is able to balance the baroclinic energy budget and accurately model turbulence losses at only 10% of the computational cost required by a non-adaptive solver running at SOMAR-LES's fine grid resolution. This relative cost is significantly reduced in situations with intermittent turbulence or where the location of the turbulence is not known a priori because SOMAR-LES does not require persistent, global, high resolution. To illustrate this point, we consider a three-dimensional bathymetry with grids adaptively refined along the tidally generated internal waves to capture remote mixing in regions of wave focusing. The computational cost in this case is found to be nearly 25 times smaller than that of a non-adaptive solver at comparable resolution. In the final test case, we consider the scattering of a mode-1 internal wave at an isolated two-dimensional and three-dimensional topography, and we compare the results with Legg (2014) numerical experiments. We find good agreement with theoretical estimates. SOMAR-LES is less dissipative than the closure scheme employed by Legg (2014) near the bathymetry. Depending on the flow configuration and resolution employed, a reduction of more than an order of magnitude in computational costs is expected, relative to traditional existing solvers.

Further Studies on Oceanic Biogeochemistry and Carbon Cycling

NASA Technical Reports Server (NTRS)

Signorini, S. R.; McClain, C. R.

2003-01-01

This TM consists of two chapters. Chapter I describes the development of a coupled, one-dimensional biogeochemical model using turbulence closure mixed layer (TCMLM) dynamics. The model is applied to the Sargasso Sea at the BATS (Bermuda Atlantic Time Series) site and the results are compared with a previous model study in the same region described in NASNTP-2001-209991. The use of the TCMLM contributed to some improvements in the model simulation of chlorophyll, PAR, nitrate, phosphate, and oxygen, but most importantly, the current model achieved good agreement with the data with much more realistic background eddy diffusivity. However, off-line calculations of horizontal transport of biogeochemical properties revealed that one-dimensional dynamics can only provide a limited assessment of the nutrient and carbon balances at BATS. Future studies in the BATS region will require comprehensive three-dimensional field studies, combined with three-dimensional eddy resolving numerical experiments, to adequately quantify the impact of the local and remote forcing on ecosystem dynamics and carbon cycling. Chapter II addresses the sensitivity of global sea-air CO, flux estimates to wind speed, temperature, and salinity. Sensitivity analyses of sea-air CO, flux to wind speed climatologies, gas transfer algorithms, SSS and SST were conducted for the global oceans and regional domains. Large uncertainties in the global sea-air flux are identified, primarily due to the different gas transfer algorithms used. The sensitivity of the sea-air flux to SST and SSS is similar in magnitude to the effect of using different wind climatologies. Globally, the mean ocean uptake of CO, changes by 5 to 16%, depending upon the combination of SST and SSS used.

Reduction of initial shock in decadal predictions using a new initialization strategy

NASA Astrophysics Data System (ADS)

He, Yujun; Wang, Bin; Liu, Mimi; Liu, Li; Yu, Yongqiang; Liu, Juanjuan; Li, Ruizhe; Zhang, Cheng; Xu, Shiming; Huang, Wenyu; Liu, Qun; Wang, Yong; Li, Feifei

2017-08-01

A novel full-field initialization strategy based on the dimension-reduced projection four-dimensional variational data assimilation (DRP-4DVar) is proposed to alleviate the well-known initial shock occurring in the early years of decadal predictions. It generates consistent initial conditions, which best fit the monthly mean oceanic analysis data along the coupled model trajectory in 1 month windows. Three indices to measure the initial shock intensity are also proposed. Results indicate that this method does reduce the initial shock in decadal predictions by Flexible Global Ocean-Atmosphere-Land System model, Grid-point version 2 (FGOALS-g2) compared with the three-dimensional variational data assimilation-based nudging full-field initialization for the same model and is comparable to or even better than the different initialization strategies for other fifth phase of the Coupled Model Intercomparison Project (CMIP5) models. Better hindcasts of global mean surface air temperature anomalies can be obtained than in other FGOALS-g2 experiments. Due to the good model response to external forcing and the reduction of initial shock, higher decadal prediction skill is achieved than in other CMIP5 models.

Resonance phenomena in a time-dependent, three-dimensional model of an idealized eddy

NASA Astrophysics Data System (ADS)

Rypina, I. I.; Pratt, L. J.; Wang, P.; Äe; -zgökmen, T. M.; Mezic, I.

2015-08-01

We analyze the geometry of Lagrangian motion and material barriers in a time-dependent, three-dimensional, Ekman-driven, rotating cylinder flow, which serves as an idealization for an isolated oceanic eddy and other overturning cells with cylindrical geometry in the ocean and atmosphere. The flow is forced at the top through an oscillating upper lid, and the response depends on the frequency and amplitude of lid oscillations. In particular, the Lagrangian geometry changes near the resonant tori of the unforced flow, whose frequencies are rationally related to the forcing frequencies. Multi-scale analytical expansions are used to simplify the flow in the vicinity of resonant trajectories and to investigate the resonant flow geometries. The resonance condition and scaling can be motivated by simple physical argument. The theoretically predicted flow geometries near resonant trajectories have then been confirmed through numerical simulations in a phenomenological model and in a full solution of the Navier-Stokes equations.

Marine infectious disease ecology

USGS Publications Warehouse

Lafferty, Kevin D.

2017-01-01

To put marine disease impacts in context requires a broad perspective on the roles infectious agents have in the ocean. Parasites infect most marine vertebrate and invertebrate species, and parasites and predators can have comparable biomass density, suggesting they play comparable parts as consumers in marine food webs. Although some parasites might increase with disturbance, most probably decline as food webs unravel. There are several ways to adapt epidemiological theory to the marine environment. In particular, because the ocean represents a three-dimensional moving habitat for hosts and parasites, models should open up the spatial scales at which infective stages and host larvae travel. In addition to open recruitment and dimensionality, marine parasites are subject to fishing, filter feeders, dosedependent infection, environmental forcing, and death-based transmission. Adding such considerations to marine disease models will make it easier to predict which infectious diseases will increase or decrease in a changing ocean.

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.

Broadband and polarization reflectors in the lookdown, Selene vomer

PubMed Central

Zhao, Shulei; Brady, Parrish Clawson; Gao, Meng; Etheredge, Robert Ian; Kattawar, George W.; Cummings, Molly E.

2015-01-01

Predator evasion in the open ocean is difficult because there are no objects to hide behind. The silvery surface of fish plays an important role in open water camouflage. Various models have been proposed to account for the broadband reflectance by the fish skin that involve one-dimensional variations in the arrangement of guanine crystal reflectors, yet the three-dimensional organization of these guanine platelets have not been well characterized. Here, we report the three-dimensional organization and the optical properties of integumentary guanine platelets in a silvery marine fish, the lookdown (Selene vomer). Our structural analysis and computational modelling show that stacks of guanine platelets with random yaw angles in the fish skin produce broadband reflectance via colour mixing. Optical axes of the guanine platelets and the collagen layer are aligned closely and provide bulk birefringence properties that influence the polarization reflectance by the skin. These data demonstrate how the lookdown preserves or alters polarization states at different incident polarization angles. These optical properties resulted from the organization of these guanine platelets and the collagen layer may have implications for open ocean camouflage in varying light fields. PMID:25673301

Landfalling Tropical Cyclones: Forecast Problems and Associated Research Opportunities

USGS Publications Warehouse

Marks, F.D.; Shay, L.K.; Barnes, G.; Black, P.; Demaria, M.; McCaul, B.; Mounari, J.; Montgomery, M.; Powell, M.; Smith, J.D.; Tuleya, B.; Tripoli, G.; Xie, Lingtian; Zehr, R.

1998-01-01

The Fifth Prospectus Development Team of the U.S. Weather Research Program was charged to identify and delineate emerging research opportunities relevant to the prediction of local weather, flooding, and coastal ocean currents associated with landfalling U.S. hurricanes specifically, and tropical cyclones in general. Central to this theme are basic and applied research topics, including rapid intensity change, initialization of and parameterization in dynamical models, coupling of atmospheric and oceanic models, quantitative use of satellite information, and mobile observing strategies to acquire observations to evaluate and validate predictive models. To improve the necessary understanding of physical processes and provide the initial conditions for realistic predictions, a focused, comprehensive mobile observing system in a translating storm-coordinate system is required. Given the development of proven instrumentation and improvement of existing systems, three-dimensional atmospheric and oceanic datasets need to be acquired whenever major hurricanes threaten the United States. The spatial context of these focused three-dimensional datasets over the storm scales is provided by satellites, aircraft, expendable probes released from aircraft, and coastal (both fixed and mobile), moored, and drifting surface platforms. To take full advantage of these new observations, techniques need to be developed to objectively analyze these observations, and initialize models aimed at improving prediction of hurricane track and intensity from global-scale to mesoscale dynamical models. Multinested models allow prediction of all scales from the global, which determine long- term hurricane motion to the convective scale, which affect intensity. Development of an integrated analysis and model forecast system optimizing the use of three-dimensional observations and providing the necessary forecast skill on all relevant spatial scales is required. Detailed diagnostic analyses of these datasets will lead to improved understanding of the physical processes of hurricane motion, intensity change, the atmospheric and oceanic boundary layers, and the air- sea coupling mechanisms. The ultimate aim of this effort is the construction of real-time analyses of storm surge, winds, and rain, prior to and during landfall, to improve warnings and provide local officials with the comprehensive information required for recovery efforts in the hardest hit areas as quickly as possible.

Three-dimensional ocean sensor networks: A survey

NASA Astrophysics Data System (ADS)

Wang, Yu; Liu, Yingjian; Guo, Zhongwen

2012-12-01

The past decade has seen a growing interest in ocean sensor networks because of their wide applications in marine research, oceanography, ocean monitoring, offshore exploration, and defense or homeland security. Ocean sensor networks are generally formed with various ocean sensors, autonomous underwater vehicles, surface stations, and research vessels. To make ocean sensor network applications viable, efficient communication among all devices and components is crucial. Due to the unique characteristics of underwater acoustic channels and the complex deployment environment in three dimensional (3D) ocean spaces, new efficient and reliable communication and networking protocols are needed in design of ocean sensor networks. In this paper, we aim to provide an overview of the most recent advances in network design principles for 3D ocean sensor networks, with focuses on deployment, localization, topology design, and position-based routing in 3D ocean spaces.

El Nino-southern oscillation related fluctuations of the marine carbon cycle

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

Winguth, A.M.E.; Heimann, M.; Kurz, K.D.

The yearly increase in global atmospheric carbon dioxide concentration is not constant, fluctuating around a mean growth rate. Some previous work has been done looking at the relationship of CO2 fluctuations with the El Nino-Southern Oscillation (ENSO) events in the Pacific. This paper describes the response of the three-dimensional ocean circulation model (Hamburg LSG) coupled on-line with a oceanic carbon cycle model (HAMOCC-3) to realistic wind and air temperature field anomalies. The focus is the marine carbon cycle and the interannual variations of carbon fluxes between ocean and atmosphere during the strong El Nino of 1982/83. 53 refs., 14 figs.

Physical-Biological-Optics Model Development and Simulation for the Pacific Ocean and Monterey Bay, California

DTIC Science & Technology

2011-04-01

advanced ROMS-CoSiNE-Optics model in a full three-dimensional environment. We collaborate with Dr. Curt Mobley at Sequoia Scientific to implement...projects. Besides working closely with the modeling group at the NRL and their BioSpace project, we are collaborating with Dr. Curtis Mobley of Sequoia

CITRATE 1.0: Phytoplankton continuous trait-distribution model with one-dimensional physical transport applied to the North Pacific

NASA Astrophysics Data System (ADS)

Chen, Bingzhang; Smith, Sherwood Lan

2018-02-01

Diversity plays critical roles in ecosystem functioning, but it remains challenging to model phytoplankton diversity in order to better understand those roles and reproduce consistently observed diversity patterns in the ocean. In contrast to the typical approach of resolving distinct species or functional groups, we present a ContInuous TRAiT-basEd phytoplankton model (CITRATE) that focuses on macroscopic system properties such as total biomass, mean trait values, and trait variance. This phytoplankton component is embedded within a nitrogen-phytoplankton-zooplankton-detritus-iron model that itself is coupled with a simplified one-dimensional ocean model. Size is used as the master trait for phytoplankton. CITRATE also incorporates trait diffusion for sustaining diversity and simple representations of physiological acclimation, i.e., flexible chlorophyll-to-carbon and nitrogen-to-carbon ratios. We have implemented CITRATE at two contrasting stations in the North Pacific where several years of observational data are available. The model is driven by physical forcing including vertical eddy diffusivity imported from three-dimensional general ocean circulation models (GCMs). One common set of model parameters for the two stations is optimized using the Delayed-Rejection Adaptive Metropolis-Hasting Monte Carlo (DRAM) algorithm. The model faithfully reproduces most of the observed patterns and gives robust predictions on phytoplankton mean size and size diversity. CITRATE is suitable for applications in GCMs and constitutes a prototype upon which more sophisticated continuous trait-based models can be developed.

Development of a three-dimensional variational data assimilation system for the Seto Island Sea, Japan

NASA Astrophysics Data System (ADS)

Kurosawa, K.; Uchiyama, Y.

2016-12-01

By optimally combined ocean models with observation data, numerical oceanic reanalysis and forecast systems allow us to predict the ocean more precisely. In general, data assimilation is exploited to prepare the initial condition for the forecast. This technique has widely been employed in atmospheric prediction, whereas oceanic prediction lags behind weather forecast. Accurate oceanic prediction systems have been demanded for operational purposes such as for fisheries, vessel navigation, marine construction, offshore platform management, marine monitoring, etc. In particular, in crowded harbors and estuaries including the Seto Inland Sea (SIS), Japan, data assimilation has seldom been adapted because data from satellites and Argo floats essential to successful oceanic predictions is desperately limited. In addition, although static data assimilation, typically three-dimensional variational data assimilation (3DVAR), is computationally cheap and statistically optimal, but is not physically balanced. For instance, 3DVAR is known to modify velocity and density fields merely mathematically, yet it does not adequately consider quasi-geostrophic balance, which is generally true in most cases. In the present study, we develop a 3DVAR system for Regional Oceanic Modeling Systems (ROMS) and apply to the high-resolution SIS model in a double nested configuration (Kosako et al., 2015). The SIS is the largest estuary in Japan with a number of autonomous in-situ monitoring of vertical profiles of temperature and salinity, tens of tidal gages, along with continuous surface current measurement using HF radars. We first present a theoretical framework of the 3DVAR algorithm by considering geostrophic and thermal-wind balance to find plausible relationships among physical variables to avoid undesirable modifications. Subsequently, the developed 3DVAR is coupled with the SIS ROMS model to compare the model outcomes against some observation data. The 3DVAR ROMS model for the SIS performs much better than the SIS model without assimilation and demonstrates good model skills with reproducing quite complex flows in the SIS because of its complicated topography with more than 3,000 islands in there. Furthermore we will share technical difficulties encountered during the experiment.

The Southern Ocean's role in ocean circulation and climate transients

NASA Astrophysics Data System (ADS)

Thompson, A. F.; Stewart, A.; Hines, S.; Adkins, J. F.

2017-12-01

The ventilation of deep and intermediate density classes at the surface of the Southern Ocean impacts water mass modification and the air-sea exchange of heat and trace gases, which in turn influences the global overturning circulation and Earth's climate. Zonal variability occurs along the Antarctic Circumpolar Current and the Antarctic margins related to flow-topography interactions, variations in surface boundary conditions, and exchange with northern basins. Information about these zonal variations, and their impact on mass and tracer transport, are suppressed when the overturning is depicted as a two-dimensional (depth-latitude) streamfunction. Here we present an idealized, multi-basin, time-dependent circulation model that applies residual circulation theory in the Southern Ocean and allows for zonal water mass transfer between different ocean basins. This model efficiently determines the temporal evolution of the ocean's stratification, ventilation and overturning strength in response to perturbations in the external forcing. With this model we explore the dynamics that lead to transitions in the circulation structure between multiple, isolated cells and a three-dimensional, "figure-of-eight," circulation in which traditional upper and lower cells are interleaved. The transient model is also used to support a mechanistic explanation of the hemispheric asymmetry and phase lag associated with Dansgaard-Oeschger (DO) events during the last glacial period. In particular, the 200 year lag in southern hemisphere temperatures, following a perturbation in North Atlantic deep water formation, depends critically on the migration of Southern Ocean isopycnal outcropping in response to low-latitude stratification changes. Our results provide a self-consistent dynamical framework to explain various ocean overturning transitions that have occurred over the Earth's last 100,000 years, and motivate an exploration of these mechanisms in more sophisticated climate models.

Three-dimensional Gravity Modeling of Ocean Core Complexes at the Central Indian Ridge

NASA Astrophysics Data System (ADS)

Kim, S. S.; Chandler, M. T.; Pak, S. J.; Son, S. K.

2017-12-01

The spatial distribution of ocean core complexes (OCCs) on mid-ocean ridge flanks can indicate the variation of magmatism and tectonic extension at a given spreading center. A recent study revealed 11 prominent OCCs developed along the middle portion of the Central Indian Ridge (CIR) based on the high-resolution shipboard bathymetry. The CIR is located between the Carlsberg Ridge and the Indian Ocean triple junction. The detailed morphotectonic interpretations from the recent study suggested that the middle ridge segments of the CIR were mainly developed through tectonic extension with little magmatism. Furthermore, the OCCs exposed by detachment faults appear to the main host for active off-axis hydrothermal circulations. Here we form a three-dimensional gravity model to investigate the crustal structures of OCCs developed between 12oS and 14oS at the CIR. These OCCs exhibit domal topographic highs with corrugated surface. The rock samples from these areas include deep-seated rocks such as serpentinized harzburgite and gabbro. A typical gravity study on mid-ocean ridges assumes a constant density contrast along the water-crust interface and constant crustal thickness and removes its gravitational contributions and thermal effects of lithospheric cooling from the free-air gravity anomaly. This approach is effective to distinguish anomalous regions that deviate from the applied crustal and thermal models. The oceanic crust around the OCCs, however, tends to be thinned due to detachment faulting and tectonic extension. In this study, we include multi-layers with different density contrast and variable thickness to approximate gravity anomalies resulting from the OCCs. In addition, we aim to differentiate the geophysical characteristics of the OCCs from the nearby ridge segments and infer tectonic relationship between the OCCs and ridges.

An Examination of a Multi-Scale Three-Dimensional Variational Data Assimilation Scheme in the Kuroshio Extension Using the Naval Coastal Ocean Model

DTIC Science & Technology

2014-01-01

forecasts. Oceanic applications of the MS3DVAR have been Contents lists available at ScienceDirect journal homepage: www.elsevier.com/locate/ csr ...intervals of 12 h. The model grid location and depth can be seen in Fig. 2. In this region the Kuroshio appears as a narrow large magni- tude current in...1997. Introduction to high-frequency radar: reality and myth. Oceanography 10, 36–39. Qu, T., Mitsudera, H., Qui, B., 2001. A climatological view

Global Ocean Prediction with the HYbrid Coordinate Ocean Model, HYCOM

NASA Astrophysics Data System (ADS)

Chassignet, E.

A broad partnership of institutions is collaborating in developing and demonstrating the performance and application of eddy-resolving, real-time global and Atlantic ocean prediction systems using the the HYbrid Coordinate Ocean Model (HYCOM). These systems will be transitioned for operational use by both the U.S. Navy at the Naval Oceanographic Office (NAVOCEANO), Stennis Space Center, MS, and the Fleet Numerical Meteorology and Oceanography Centre (FNMOC), Monterey, CA, and by NOAA at the National Centers for Environmental Prediction (NCEP), Washington, D.C. These systems will run efficiently on a variety of massively parallel computers and will include sophisticated data assimilation techniques for assimilation of satellite altimeter sea surface height and sea surface temperature as well as in situ temperature, salinity, and float displacement. The Partnership addresses the Global Ocean Data Assimilation Experiment (GODAE) goals of three-dimensional (3D) depiction of the ocean state at fine resolution in real-time and provision of boundary conditions for coastal and regional models. An overview of the effort will be presented.

Broadband and polarization reflectors in the lookdown, Selene vomer.

PubMed

Zhao, Shulei; Brady, Parrish Clawson; Gao, Meng; Etheredge, Robert Ian; Kattawar, George W; Cummings, Molly E

2015-03-06

Predator evasion in the open ocean is difficult because there are no objects to hide behind. The silvery surface of fish plays an important role in open water camouflage. Various models have been proposed to account for the broadband reflectance by the fish skin that involve one-dimensional variations in the arrangement of guanine crystal reflectors, yet the three-dimensional organization of these guanine platelets have not been well characterized. Here, we report the three-dimensional organization and the optical properties of integumentary guanine platelets in a silvery marine fish, the lookdown (Selene vomer). Our structural analysis and computational modelling show that stacks of guanine platelets with random yaw angles in the fish skin produce broadband reflectance via colour mixing. Optical axes of the guanine platelets and the collagen layer are aligned closely and provide bulk birefringence properties that influence the polarization reflectance by the skin. These data demonstrate how the lookdown preserves or alters polarization states at different incident polarization angles. These optical properties resulted from the organization of these guanine platelets and the collagen layer may have implications for open ocean camouflage in varying light fields. © 2015 The Author(s) Published by the Royal Society. All rights reserved.

Ecosystem behavior at Bermuda Station [open quotes]S[close quotes] and ocean weather station [open quotes]India[close quotes]: A general circulation model and observational analysis

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

Fasham, M.J.R.; Sarmiento, J.L.; Slater, R.D.

1993-06-01

One important theme of modern biological oceanography has been the attempt to develop models of how the marine ecosystem responds to variations in the physical forcing functions such as solar radiation and the wind field. The authors have addressed the problem by embedding simple ecosystem models into a seasonally forced three-dimensional general circulation model of the North Atlantic ocean. In this paper first, some of the underlying biological assumptions of the ecosystem model are presented, followed by an analysis of how well the model predicts the seasonal cycle of the biological variables at Bermuda Station s' and Ocean Weather Stationmore » India. The model gives a good overall fit to the observations but does not faithfully model the whole seasonal ecosystem model. 57 refs., 25 figs., 5 tabs.« less

Lateral variations in upper-mantle seismic anisotropy in the Pacific from inversion of a surface-wave dispersion dataset

NASA Astrophysics Data System (ADS)

Eddy, C. L.; Ekstrom, G.; Nettles, M.; Gaherty, J. B.

2017-12-01

We present a three-dimensional model of the anisotropic velocity structure of the Pacific lithosphere and asthenosphere. The presence of seismic anisotropy in the oceanic upper mantle provides information about the geometry of flow in the mantle, the nature of the lithosphere-asthenosphere boundary, and the possible presence of partial melt in the asthenosphere. Our dataset consists of fundamental-mode dispersion for Rayleigh and Love waves measured between 25-250 s with paths crossing the Pacific Ocean. We invert the phase anomaly measurements directly for three-dimensional anisotropic velocity structure. Our models are radially anisotropic and include the full set of elastic parameters that describe azimuthal variations in velocity (e.g. Gc, Gs). We investigate the age dependence of seismic velocity and radial anisotropy and find that there are significant deviations from the velocities predicted by a simple oceanic plate cooling model. We observe strong radial anisotropy with vsh > vsv in the asthenosphere of the central Pacific. We investigate the radial anisotropy in the shallow lithosphere, where previous models have reported conflicting results. There is a contrast in both upper-mantle isotropic velocities and radial anisotropy between the Pacific and Nazca plates, across the East Pacific Rise. We also investigate lateral variations in azimuthal anisotropy throughout the Pacific upper mantle and find that there are large areas over which the anisotropy fast axis does not align with absolute plate motion, suggesting the presence of small-scale convection or pressure-driven flow beneath the base of the oceanic plate.

Studies on Three-Dimensional Slamming on Slender Ships

DTIC Science & Technology

2007-05-10

Objectives This study is a joint effort among MIT, Seoul National University (SNU), and Korean Research Institute of Ships and Ocean Engineering ( KRISO ...model and laboratory experiments were carried out by SNU/ KRISO . The key responsibility of MIT is to investigate the feasibility and effectiveness of

Evidence and mechanism of Hurricane Fran-Induced ocean cooling in the Charleston Trough

NASA Astrophysics Data System (ADS)

Xie, Lian; Pietrafesa, L. J.; Bohm, E.; Zhang, C.; Li, X.

Evidence of enhanced sea surface cooling during and following the passage of Hurricane Fran in September 1996 over an oceanic depression located on the ocean margin offshore of Charleston, South Carolina (referred to as the Charleston Trough), [Pietrafesa, 1983] is documented. Approximately 4C° of sea surface temperature (SST) reduction within the Charleston Trough following the passage of Hurricane Fran was estimated based on SST imagery from Advanced Very High Resolution Radiometer (AVHRR) on the NOAA-14 polar orbiting satellite. Simulations using a three-dimensional coastal ocean model indicate that the largest SST reduction occurred within the Charleston Trough. This SST reduction can be explained by oceanic mixing due to storm-induced internal inertia-gravity waves.

The Pale Blue Dot: Utilizing Real World Globes in High School and Undergraduate Oceanography Classrooms

NASA Astrophysics Data System (ADS)

Rogers, D. B.

2017-12-01

Geoscience classrooms have benefitted greatly from the use of interactive, dry-erasable globes to supplement instruction on topics that require three-dimensional visualization, such as seismic wave propagation and the large-scale movements of tectonic plates. Indeed, research by Bamford (2013) demonstrates that using three-dimensional visualization to illustrate complex processes enhances student comprehension. While some geoscience courses tend to bake-in lessons on visualization, other disciplines of earth science that require three-dimensional visualization, such as oceanography, tend to rely on students' prior spatial abilities. In addition to spatial intelligence, education on the three-dimensional structure of the ocean requires knowledge of the external processes govern the behavior of the ocean, as well as the vertical and lateral distribution of water properties around the globe. Presented here are two oceanographic activities that utilize RealWorldGlobes' dry-erase globes to supplement traditional oceanography lessons on thermohaline and surface ocean circulation. While simultaneously promoting basic plotting techniques, mathematical calculations, and unit conversions, these activities touch on the processes that govern global ocean circulation, the principles of radiocarbon dating, and the various patterns exhibited by surface ocean currents. These activities challenge students to recognize inherent patterns within their data and synthesize explanations for their occurrence. Spatial visualization and critical thinking are integral to any geoscience education, and the combination of these abilities with engaging hands-on activities has the potential to greatly enhance oceanography education in both secondary and postsecondary settings

How well-connected is the surface of the global ocean?

PubMed

Froyland, Gary; Stuart, Robyn M; van Sebille, Erik

2014-09-01

The Ekman dynamics of the ocean surface circulation is known to contain attracting regions such as the great oceanic gyres and the associated garbage patches. Less well-known are the extents of the basins of attractions of these regions and how strongly attracting they are. Understanding the shape and extent of the basins of attraction sheds light on the question of the strength of connectivity of different regions of the ocean, which helps in understanding the flow of buoyant material like plastic litter. Using short flow time trajectory data from a global ocean model, we create a Markov chain model of the surface ocean dynamics. The surface ocean is not a conservative dynamical system as water in the ocean follows three-dimensional pathways, with upwelling and downwelling in certain regions. Using our Markov chain model, we easily compute net surface upwelling and downwelling, and verify that it matches observed patterns of upwelling and downwelling in the real ocean. We analyze the Markov chain to determine multiple attracting regions. Finally, using an eigenvector approach, we (i) identify the five major ocean garbage patches, (ii) partition the ocean into basins of attraction for each of the garbage patches, and (iii) partition the ocean into regions that demonstrate transient dynamics modulo the attracting garbage patches.

Detection of Rossby Waves in Multi-Parameters in Multi-Mission Satellite Observations and HYCOM Simulations in the Indian Ocean

NASA Technical Reports Server (NTRS)

Subrahmanyam, Bulusu; Heffner, David M.; Cromwell, David; Shriver, Jay F.

2009-01-01

Rossby waves are difficult to detect with in situ methods. However, as we show in this paper, they can be clearly identified in multi-parameters in multi-mission satellite observations of sea surface height (SSH), sea surface temperature (SST) and ocean color observations of chlorophyll-a (chl-a), as well as 1/12-deg global HYbrid Coordinate Ocean Model (HYCOM) simulations of SSH, SST and sea surface salinity (SSS) in the Indian Ocean. While the surface structure of Rossby waves can be elucidated from comparisons of the signal in different sea surface parameters, models are needed to gain direct information about how these waves affect the ocean at depth. The first three baroclinic modes of the Rossby waves are inferred from the Fast Fourier Transform (FFT), and two-dimensional Radon Transform (2D RT). At many latitudes the first and second baroclinic mode Rossby wave phase speeds from satellite observations and model parameters are identified.

Evaluation of Regional Travel-Time and Location Improvement Along the Tethyan Margin Using a New Three-Dimensional Velocity Model

DTIC Science & Technology

2008-09-01

part of the Atlantic Ocean for reference. The Moho depth result is broadly consistent with CRUST2.0, except in mid-northern Africa, where the crust...plate boundaries is shown by the pink line in Figure 1. The interaction of these five major tectonic plates with each other and with several microplates ...acquired from literatures. Artificial point constraints of 10 km depth are placed to the Atlantic and Indian Oceans where measurements are absent

Semiannual progress report, April - September 1991

NASA Technical Reports Server (NTRS)

1991-01-01

Research conducted during the past year in the climate and modeling programs has concentrated on the development of appropriate atmospheric and upper ocean models, and preliminary applications of these models. Principal models are a one-dimensional radiative-convective model, a three dimensional global climate model, and an upper ocean model. Principal applications have been the study of the impact of CO2, aerosols, and the solar constant on climate. Progress was made in the 3-D model development towards physically realistic treatment of these processes. In particular, a map of soil classifications on 1 degree by 1 degree resolution has now been digitized, and soil properties have been assigned to each soil type. Using this information about soil properties, a method has been developed to simulate the hydraulic behavior of the soils of the world. This improved treatment of soil hydrology, together with the seasonally varying vegetation cover, will provide a more realistic study of the role of the terrestrial biota in climate change. A new version of the climate model was created which follows the isotopes of water and sources of water throughout the planet.

Transient rheology of the oceanic asthenosphere following the 2012 Indian Ocean Earthquake inferred from geodetic data

NASA Astrophysics Data System (ADS)

Pratama, Cecep; Ito, Takeo; Sasajima, Ryohei; Tabei, Takao; Kimata, Fumiaki; Gunawan, Endra; Ohta, Yusaku; Yamashina, Tadashi; Ismail, Nazli; Nurdin, Irwandi; Sugiyanto, Didik; Muksin, Umar; Meilano, Irwan

2017-10-01

Postseismic motion in the middle-field (100-500 km from the epicenter) geodetic data resulting from the 2012 Indian Ocean earthquake exhibited rapid change during the two months following the rupture. This pattern probably indicates multiple postseismic deformation mechanisms and might have been controlled by transient rheology. Therefore, the relative contribution of transient rheology in the oceanic asthenosphere and afterslip in the oceanic lithosphere should be incorporated to explain short- and long-term transitional features of postseismic signals. In this study, using two years of post-earthquake geodetic data from northern Sumatra, a three-dimensional spherical-earth finite-element model was constructed based on a heterogeneous structure and incorporating transient rheology. A rheology model combined with stress-driven afterslip was estimated. Our best-fit model suggests an oceanic lithosphere thickness of 75 km with oceanic asthenosphere viscosity values of 1 × 1017 Pa s and 2 × 1018 Pa s for the Kelvin and Maxwell viscosity models, respectively. The model results indicate that horizontal landward motion and vertical uplift in northern Sumatra require viscoelastic relaxation of the oceanic asthenosphere coupled with afterslip in the lithosphere. The present study demonstrates that transient rheology is essential for reproducing the rapidly changing motion of postseismic deformation in the middle-field area.

A contrastive study on the influences of radial and three-dimensional satellite gravity gradiometry on the accuracy of the Earth's gravitational field recovery

NASA Astrophysics Data System (ADS)

Zheng, Wei; Hsu, Hou-Tse; Zhong, Min; Yun, Mei-Juan

2012-10-01

The accuracy of the Earth's gravitational field measured from the gravity field and steady-state ocean circulation explorer (GOCE), up to 250 degrees, influenced by the radial gravity gradient Vzz and three-dimensional gravity gradient Vij from the satellite gravity gradiometry (SGG) are contrastively demonstrated based on the analytical error model and numerical simulation, respectively. Firstly, the new analytical error model of the cumulative geoid height, influenced by the radial gravity gradient Vzz and three-dimensional gravity gradient Vij are established, respectively. In 250 degrees, the GOCE cumulative geoid height error measured by the radial gravity gradient Vzz is about 2½ times higher than that measured by the three-dimensional gravity gradient Vij. Secondly, the Earth's gravitational field from GOCE completely up to 250 degrees is recovered using the radial gravity gradient Vzz and three-dimensional gravity gradient Vij by numerical simulation, respectively. The study results show that when the measurement error of the gravity gradient is 3 × 10-12/s2, the cumulative geoid height errors using the radial gravity gradient Vzz and three-dimensional gravity gradient Vij are 12.319 cm and 9.295 cm at 250 degrees, respectively. The accuracy of the cumulative geoid height using the three-dimensional gravity gradient Vij is improved by 30%-40% on average compared with that using the radial gravity gradient Vzz in 250 degrees. Finally, by mutual verification of the analytical error model and numerical simulation, the orders of magnitude from the accuracies of the Earth's gravitational field recovery make no substantial differences based on the radial and three-dimensional gravity gradients, respectively. Therefore, it is feasible to develop in advance a radial cold-atom interferometric gradiometer with a measurement accuracy of 10-13/s2-10-15/s2 for precisely producing the next-generation GOCE Follow-On Earth gravity field model with a high spatial resolution.

Visualizing characteristics of ocean data collected during the Shuttle Imaging Radar-B experiment

NASA Technical Reports Server (NTRS)

Tilley, David G.

1991-01-01

Topographic measurements of sea surface elevation collected by the Surface Contour Radar (SCR) during NASA's Shuttle Imaging Radar (SIR-B) experiment are plotted as three dimensional surface plots to observe wave height variance along the track of a P-3 aircraft. Ocean wave spectra were computed from rotating altimeter measurements acquired by the Radar Ocean Wave Spectrometer (ROWS). Fourier power spectra computed from SIR-B synthetic aperture radar (SAR) images of the ocean are compared to ROWS surface wave spectra. Fourier inversion of SAR spectra, after subtraction of spectral noise and modeling of wave height modulation, yields topography similar to direct measurements made by SCR. Visual perspectives on the SCR and SAR ocean data are compared. Threshold distinctions between surface elevation and texture modulations of SAR data are considered within the context of a dynamic statistical model of rough surface scattering. The result of these endeavors is insight as to the physical mechanism governing the imaging of ocean waves with SAR.

West Florida shelf circulation and temperature budget for the 1998 fall transition

NASA Astrophysics Data System (ADS)

He, Ruoying; Weisberg, Robert H.

2003-05-01

Mid-latitude continental shelves undergo a fall transition as the net heat flux changes from warming to cooling. Using in situ data and a numerical model we investigate the circulation on the west Florida shelf (WFS) for the fall transition of 1998. The model is a regional adaptation of the primitive equation, Princeton Ocean Model forced by NCEP reanalysis wind, air pressure, and heat flux fields, plus river inflows. After comparison with observations the model is used to draw inferences on the seasonal and synoptic scale features of the shelf circulation. By running twin experiments, one without and the other with an idealized Loop Current (LC), we explore the relative importance of local versus deep-ocean forcing. We find that local forcing largely controls the inner-shelf circulation, including changes from the Florida Panhandle in the north to regions farther south. The effects of the LC in fall 1998 are to reinforce the mid-shelf currents and to increase the across-shelf transports in the bottom Ekman layer, thereby accentuating the shoreward transport of cold, nutrient rich water of deep-ocean origin. A three-dimensional analysis of the temperature budget reveals that surface heat flux largely controls both the seasonal and synoptic scale temperature variations. Surface cooling leads to convective mixing that rapidly alters temperature gradients. One interesting consequence is that upwelling can result in near-shore warming as warmer offshore waters are advected landward. The temperature balances on the shelf are complex and fully three-dimensional.

MEASUREMENTS OF PAST 14C LEVELS AND 13C/12C RATIOS IN THE SURFACE WATERS OF THE WORLD'S SUBPOLAR OCEANS.

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

Brown, T A

2010-04-22

Under this project we have developed methods that allow the reconstruction of past {sup 14}C levels of the surface waters of the subpolar North Pacific Ocean by measuring the {sup 14}C contents of archived salmon scales. The overall goal of this research was to reduce of the uncertainty in the uptake of fossil CO{sub 2} by the oceans and thereby improve the quantification of the global carbon cycle and to elucidate the fate of anthropogenic CO{sub 2}. Ocean General Circulation Models (OGCMs), with their three dimensional global spatial coverage and temporal modeling capabilities, provide the best route to accurately calculatingmore » the total uptake of CO{sub 2} by the oceans and, hence, to achieving the desired reduction in uncertainty. {sup 14}C has played, and continues to play, a central role in the validation of the OGCMs calculations, particularly with respect to those model components which govern the uptake of CO{sub 2} from the atmosphere and the transport of this carbon within the oceans. Under this project, we have developed time-series records of the {sup 14}C levels of the surface waters of three areas of the subpolar North Pacific Ocean. As the previously available data on the time-history of oceanic surface water {sup 14}C levels are very limited, these time-series records provide significant new {sup 14}C data to constrain and validate the OGCMs.« less

Pliocene three-dimensional global ocean temperature reconstruction

USGS Publications Warehouse

Dowsett, H.J.; Robinson, M.M.; Foley, K.M.

2009-01-01

The thermal structure of the mid-Piacenzian ocean is obtained by combining the Pliocene Research, Interpretation and Synoptic Mapping Project (PRISM3) multiproxy sea-surface temperature (SST) reconstruction with bottom water temperature estimates from 27 locations produced using Mg/Ca paleothermometry based upon the ostracod genus Krithe. Deep water temperature estimates are skewed toward the Atlantic Basin (63% of the locations) and represent depths from 1000m to 4500 m. This reconstruction, meant to serve as a validation data set as well as an initialization for coupled numerical climate models, assumes a Pliocene water mass framework similar to that which exists today, with several important modifications. The area of formation of present day North Atlantic Deep Water (NADW) was expanded and extended further north toward the Arctic Ocean during the mid-Piacenzian relative to today. This, combined with a deeper Greenland-Scotland Ridge, allowed a greater volume of warmer NADW to enter the Atlantic Ocean. In the Southern Ocean, the Polar Front Zone was expanded relative to present day, but shifted closer to the Antarctic continent. This, combined with at least seasonal reduction in sea ice extent, resulted in decreased Antarctic Bottom Water (AABW) production (relative to present day) as well as possible changes in the depth of intermediate waters. The reconstructed mid-Piacenzian three-dimensional ocean was warmer overall than today, and the hypothesized aerial extent of water masses appears to fit the limited stable isotopic data available for this time period. ?? Author(s) 2009.

Spiraling pathways of global deep waters to the surface of the Southern Ocean.

PubMed

Tamsitt, Veronica; Drake, Henri F; Morrison, Adele K; Talley, Lynne D; Dufour, Carolina O; Gray, Alison R; Griffies, Stephen M; Mazloff, Matthew R; Sarmiento, Jorge L; Wang, Jinbo; Weijer, Wilbert

2017-08-02

Upwelling of global deep waters to the sea surface in the Southern Ocean closes the global overturning circulation and is fundamentally important for oceanic uptake of carbon and heat, nutrient resupply for sustaining oceanic biological production, and the melt rate of ice shelves. However, the exact pathways and role of topography in Southern Ocean upwelling remain largely unknown. Here we show detailed upwelling pathways in three dimensions, using hydrographic observations and particle tracking in high-resolution models. The analysis reveals that the northern-sourced deep waters enter the Antarctic Circumpolar Current via southward flow along the boundaries of the three ocean basins, before spiraling southeastward and upward through the Antarctic Circumpolar Current. Upwelling is greatly enhanced at five major topographic features, associated with vigorous mesoscale eddy activity. Deep water reaches the upper ocean predominantly south of the Antarctic Circumpolar Current, with a spatially nonuniform distribution. The timescale for half of the deep water to upwell from 30° S to the mixed layer is ~60-90 years.Deep waters of the Atlantic, Pacific and Indian Oceans upwell in the Southern Oceanbut the exact pathways are not fully characterized. Here the authors present a three dimensional view showing a spiralling southward path, with enhanced upwelling by eddy-transport at topographic hotspots.

Geosat Data Assimilation with Application to the Eastern North Atlantic

NASA Technical Reports Server (NTRS)

Stammer, Detlef

1997-01-01

An attempt is made to determine the three-dimensional ocean circulation from satellite altimeter measurements by assimilating Geosat sea surface height data into an eddy-resolving QuasiGeostrophic (QG) model of the eastern North Atlantic Ocean. Results are tested against independent information from hydrographic field observations and moored current meter data collected during the Geosat ERM. The comparison supports the concept of inferring aspects of the three-dimensional flow field from sea surface height observations by combining altimetric measurements with the dynamics of ocean circulation models. A Holland-type QG model with open boundaries is set up on a 2000 km X 2000 km domain of the eastern North Atlantic between 25 deg. and 45 deg. N, 32 deg. and 8 deg. W. By using a simple nudging technique, about two years of Geosat altimeter data are assimilated into the model every five days as space-time objective analyses on the model grid. The error information resulting from the analysis is used during the assimilation procedure to account for data uncertainties. Results show an intense eddy field, which in the surface layer interacts with a meandering Azores Front. Compared to Geosat, the model leads to smoothed fields that follow the observations. Model simulations are significantly correlated with hydrographic data from March 1988 and June 1989, both close to the surface and in the subsurface. Good agreement is also found between the model velocity fields and moored current meter data in the top two model layers. The agreement is visually weak in the bottom layer, although a coherence analysis reveals an agreement between the model simulation and current meter data over the full water column at periods exceeding 80 days.

Comparison of global sst analyses for atmospheric data assimilation

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

Phoebus, P.A.; Cummings, J.A.

1995-03-17

Traditionally, atmospheric models were executed using a climatological estimate of the sea surface temperature (SST) to define the marine boundary layer. More recently, particularly since the deployment of remote sensing instruments and the advent of multichannel SST observations atmospheric models have been improved by using more timely estimates of the actual state of the ocean. Typically, some type of objective analysis is performed using the data from satellites along with ship, buoy, and bathythermograph observations, and perhaps even climatology, to produce a weekly or daily analysis of global SST. Some of the earlier efforts to produce real-time global temperature analysesmore » have been described by Clancy and Pollak (1983) and Reynolds (1988). However, just as new techniques have been developed for atmospheric data assimilation, improvements have been made to ocean data assimilation systems as well. In 1988, the U.S. Navy`s Fleet Numerical Meteorology and Oceanography Center (FNMOC) implemented a global three-dimensional ocean temperature analysis that was based on the optimum interpolation methodology (Clancy et al., 1990). This system, the Optimum Thermal Interpolation System (OTIS 1.0), was initially distributed on a 2.50 resolution grid, and was later modified to generate fields on a 1.250 grid (OTIS 1.1; Clancy et al., 1992). Other optimum interpolation-based analyses (OTIS 3.0) were developed by FNMOC to perform high-resolution three-dimensional ocean thermal analyses in areas with strong frontal gradients and clearly defined water mass characteristics.« less

An Integrated Modeling and Observational Study of Three-Dimensional Upper Ocean Boundary Layer Dynamics and Parameterizations

DTIC Science & Technology

2012-04-26

for public release ; distribution is unlimited. 4 Figure 3. First Kuroshio survey. Colors show MODIS SST image. Inset shows the survey...Institute of Diabetes and Digestive and Kidney Diseases, Bethesda, MD 20892, USA. *To whom correspondence should be addressed. E-mail: stevens@scripps.edu 15

Real-time Assimilation of Altimeter Derived Synthetic Profiles Into a Global version of the Naval Research Laboratory's Coastal Ocean Model (NCOM)

NASA Astrophysics Data System (ADS)

Rhodes, R. C.; Barron, C. N.; Fox, D. N.; Smedstad, L. F.

2001-12-01

A global implementation of the Navy Coastal Ocean Model (NCOM), developed by the Naval Research Laboratory (NRL) at Stennis Space Center is currently running in real-time and is planned for transition to the Naval Oceanographic Office (NAVOCEANO) in 2002. The model encompasses the open ocean to 5 m depth on a curvilinear global model grid with 1/8 degree grid spacing at 45N, extending from 80 S to a complete arctic cap with grid singularities mapped into Canada and Russia. Vertically, the model employs 41 sigma-z levels with sigma in the upper-ocean and coastal regions and z in the deeper ocean. The Navy Operational Global Atmospheric Prediction System (NOGAPS) provides 6-hourly wind stresses and heat fluxes for forcing, while the operational Modular Ocean Data Assimilation System (MODAS) provides the background climatology and tools for data pre-processing. Operationally available sea surface temperature (SST) and altimetry (SSH) data are assimilated into the NAVOCEANO global 1/8 degree MODAS 2-D analysis and the 1/16 degree Navy Layered Ocean Model (NLOM) to provide analyses and forecasts of SSH and SST. The 2-D SSH and SST nowcast fields are used as input to the MODAS synthetic climatology database to yield three-dimensional fields of synthetic temperature and salinity for assimilation into global NCOM. The synthetic profiles are weighted higher at depth in the assimilation process to allow the numerical model to properly develop the mixed-layer structure driven by the real-time atmospheric forcing. Global NCOM nowcasts and forecasts provide a valuable resource for rapid response to the varied and often unpredictable operational requests for 3-dimensional fields of ocean temperature, salinity, and currents. In some cases, the resolution of the global product is sufficient for guidance. In cases requiring higher resolution, the global product offers a quick overview of local circulation and provides initial and boundary conditions for higher resolution coastal models that may be more specialized for a particular task or domain. Nowcast and forecast results are presented globally and in selected areas of interest and model results are compared with historical and concurrent observations and analyses.

Assessing the importance of internal tide scattering in the deep ocean

NASA Astrophysics Data System (ADS)

Haji, Maha; Peacock, Thomas; Carter, Glenn; Johnston, T. M. Shaun

2014-11-01

Tides are one of the main sources of energy input to the deep ocean, and the pathways of energy transfer from barotropic tides to turbulent mixing scales via internal tides are not well understood. Large-scale (low-mode) internal tides account for the bulk of energy extracted from barotropic tides and have been observed to propagate over 1000 km from their generation sites. We seek to examine the fate of these large-scale internal tides and the processes by which their energy is transferred, or ``scattered,'' to small-scale (high-mode) internal tides, which dissipate locally and are responsible for internal tide driven mixing. The EXperiment on Internal Tide Scattering (EXITS) field study conducted in 2010-2011 sought to examine the role of topographic scattering at the Line Islands Ridge. The scattering process was examined via data from three moorings equipped with moored profilers, spanning total depths of 3000--5000 m. The results of our field data analysis are rationalized via comparison to data from two- and three-dimensional numerical models and a two-dimensional analytical model based on Green function theory.

Dynamical analysis of a satellite-observed anticyclonic eddy in the northern Bering Sea

NASA Astrophysics Data System (ADS)

Li, Yineng; Li, Xiaofeng; Wang, Jia; Peng, Shiqiu

2016-05-01

The characteristics and evolution of a satellite-observed anticyclonic eddy in the northern Bering Sea during March and April 1999 are investigated using a three-dimensional Princeton Ocean Model (POM). The anticyclonic-like current pattern and asymmetric feature of the eddy were clearly seen in the synthetic aperture radar (SAR), sea surface temperature, and ocean color images in April 1999. The results from model simulation reveal the three-dimensional structure of the anticyclonic eddy, its movement, and dissipation. Energy analysis indicates that the barotropic instability (BTI) is the main energy source for the growth of the anticyclonic eddy. The momentum analysis further reveals that the larger magnitude of the barotropic pressure gradient in the meridional direction causes the asymmetry of the anticyclonic eddy in the zonal and meridional directions, while the different magnitudes of the meridional baroclinic pressure gradient are responsible for the different intensity of currents between the northern and southern parts of the anticyclonic eddy. This article was corrected on 23 JUL 2016. See the end of the full text for details.

Ocean-driven heating of Europa's icy shell at low latitudes

NASA Astrophysics Data System (ADS)

Soderlund, K. M.; Schmidt, B. E.; Wicht, J.; Blankenship, D. D.

2014-01-01

The ice shell of Jupiter's moon Europa is marked by regions of disrupted ice known as chaos terrains that cover up to 40% of the satellite's surface, most commonly occurring within 40° of the equator. Concurrence with salt deposits implies a coupling between the geologically active ice shell and the underlying liquid water ocean at lower latitudes. Europa's ocean dynamics have been assumed to adopt a two-dimensional pattern, which channels the moon's internal heat to higher latitudes. Here we present a numerical model of thermal convection in a thin, rotating spherical shell where small-scale convection instead adopts a three-dimensional structure and is more vigorous at lower latitudes. Global-scale currents are organized into three zonal jets and two equatorial Hadley-like circulation cells. We find that these convective motions transmit Europa's internal heat towards the surface most effectively in equatorial regions, where they can directly influence the thermo-compositional state and structure of the ice shell. We suggest that such heterogeneous heating promotes the formation of chaos features through increased melting of the ice shell and subsequent deposition of marine ice at low latitudes. We conclude that Europa's ocean dynamics can modulate the exchange of heat and materials between the surface and interior and explain the observed distribution of chaos terrains.

Challenges in Ocean Data Assimilation for the US West Coast

NASA Astrophysics Data System (ADS)

Li, Z.; Chao, Y.; Farrara, J.; Wang, X.

2006-12-01

A three-dimensional variational data assimilation (3DVAR) system has been developed for the Regional Ocean Modeling System (ROMS), and it is called ROMS-DAS. This system provides a capability of predicting meso- to small-scale variations with temporal scales from hours to days in the coastal oceans. To cope with the particular difficulties that result from complex coastlines and bottom topography, unbalanced flows and sparse observations, ROMS-DAS utilizes several novel strategies. These strategies include the implementation of three-dimensional anisotropic and inhomogeneous error correlations, application of particular weak dynamic constraints, and implementation of efficient and reliable algorithms for minimizing the cost function. The ROMS-DAS system was applied in field experiments for Monterey Bay during both 2003 (Autonomous Ocean Sampling Network - AOSN) and 2006 (MB06). These two experiments included intensive data collection from a variety of observational platforms, including satellites, airplanes, High Frequency radars, Acoustic Doppler Current Profilers, ships, drifters, buoys, autonomous underwater vehicles (AUV), and particularly a fleet of undersea gliders. Using these data sets, various data assimilation experiments were performed to address several major data assimilation challenges that arise from multi-scales structures, inhomogeneous properties, dynamical imbalance of the flow, and tides. Basing on these experiments, a set of strategies were formulated to deal with those challenges.

Ocean forecasting in terrain-following coordinates: Formulation and skill assessment of the Regional Ocean Modeling System

USGS Publications Warehouse

Haidvogel, D.B.; Arango, H.; Budgell, W.P.; Cornuelle, B.D.; Curchitser, E.; Di, Lorenzo E.; Fennel, K.; Geyer, W.R.; Hermann, A.J.; Lanerolle, L.; Levin, J.; McWilliams, J.C.; Miller, A.J.; Moore, A.M.; Powell, T.M.; Shchepetkin, A.F.; Sherwood, C.R.; Signell, R.P.; Warner, J.C.; Wilkin, J.

2008-01-01

Systematic improvements in algorithmic design of regional ocean circulation models have led to significant enhancement in simulation ability across a wide range of space/time scales and marine system types. As an example, we briefly review the Regional Ocean Modeling System, a member of a general class of three-dimensional, free-surface, terrain-following numerical models. Noteworthy characteristics of the ROMS computational kernel include: consistent temporal averaging of the barotropic mode to guarantee both exact conservation and constancy preservation properties for tracers; redefined barotropic pressure-gradient terms to account for local variations in the density field; vertical interpolation performed using conservative parabolic splines; and higher-order, quasi-monotone advection algorithms. Examples of quantitative skill assessment are shown for a tidally driven estuary, an ice-covered high-latitude sea, a wind- and buoyancy-forced continental shelf, and a mid-latitude ocean basin. The combination of moderate-order spatial approximations, enhanced conservation properties, and quasi-monotone advection produces both more robust and accurate, and less diffusive, solutions than those produced in earlier terrain-following ocean models. Together with advanced methods of data assimilation and novel observing system technologies, these capabilities constitute the necessary ingredients for multi-purpose regional ocean prediction systems. 

Preparing to Predict: The Second Autonomous Ocean Sampling Network (AOSN-II) Experiment in the Monterey Bay

DTIC Science & Technology

2008-06-06

MBARI buoy M2, profiling to 250 m depth (Figure 1). Instruments on board included a CTD, fluorometer, oxygen and nitrate sensors, bioluminescence, and...dimensional multiscale ocean variability: Massachusetts Bay. Journal of Marine Systems, Special issue on “Three-dimensional ocean circulation: Lagrangian...Oceanography”, T. Paluszkiewicz and S. Harper, Eds., Vol. 19, 1, 172-183. Liang, X.S. and Anderson, D.G.M. (2007) Multiscale Window Transform, SIAM J

Global environmental effects of impact-generated aerosols: Results from a general circulation model, revision 1

NASA Technical Reports Server (NTRS)

Covey, Curt; Ghan, Steven J.; Walton, John J.; Weissman, Paul R.

1989-01-01

Interception of sunlight by the high altitude worldwide dust cloud generated by impact of a large asteroid or comet would lead to substantial land surface cooling, according to our three-dimensional atmospheric general circulation model (GCM). This result is qualitatively similar to conclusions drawn from an earlier study that employed a one-dimensional atmospheric model, but in the GCM simulation the heat capacity of the oceans substantially mitigates land surface cooling, an effect that one-dimensional models cannot quantify. On the other hand, the low heat capacity of the GCM's land surface allows temperatures to drop more rapidly in the initial stage of cooling than in the one-dimensional model study. These two differences between three-dimensional and one-dimensional model simulations were noted previously in studies of nuclear winter; GCM-simulated climatic changes in the Alvarez-inspired scenario of asteroid/comet winter, however, are more severe than in nuclear winter because the assumed aerosol amount is large enough to intercept all sunlight falling on earth. Impacts of smaller objects could also lead to dramatic, though less severe, climatic changes, according to our GCM. Our conclusion is that it is difficult to imagine an asteroid or comet impact leading to anything approaching complete global freezing, but quite reasonable to assume that impacts at the Alvarez level, or even smaller, dramatically alter the climate in at least a patchy sense.

Natural biogeochemical cycle of mercury in a global three-dimensional ocean tracer model

NASA Astrophysics Data System (ADS)

Zhang, Yanxu; Jaeglé, Lyatt; Thompson, LuAnne

2014-05-01

We implement mercury (Hg) biogeochemistry in the offline global 3-D ocean tracer model (OFFTRAC) to investigate the natural Hg cycle, prior to any anthropogenic input. The simulation includes three Hg tracers: dissolved elemental (Hg0aq), dissolved divalent (HgIIaq), and particle-bound mercury (HgPaq). Our Hg parameterization takes into account redox chemistry in ocean waters, air-sea exchange of Hg0, scavenging of HgIIaq onto sinking particles, and resupply of HgIIaq at depth by remineralization of sinking particles. Atmospheric boundary conditions are provided by a global simulation of the natural atmospheric Hg cycle in the GEOS-Chem model. In the surface ocean, the OFFTRAC model predicts global mean concentrations of 0.16 pM for total Hg, partitioned as 80% HgIIaq, 14% Hg0aq, and 6% HgPaq. Total Hg concentrations increase to 0.38 pM in the thermocline/intermediate waters (between the mixed layer and 1000 m depth) and 0.82 pM in deep waters (below 1000 m), reflecting removal of Hg from the surface to the subsurface ocean by particle sinking followed by remineralization at depth. Our model predicts that Hg concentrations in the deep North Pacific Ocean (>2000 m) are a factor of 2-3 higher than in the deep North Atlantic Ocean. This is the result of cumulative input of Hg from particle remineralization as deep waters transit from the North Atlantic to the North Pacific on their ~2000 year journey. The model is able to reproduce the relatively uniform concentrations of total Hg observed in the old deep waters of the North Pacific Ocean (observations: 1.2 ± 0.4 pM; model: 1.1 ± 0.04 pM) and Southern Ocean (observations: 1.1 ± 0.2 pM; model: 0.8 ± 0.02 pM). However, the modeled concentrations are factors of 5-6 too low compared to observed concentrations in the surface ocean and in the young water masses of the deep North Atlantic Ocean. This large underestimate for these regions implies a factor of 5-6 anthropogenic enhancement in Hg concentrations.

Observationally constrained modeling of sound in curved ocean internal waves: examination of deep ducting and surface ducting at short range.

PubMed

Duda, Timothy F; Lin, Ying-Tsong; Reeder, D Benjamin

2011-09-01

A study of 400 Hz sound focusing and ducting effects in a packet of curved nonlinear internal waves in shallow water is presented. Sound propagation roughly along the crests of the waves is simulated with a three-dimensional parabolic equation computational code, and the results are compared to measured propagation along fixed 3 and 6 km source/receiver paths. The measurements were made on the shelf of the South China Sea northeast of Tung-Sha Island. Construction of the time-varying three-dimensional sound-speed fields used in the modeling simulations was guided by environmental data collected concurrently with the acoustic data. Computed three-dimensional propagation results compare well with field observations. The simulations allow identification of time-dependent sound forward scattering and ducting processes within the curved internal gravity waves. Strong acoustic intensity enhancement was observed during passage of high-amplitude nonlinear waves over the source/receiver paths, and is replicated in the model. The waves were typical of the region (35 m vertical displacement). Two types of ducting are found in the model, which occur asynchronously. One type is three-dimensional modal trapping in deep ducts within the wave crests (shallow thermocline zones). The second type is surface ducting within the wave troughs (deep thermocline zones). © 2011 Acoustical Society of America

Comparison of Shuttle Imaging Radar-B ocean wave image spectra with linear model predictions based on aircraft measurements

NASA Technical Reports Server (NTRS)

Monaldo, Frank M.; Lyzenga, David R.

1988-01-01

During October 1984, coincident Shuttle Imaging Radar-B synthetic aperture radar (SAR) imagery and wave measurements from airborne instrumentation were acquired. The two-dimensional wave spectrum was measured by both a radar ocean-wave spectrometer and a surface-contour radar aboard the aircraft. In this paper, two-dimensional SAR image intensity variance spectra are compared with these independent measures of ocean wave spectra to verify previously proposed models of the relationship between such SAR image spectra and ocean wave spectra. The results illustrate both the functional relationship between SAR image spectra and ocean wave spectra and the limitations imposed on the imaging of short-wavelength, azimuth-traveling waves.

Observations and simulations of microplastic marine debris in the ocean surface boundary layer

NASA Astrophysics Data System (ADS)

Kukulka, T.; Brunner, K.; Proskurowski, G. K.; Lavender Law, K. L.

2016-02-01

Motivated by observations of buoyant microplastic marine debris (MPMD) in the ocean surface boundary layer (OSBL), this study applies a large eddy simulation model and a parametric one-dimensional column model to examine vertical distributions of MPMD. MPMD is widely distributed in vast regions of the subtropical gyres and has emerged as a major open ocean pollutant whose distribution is subject to upper ocean turbulence. The models capture wind-driven turbulence, Langmuir turbulence (LT), and enhanced turbulent kinetic energy input due to breaking waves (BW). Model results are only consistent with MPMD observations if LT effects are included. Neither BW nor shear-driven turbulence is capable of deeply submerging MPMD, suggesting that the observed vertical MPMD distributions are a characteristic signature of wave-driven LT. Thus, this study demonstrates that LT substantially increases turbulent transport in the OSBL, resulting in deep submergence of buoyant tracers. The parametric model is applied to eleven years of observations in the North Atlantic and North Pacific subtropical gyres to show that surface measurements substantially underestimate MPMD concentrations by a factor of three to thirteen.

A spectral formalism for computing three-dimensional deformations due to surface loads. 2: Present-day glacial isostatic adjustment

NASA Technical Reports Server (NTRS)

Mitrovica, J. X.; Davis, J. L.; Shapiro, I. I.

1994-01-01

Using a spherically symmetric, self-gravitating, linear viscoelastic Earth model, we predict present-day three-dimensional surface deformation rates and baseline evolutions arising as a consequence of the late Pleistocene glacial cycles. In general, we use realistic models for the space-time geometry of the final late Pleistocene deglaciation event and incorporate a gravitationally self-consistent ocean meltwater redistribution. The predictions of horizontal velocity presented differ significantly, in both their amplitude and their spatial variation, from those presented in earlier analysis of others which adopted simplified models of both the late Pleistocene ice history and the Earth rheology. An important characteristic of our predicted velocity fields is that the melting of the Laurentide ice sheet over Canada is capable of contributing appreciably to the adjustment in Europe. The sensitivity of the predictions to variations in mantle rheology is investigated by considering a number of different Earth models, and by computing appropriate Frechet kernels. These calculations suggest that the sensitivity of the deformations to the Earth's rheology is significant and strongly dependent on the location of the site relative to the ancient ice sheet. The effects on the predictions of three-dimensional deformation rates of altering the ice history or adopting approximate models for the ocean meltwater redistribution have also been considered and found to be important (the former especially so). Finally, for a suite of Earth models we provide predictions of the velocity of a number of baselines in North America and Europe. We find that, in general, both radial and tangential motions contribute significantly to baseline length changes, and that these contributions are a strong function of the Earth model. We have, furthermore, found a set of Earth models which, together with the ICE-3G deglaciation chronology, produce predictions of baseline length changes that are consistent with very long baseline interferometry measurements of baselines within Europe.

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

NASA Technical Reports Server (NTRS)

Song, Tony Y.; Ji, Chen; Fu, L. -L.; Zlotnicki, Victor; Shum, C. K.; Yi, Yuchan; Hjorleifsdottir, Vala

2005-01-01

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

A world ocean model for greenhouse sensitivity studies: resolution intercomparison and the role of diagnostic forcing

NASA Astrophysics Data System (ADS)

Washington, Warren M.; Meehl, Gerald A.; Verplank, Lynda; Bettge, Thomas W.

1994-05-01

We have developed an improved version of a world ocean model with the intention of coupling to an atmospheric model. This article documents the simulation capability of this 1° global ocean model, shows improvements over our earlier 5° version, and compares it to features simulated with a 0.5° model. These experiments use a model spin-up methodology whereby the ocean model can subsequently be coupled to an atmospheric model and used for order 100-year coupled model integrations. With present-day computers, 1° is a reasonable compromise in resolution that allows for century-long coupled experiments. The 1° ocean model is derived from a 0.5°-resolution model developed by A. Semtner (Naval Postgraduate School) and R. Chervin (National Center for Atmospheric Research) for studies of the global eddy-resolving world ocean circulation. The 0.5° bottom topography and continental outlines have been altered to be compatible with the 1° resolution, and the Arctic Ocean has been added. We describe the ocean simulation characteristics of the 1° version and compare the result of weakly constraining (three-year time scale) the three-dimensional temperature and salinity fields to the observations below the thermocline (710 m) with the model forced only at the top of the ocean by observed annual mean wind stress, temperature, and salinity. The 1° simulations indicate that major ocean circulation patterns are greatly improved compared to the 5° version and are qualitatively reproduced in comparison to the 0.5° version. Using the annual mean top forcing alone in a 100-year simulation with the 1° version preserves the general features of the major observed temperature and salinity structure with most climate drift occurring mainly beneath the thermocline in the first 50 75 years. Because the thermohaline circulation in the 1° version is relatively weak with annual mean forcing, we demonstrate the importance of the seasonal cycle by performing two sensitivity experiments. Results show a dramatic intensification of the meridional overturning circulation (order of magnitude) with perpetual winter surface temperature forcing in the North Atlantic and strong intensification (factor of three) with perpetual early winter temperatures in that region. These effects are felt throughout the Atlantic (particularly an intensified and northward-shifted Gulf Stream outflow). In the Pacific, the temperature gradient strengthens in the thermocline, thus helping counter the systematic error of a thermocline that is too diffuse.

Pathways of PFOA to the Arctic: variabilities and contributions of oceanic currents and atmospheric transport and chemistry sources

NASA Astrophysics Data System (ADS)

Stemmler, I.; Lammel, G.

2010-10-01

Perfluorooctanoic acid (PFOA) and other perfluorinated compounds are industrial chemicals in use for decades which resist degradation in the environment and seem to accumulate in polar regions. Transport of PFOA was modeled using a spatially resolved global multicompartment model including fully coupled three-dimensional ocean and atmosphere general circulation models, and two-dimensional top soil, vegetation surfaces, and sea ice compartments. In addition to primary emissions, the formation of PFOA in the atmosphere from degradation of 8:2 fluorotelomer alcohol was included as a PFOA source. Oceanic transport, delivered 14.8±5.0 (8-23) t a-1 to the Arctic, strongly influenced by changes in water transport, which determined its interannual variability. This pathway constituted the dominant source of PFOA to the Arctic. Formation of PFOA in the atmosphere led to episodic transport events (timescale of days) into the Arctic with small spatial extent. Deposition in the polar region was found to be dominated by wet deposition over land, and shows maxima in boreal winter. The total atmospheric deposition of PFOA in the Arctic in the 1990s was ≈1 t a-1, much higher than previously estimated, and is dominated by primary emissions rather than secondary formation.

Pathways of PFOA to the Arctic: variabilities and contributions of oceanic currents and atmospheric transport and chemistry sources

NASA Astrophysics Data System (ADS)

Stemmler, I.; Lammel, G.

2010-05-01

Perfluorooctanoic acid (PFOA) and other perfluorinated compounds are industrial chemicals in use since decades which resist degradation in the environment and seem to accumulate in polar regions. Transport of PFOA was modeled using a spatially resolved global multicompartment model including fully coupled three-dimensional ocean and atmosphere general circulation models, and two-dimensional top soil, vegetation surfaces, and sea ice compartments. In addition to primary emissions, the formation of PFOA in the atmosphere from degradation of 8:2 fluorotelomer alcohol was included as a PFOA source. Oceanic transport, delivered 14.8±5.0 (8-23) t a-1 to the Arctic, strongly influenced by changes in water transport, which determined its interannual variability. This pathway constituted the dominant source of PFOA to the Arctic. Formation of PFOA in the atmosphere lead to episodic transport events (timescale of days) into the Arctic with small spatial extent. Deposition in the polar region was found to be dominated by wet deposition over land, and shows maxima in boreal winter. The total atmospheric deposition of PFOA in the Arctic in the 1990s was ≍1 t a-1, much higher than previously estimated, and is dominated by primary emissions rather than secondarily formed.

Forcing a three-dimensional, hydrostatic, primitive-equation model for application in the surf zone: 2. Application to DUCK94

NASA Astrophysics Data System (ADS)

Newberger, P. A.; Allen, J. S.

2007-08-01

A three-dimensional primitive-equation model for application to the nearshore surf zone has been developed. This model, an extension of the Princeton Ocean Model (POM), predicts the wave-averaged circulation forced by breaking waves. All of the features of the original POM are retained in the extended model so that applications can be made to regions where breaking waves, stratification, rotation, and wind stress make significant contributions to the flow behavior. In this study we examine the effects of breaking waves and wind stress. The nearshore POM circulation model is embedded within the NearCom community model and is coupled with a wave model. This combined modeling system is applied to the nearshore surf zone off Duck, North Carolina, during the DUCK94 field experiment of October 1994. Model results are compared to observations from this experiment, and the effects of parameter choices are examined. A process study examining the effects of tidal depth variation on depth-dependent wave-averaged currents is carried out. With identical offshore wave conditions and model parameters, the strength and spatial structure of the undertow and of the alongshore current vary systematically with water depth. Some three-dimensional solutions show the development of shear instabilities of the alongshore current. Inclusion of wave-current interactions makes an appreciable difference in the characteristics of the instability.

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.

Hierarchical algorithms for modeling the ocean on hierarchical architectures

NASA Astrophysics Data System (ADS)

Hill, C. N.

2012-12-01

This presentation will describe an approach to using accelerator/co-processor technology that maps hierarchical, multi-scale modeling techniques to an underlying hierarchical hardware architecture. The focus of this work is on making effective use of both CPU and accelerator/co-processor parts of a system, for large scale ocean modeling. In the work, a lower resolution basin scale ocean model is locally coupled to multiple, "embedded", limited area higher resolution sub-models. The higher resolution models execute on co-processor/accelerator hardware and do not interact directly with other sub-models. The lower resolution basin scale model executes on the system CPU(s). The result is a multi-scale algorithm that aligns with hardware designs in the co-processor/accelerator space. We demonstrate this approach being used to substitute explicit process models for standard parameterizations. Code for our sub-models is implemented through a generic abstraction layer, so that we can target multiple accelerator architectures with different programming environments. We will present two application and implementation examples. One uses the CUDA programming environment and targets GPU hardware. This example employs a simple non-hydrostatic two dimensional sub-model to represent vertical motion more accurately. The second example uses a highly threaded three-dimensional model at high resolution. This targets a MIC/Xeon Phi like environment and uses sub-models as a way to explicitly compute sub-mesoscale terms. In both cases the accelerator/co-processor capability provides extra compute cycles that allow improved model fidelity for little or no extra wall-clock time cost.

Potential utility of three-dimensional temperature and salinity fields estimated from satellite altimetry and Argo data for improving mesoscale reproducibility in regional ocean modeling

NASA Astrophysics Data System (ADS)

Kanki, R.; Uchiyama, Y.; Miyazaki, D.; Takano, A.; Miyazawa, Y.; Yamazaki, H.

2014-12-01

Mesoscale oceanic structure and variability are required to be reproduced as accurately as possible in realistic regional ocean modeling. Uchiyama et al. (2012) demonstrated with a submesoscale eddy-resolving JCOPE2-ROMS downscaling oceanic modeling system that the mesoscale reproducibility of the Kuroshio meandering along Japan is significantly improved by introducing a simple restoration to data which we call "TS nudging" (a.k.a. robust diagnosis) where the prognostic temperature and salinity fields are weakly nudged four-dimensionally towards the assimilative JCOPE2 reanalysis (Miyazawa et al., 2009). However, there is not always a reliable reanalysis for oceanic downscaling in an arbitrary region and at an arbitrary time, and therefore alternative dataset should be prepared. Takano et al. (2009) proposed an empirical method to estimate mesoscale 3-D thermal structure from the near real-time AVISO altimetry data along with the ARGO float data based on the two-layer model of Goni et al. (1996). In the present study, we consider the TS data derived from this method as a candidate. We thus conduct a synoptic forward modeling of the Kuroshio using the JCOPE2-ROMS downscaling system to explore potential utility of this empirical TS dataset (hereinafter TUM-TS) by carrying out two runs with the T-S nudging towards 1) the JCOPE2-TS and 2) TUM-TS fields.　An example of the comparison between the two ROMS test runs is shown in the attached figure showing the annually averaged surface EKE. Both of TUM-TS and JCOPE2-TS are found to help reproducing the mesoscale variance of the Koroshio and its extension as well as its mean paths, surface KE and EKE reasonably well. Therefore, the AVISO-ARGO derived empirical 3-D TS estimation is potentially exploitable for the dataset to conduct the T-S nudging to reproduce mesoscale oceanic structure.

Effect of 3-D heterogeneous-earth on rheology inference of postseismic model following the 2012 Indian Ocean earthquake

NASA Astrophysics Data System (ADS)

Pratama, C.; Ito, T.; Sasajima, R.; Tabei, T.; Kimata, F.; Gunawan, E.; Ohta, Y.; Yamashina, T.; Ismail, N.; Muksin, U.; Maulida, P.; Meilano, I.; Nurdin, I.; Sugiyanto, D.; Efendi, J.

2017-12-01

Postseismic deformation following the 2012 Indian Ocean earthquake has been modeled by several studies (Han et al. 2015, Hu et al. 2016, Masuti et al. 2016). Although each study used different method and dataset, the previous studies constructed a significant difference of earth structure. Han et al. (2015) ignored subducting slab beneath Sumatra while Masuti et al. (2016) neglect sphericity of the earth. Hu et al. (2016) incorporated elastic slab and spherical earth but used uniform rigidity in each layer of the model. As a result, Han et al. (2015) model estimated one order higher Maxwell viscosity than the Hu et al. (2016) and half order lower Kelvin viscosity than the Masuti et al. (2016) model predicted. In the present study, we conduct a quantitative analysis of each heterogeneous geometry and parameter effect on rheology inference. We develop heterogeneous three-dimensional spherical-earth finite element models. We investigate the effect of subducting slab, spherical earth, and three-dimensional earth rigidity on estimated lithosphere-asthenosphere rheology beneath the Indian Ocean. A wide range of viscosity structure from time constant rheology to time dependent rheology was chosen as previous studies have been modeled. In order to evaluate actual displacement, we compared the model to the Global Navigation Satellite System (GNSS) observation. We incorporate the GNSS data from previous studies and introduce new GNSS site as a part of the Indonesian Continuously Operating Reference Stations (InaCORS) located in Sumatra that has not been used in the last analysis. As a preliminary result, we obtained the effect of the spherical earth and elastic slab when we assumed burgers rheology. The model that incorporates the sphericity of the earth needs a one third order lower viscosity than the model that neglects earth curvature. The model that includes elastic slab needs half order lower viscosity than the model that excluding the elastic slab.

Three-Dimensional Spherical Models of Convection in the Earth's Mantle.

PubMed

Bercovici, D; Schubert, G; Glatzmaier, G A

1989-05-26

Three-dimensional, spherical models of mantle convection in the earth reveal that upwelling cylindrical plumes and downwelling planar sheets are the primary features of mantle circulation. Thus, subduction zones and descending sheetlike slabs in the mantle are fundamental characteristics of thermal convection in a spherical shell and are not merely the consequences of the rigidity of the slabs, which are cooler than the surrounding mantle. Cylindrical mantle plumes that cause hotspots such as Hawaii are probably the only form of active upwelling and are therefore not just secondary convective currents separate from the large-scale mantle circulation. Active sheetlike upwellings that could be associated with mid-ocean ridges did not develop in the model simulations, a result that is in agreement with evidence suggesting that ridges are passive phenomena resulting from the tearing of surface plates by the pull of descending slabs.

Sensitivity experiments with a one-dimensional coupled plume - iceflow model

NASA Astrophysics Data System (ADS)

Beckmann, Johanna; Perette, Mahé; Alexander, David; Calov, Reinhard; Ganopolski, Andrey

2016-04-01

Over the last few decades Greenland Ice sheet mass balance has become increasingly negative, caused by enhanced surface melting and speedup of the marine-terminating outlet glaciers at the ice sheet margins. Glaciers speedup has been related, among other factors, to enhanced submarine melting, which in turn is caused by warming of the surrounding ocean and less obviously, by increased subglacial discharge. While ice-ocean processes potentially play an important role in recent and future mass balance changes of the Greenland Ice Sheet, their physical understanding remains poorly understood. In this work we performed numerical experiments with a one-dimensional plume model coupled to a one-dimensional iceflow model. First we investigated the sensitivity of submarine melt rate to changes in ocean properties (ocean temperature and salinity), to the amount of subglacial discharge and to the glacier's tongue geometry itself. A second set of experiments investigates the response of the coupled model, i.e. the dynamical response of the outlet glacier to altered submarine melt, which results in new glacier geometry and updated melt rates.

Design of the MISMIP+, ISOMIP+, and MISOMIP ice-sheet, ocean, and coupled ice sheet-ocean intercomparison projects

NASA Astrophysics Data System (ADS)

Asay-Davis, Xylar; Cornford, Stephen; Martin, Daniel; Gudmundsson, Hilmar; Holland, David; Holland, Denise

2015-04-01

The MISMIP and MISMIP3D marine ice sheet model intercomparison exercises have become popular benchmarks, and several modeling groups have used them to show how their models compare to both analytical results and other models. Similarly, the ISOMIP (Ice Shelf-Ocean Model Intercomparison Project) experiments have acted as a proving ground for ocean models with sub-ice-shelf cavities.As coupled ice sheet-ocean models become available, an updated set of benchmark experiments is needed. To this end, we propose sequel experiments, MISMIP+ and ISOMIP+, with an end goal of coupling the two in a third intercomparison exercise, MISOMIP (the Marine Ice Sheet-Ocean Model Intercomparison Project). Like MISMIP3D, the MISMIP+ experiments take place in an idealized, three-dimensional setting and compare full 3D (Stokes) and reduced, hydrostatic models. Unlike the earlier exercises, the primary focus will be the response of models to sub-shelf melting. The chosen configuration features an ice shelf that experiences substantial lateral shear and buttresses the upstream ice, and so is well suited to melting experiments. Differences between the steady states of each model are minor compared to the response to melt-rate perturbations, reflecting typical real-world applications where parameters are chosen so that the initial states of all models tend to match observations. The three ISOMIP+ experiments have been designed to to make use of the same bedrock topography as MISMIP+ and using ice-shelf geometries from MISMIP+ results produced by the BISICLES ice-sheet model. The first two experiments use static ice-shelf geometries to simulate the evolution of ocean dynamics and resulting melt rates to a quasi-steady state when far-field forcing changes in either from cold to warm or from warm to cold states. The third experiment prescribes 200 years of dynamic ice-shelf geometry (with both retreating and advancing ice) based on a BISICLES simulation along with similar flips between warm and cold states in the far-field ocean forcing. The MISOMIP experiment combines the MISMIP+ experiments with the third ISOMIP+ experiment. Changes in far-field ocean forcing lead to a rapid (over ~1-2 years) increase in sub-ice-shelf melting, which is allowed to drive ice-shelf retreat for ~100 years. Then, the far-field forcing is switched to a cold state, leading to a rapid decrease in melting and a subsequent advance over ~100 years. To illustrate, we present results from BISICLES and POP2x experiments for each of the three intercomparison exercises.

Global isostatic geoid anomalies for plate and boundary layer models of the lithosphere

NASA Technical Reports Server (NTRS)

Hager, B. H.

1981-01-01

Commonly used one dimensional geoid models predict that the isostatic geoid anomaly over old ocean basins for the boundary layer thermal model of the lithosphere is a factor of two greater than that for the plate model. Calculations presented, using the spherical analogues of the plate and boundary layer thermal models, show that for the actual global distribution of plate ages, one dimensional models are not accurate and a spherical, fully three dimensional treatment is necessary. The maximum difference in geoid heights predicted for the two models is only about two meters. The thermal structure of old lithosphere is unlikely to be resolvable using global geoid anomalies. Stripping the effects of plate aging and a hypothetical uniform, 35 km, isostatically-compensated continental crust from the observed geoid emphasizes that the largest-amplitude geoid anomaly is the geoid low of almost 120 m over West Antarctica, a factor of two greater than the low of 60 m over Ceylon.

Three-dimensional shear wave velocity structure in the Atlantic upper mantle

NASA Astrophysics Data System (ADS)

James, Esther Kezia Candace

Oceanic lithosphere constitutes the upper boundary layer of the Earth's convecting mantle. Its structure and evolution provide a vital window on the dynamics of the mantle and important clues to how the motions of Earth's surface plates are coupled to convection in the mantle below. The three-dimensional shear-velocity structure of the upper mantle beneath the Atlantic Ocean is investigated to gain insight into processes that drive formation of oceanic lithosphere. Travel times are measured for approximately 10,000 fundamental-mode Rayleigh waves, in the period range 30-130 seconds, traversing the Atlantic basin. Paths with >30% of their length through continental upper mantle are excluded to maximize sensitivity to the oceanic upper mantle. The lateral distribution of Rayleigh wave phase velocity in the Atlantic upper mantle is explored with two approaches. One, phase velocity is allowed to vary only as a function of seafloor age. Two, a general two-dimensional parameterization is utilized in order to capture perturbations to age-dependent structure. Phase velocity shows a strong dependence on seafloor age, and removing age-dependent velocity from the 2-D maps highlights areas of anomalously low velocity, almost all of which are proximal to locations of hotspot volcanism. Depth-dependent variations in vertically-polarized shear velocity (Vsv) are determined with two sets of 3-D models: a layered model that requires constant VSV in each depth layer, and a splined model that allows VSV to vary continuously with depth. At shallow depths (˜75 km) the seismic structure shows the expected dependence on seafloor age. At greater depths (˜200 km) high-velocity lithosphere is found only beneath the oldest seafloor; velocity variations beneath younger seafloor may result from temperature or compositional variations within the asthenosphere. The age-dependent phase velocities are used to constrain temperature in the mantle and show that, in contrast to previous results for the Pacific, phase velocities for the Atlantic are not consistent with a half-space cooling model but are best explained by a plate-cooling model with thickness of 75 km and mantle temperature of 1400°C. Comparison with data such as basalt chemistry and seafloor elevation helps to separate thermal and compositional effects on shear velocity.

Implementation of the vortex force formalism in the coupled ocean-atmosphere-wave-sediment transport (COAWST) modeling system for inner shelf and surf zone applications

NASA Astrophysics Data System (ADS)

Kumar, Nirnimesh; Voulgaris, George; Warner, John C.; Olabarrieta, Maitane

The coupled ocean-atmosphere-wave-sediment transport modeling system (COAWST) enables simulations that integrate oceanic, atmospheric, wave and morphological processes in the coastal ocean. Within the modeling system, the three-dimensional ocean circulation module (ROMS) is coupled with the wave generation and propagation model (SWAN) to allow full integration of the effect of waves on circulation and vice versa. The existing wave-current coupling component utilizes a depth dependent radiation stress approach. In here we present a new approach that uses the vortex force formalism. The formulation adopted and the various parameterizations used in the model as well as their numerical implementation are presented in detail. The performance of the new system is examined through the presentation of four test cases. These include obliquely incident waves on a synthetic planar beach and a natural barred beach (DUCK' 94); normal incident waves on a nearshore barred morphology with rip channels; and wave-induced mean flows outside the surf zone at the Martha's Vineyard Coastal Observatory (MVCO). Model results from the planar beach case show good agreement with depth-averaged analytical solutions and with theoretical flow structures. Simulation results for the DUCK' 94 experiment agree closely with measured profiles of cross-shore and longshore velocity data from Garcez Faria et al. (1998, 2000). Diagnostic simulations showed that the nonlinear processes of wave roller generation and wave-induced mixing are important for the accurate simulation of surf zone flows. It is further recommended that a more realistic approach for determining the contribution of wave rollers and breaking induced turbulent mixing can be formulated using non-dimensional parameters which are functions of local wave parameters and the beach slope. Dominant terms in the cross-shore momentum balance are found to be the quasi-static pressure gradient and breaking acceleration. In the alongshore direction, bottom stress, breaking acceleration, horizontal advection and horizontal vortex forces dominate the momentum balance. The simulation results for the bar/rip channel morphology case clearly show the ability of the modeling system to reproduce horizontal and vertical circulation patterns similar to those found in laboratory studies and to numerical simulations using the radiation stress representation. The vortex force term is found to be more important at locations where strong flow vorticity interacts with the wave-induced Stokes flow field. Outside the surf zone, the three-dimensional model simulations of wave-induced flows for non-breaking waves closely agree with flow observations from MVCO, with the vertical structure of the simulated flow varying as a function of the vertical viscosity as demonstrated by Lentz et al. (2008).

Pathways of upwelling deep waters to the surface of the Southern Ocean

NASA Astrophysics Data System (ADS)

Tamsitt, Veronica; Drake, Henri; Morrison, Adele; Talley, Lynne; Dufour, Carolina; Gray, Alison; Griffies, Stephen; Mazloff, Matthew; Sarmiento, Jorge; Wang, Jinbo; Weijer, Wilbert

2017-04-01

Upwelling of Atlantic, Indian and Pacific deep waters to the sea surface in the Southern Ocean closes the global overturning circulation and is fundamentally important for oceanic uptake of anthropogenic carbon and heat, nutrient resupply for sustaining oceanic biological production, and the melt rate of ice shelves. Here we go beyond the two-dimensional view of Southern Ocean upwelling, to show detailed Southern Ocean upwelling pathways in three dimensions, using hydrographic observations and particle tracking in high-resolution ocean and climate models. The northern deep waters enter the Antarctic Circumpolar Current (ACC) via narrow southward currents along the boundaries of the three ocean basins, before spiraling southeastward and upward through the ACC. Upwelling is greatly enhanced at five major topographic features, associated with vigorous mesoscale eddy activity. Deep water reaches the upper ocean predominantly south of the southern ACC boundary, with a spatially nonuniform distribution, regionalizing warm water supply to Antarctic ice shelves and the delivery of nutrient and carbon-rich water to the sea surface. The timescale for half of the deep water to upwell from 30°S to the mixed layer is on the order of 60-90 years, which has important implications for the timescale for signals to propagate through the deep ocean. In addition, we quantify the diabatic transformation along particle trajectories, to identify where diabatic processes are important along the upwelling pathways.

Lagrangian water mass tracing from pseudo-Argo, model-derived salinity, tracer and velocity data: An application to Antarctic Intermediate Water in the South Atlantic Ocean

NASA Astrophysics Data System (ADS)

Blanke, Bruno; Speich, Sabrina; Rusciano, Emanuela

2015-01-01

We use the tracer and velocity fields of a climatological ocean model to investigate the ability of Argo-like data to estimate accurately water mass movements and transformations, in the style of analyses commonly applied to the output of ocean general circulation model. To this end, we introduce an algorithm for the reconstruction of a fully non-divergent three-dimensional velocity field from the simple knowledge of the model vertical density profiles and 1000-m horizontal velocity components. The validation of the technique consists in comparing the resulting pathways for Antarctic Intermediate Water in the South Atlantic Ocean to equivalent reference results based on the full model information available for velocity and tracers. We show that the inclusion of a wind-induced Ekman pumping and of a well-thought-out expression for vertical velocity at the level of the intermediate waters is essential for the reliable reproduction of quantitative Lagrangian analyses. Neglecting the seasonal variability of the velocity and tracer fields is not a significant source of errors, at least well below the permanent thermocline. These results give us confidence in the success of the adaptation of the algorithm to true gridded Argo data for investigating the dynamics of flows in the ocean interior.

Oceanographic and atmospheric conditions on the continental shelf north of the Monterey Bay during August 2006

NASA Astrophysics Data System (ADS)

Ramp, Steven R.; Lermusiaux, Pierre F. J.; Shulman, Igor; Chao, Yi; Wolf, Rebecca E.; Bahr, Frederick L.

2011-09-01

A comprehensive data set from the ocean and atmosphere was obtained just north of the Monterey Bay as part of the Monterey Bay 2006 (MB06) field experiment. The wind stress, heat fluxes, and sea surface temperature were sampled by the Naval Postgraduate School's TWIN OTTER research aircraft. In situ data were collected using ships, moorings, gliders and AUVs. Four data-assimilating numerical models were additionally run, including the Coupled Ocean/Atmosphere Mesoscale Prediction System (COAMPS ®) model for the atmosphere and the Harvard Ocean Prediction System (HOPS), the Regional Ocean Modeling System (ROMS), and the Navy Coastal Ocean Model (NCOM) for the ocean. The scientific focus of the Adaptive Sampling and Prediction Experiment (ASAP) was on the upwelling/relaxation cycle and the resulting three-dimensional coastal circulation near a coastal promontory, in this case Point Año Nuevo, CA. The emphasis of this study is on the circulation over the continental shelf as estimated from the wind forcing, two ADCP moorings, and model outputs. The wind stress during August 2006 consisted of 3-10 day upwelling favorable events separated by brief 1-3 day relaxations. During the first two weeks there was some correlation between local winds and currents and the three models' capability to reproduce the events. During the last two weeks, largely equatorward surface wind stress forced the sea surface and barotropic poleward flow occurred over the shelf, reducing model skill at predicting the circulation. The poleward flow was apparently remotely forced by mesoscale eddies and alongshore pressure gradients, which were not well simulated by the models. The small, high-resolution model domains were highly reliant on correct open boundary conditions to drive these larger-scale poleward flows. Multiply-nested models were no more effective than well-initialized local models in this respect.

TOPEX/Poseidon - An international satellite oceanography mission

NASA Technical Reports Server (NTRS)

Townsend, W. F.; Fellous, J.-L.

1986-01-01

The TOPEX/Poseidon mission, a joint NASA-CNES effort, strives to provide highly accurate global ocean topography measurements over a three year period utilizing highly advanced satellite radar altimetry techniques. Scheduled for launch in late 1991, the TOPEX/Poseidon satellite, together with ESA's first European remote sensing satellite and NASA's scatterometer, promises to provide a fundamental breakthrough in the present knowledge of how the oceans work as a global system. As part of the World Ocean Circulation Experiment, TOPEX/Poseidon measurements will aid in the determination of the three-dimensional current structure of the global oceans.

Ocean Circulation Modeling for Aquatic Dispersion of Liquid Radioactive Effluents from Nuclear Power Plants

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

Chung, Y.G.; Lee, G.B.; Bang, S.Y.

2006-07-01

Recently, three-dimensional models have been used for aquatic dispersion of radioactive effluents in relation to nuclear power plant siting based on the Notice No. 2003-12 'Guideline for investigating and assessing hydrological and aquatic characteristics of nuclear facility site' of the Ministry of Science and Technology (MOST) in Korea. Several nuclear power plants have been under construction or planed, which are Shin-Kori Unit 1 and 2, Shin-Wolsong Unit 1 and 2, and Shin-Ulchin Unit 1 and 2. For assessing the aquatic dispersion of radionuclides released from the above nuclear power plants, it is necessary to know the coastal currents around sitesmore » which are affected by circulation of East Sea. In this study, a three dimensional hydrodynamic model for the circulation of the East Sea of Korea has been developed as the first phase, which is based on the RIAMOM (Research Institute of Applied Mechanics' Ocean Model, Kyushu University, Japan). The model uses the primitive equation with hydrostatic approximation, and uses Arakawa-B grid system horizontally and Z coordinate vertically. Model domain is 126.5 deg. E to 142.5 deg. E of east longitude and 33 deg. N and 52 deg. N of the north latitude. The space of the horizontal grid was 1/12 deg. to longitude and latitude direction and vertical level was divided to 20. This model uses Generalized Arakawa Scheme, Slant Advection, and Mode-Splitting Method. The input data were from JODC (Japan Oceanographic Data Center), KNFRDI (Korea National Fisheries Research and Development Institute), and ECMWF (European Center for Medium-Range Weather Forecasts). The modeling results are in fairly good agreement with schematic patterns of the surface circulation in the East Sea/Japan Sea. The local current model and aquatic dispersion model of the coastal region will be developed as the second phase. The oceanic dispersion experiments will be also carried out by using ARGO Drifter around a nuclear power plant site. (authors)« less

An Oceanic General Circulation Model (OGCM) investigation of the Red Sea circulation: 2. Three-dimensional circulation in the Red Sea

NASA Astrophysics Data System (ADS)

Sofianos, Sarantis S.; Johns, William E.

2003-03-01

The three-dimensional circulation of the Red Sea is studied using a set of Miami Isopycnic Coordinate Ocean Model (MICOM) simulations. The model performance is tested against the few available observations in the basin and shows generally good agreement with the main observed features of the circulation. The main findings of this analysis include an intensification of the along-axis flow toward the coasts, with a transition from western intensified boundary flow in the south to eastern intensified flow in the north, and a series of strong seasonal or permanent eddy-like features. Model experiments conducted with different forcing fields (wind-stress forcing only, surface buoyancy forcing only, or both forcings combined) showed that the circulation produced by the buoyancy forcing is stronger overall and dominates the wind-driven part of the circulation. The main circulation pattern is related to the seasonal buoyancy flux (mostly due to the evaporation), which causes the density to increase northward in the basin and produces a northward surface pressure gradient associated with the downward sloping of the sea surface. The response of the eastern boundary to the associated mean cross-basin geostrophic current depends on the stratification and β-effect. In the northern part of the basin this results in an eastward intensification of the northward surface flow associated with the presence of Kelvin waves while in the south the traditional westward intensification due to Rossby waves takes place. The most prominent gyre circulation pattern occurs in the north where a permanent cyclonic gyre is present that is involved in the formation of Red Sea Outflow Water (RSOW). Beneath the surface boundary currents are similarly intensified southward undercurrents that carry the RSOW to the sill to flow out of the basin into the Indian Ocean.

Ocean products delivered by the Mercator Ocean Service Department

NASA Astrophysics Data System (ADS)

Crosnier, L.; Durand, E.; Soulat, F.; Messal, F.; Buarque, S.; Toumazou, V.; Landes, V.; Drevillon, M.; Lellouche, J.

2008-12-01

The newly created Service Department at Mercator Ocean is now offering various services for academic and private ocean applications. Mercator Ocean runs operationally ocean forecast systems for the Global and North Atlantic Ocean. These systems are based on an ocean general circulation model NEMO as well as on data assimilation of sea level anomalies, sea surface temperature and temperature and salinity vertical profiles. Three dimensional ocean fields of temperature, salinity and currents are updated and available weekly, including analysis and 2 weeks forecast fields. The Mercator Ocean service department is now offering a wide range of ocean derived products. This presentation will display some of the various products delivered in the framework of academic and private ocean applications: " Monitoring of the ocean current at the surface and at depth in several geographical areas for offshore oil platform, for offshore satellite launch platform, for transatlantic sailing or rowing boat races. " Monitoring of ocean climate indicators (Coral bleaching...) for marine reserve survey; " Monitoring of upwelling systems for fisheries; " Monitoring of the ocean heat content for tropical cyclone monitoring. " Monitoring of the ocean temperature/salinity and currents to guide research vessels during scientific cruises. The Mercator Ocean products catalogue will grow wider in the coming years, especially in the framework of the European GMES MyOcean project (FP7).

Experiments with a three-dimensional statistical objective analysis scheme using FGGE data

NASA Technical Reports Server (NTRS)

Baker, Wayman E.; Bloom, Stephen C.; Woollen, John S.; Nestler, Mark S.; Brin, Eugenia

1987-01-01

A three-dimensional (3D), multivariate, statistical objective analysis scheme (referred to as optimum interpolation or OI) has been developed for use in numerical weather prediction studies with the FGGE data. Some novel aspects of the present scheme include: (1) a multivariate surface analysis over the oceans, which employs an Ekman balance instead of the usual geostrophic relationship, to model the pressure-wind error cross correlations, and (2) the capability to use an error correlation function which is geographically dependent. A series of 4-day data assimilation experiments are conducted to examine the importance of some of the key features of the OI in terms of their effects on forecast skill, as well as to compare the forecast skill using the OI with that utilizing a successive correction method (SCM) of analysis developed earlier. For the three cases examined, the forecast skill is found to be rather insensitive to varying the error correlation function geographically. However, significant differences are noted between forecasts from a two-dimensional (2D) version of the OI and those from the 3D OI, with the 3D OI forecasts exhibiting better forecast skill. The 3D OI forecasts are also more accurate than those from the SCM initial conditions. The 3D OI with the multivariate oceanic surface analysis was found to produce forecasts which were slightly more accurate, on the average, than a univariate version.

One-Dimensional Coupled Ecosystem-Carbon Flux Model for the Simulation of Biogeochemical Parameters at Ocean Weather Station P

NASA Technical Reports Server (NTRS)

Signorini, S.; McClain, C.; Christian, J.; Wong, C. S.

2000-01-01

In this Technical Publication, we describe the model functionality and analyze its application to the seasonal and interannual variations of phytoplankton, nutrients, pCO2 and CO2 concentrations in the eastern subarctic Pacific at Ocean Weather Station P (OWSP, 50 deg. N 145 deg. W). We use a verified one-dimensional ecosystem model, coupled with newly incorporated carbon flux and carbon chemistry components, to simulate 22 years (1958-1980) of pCO2 and CO2 variability at Ocean Weather Station P (OWS P). This relatively long period of simulation verifies and extends the findings of previous studies using an explicit approach for the biological component and realistic coupling with the carbon flux dynamics. The slow currents and the horizontally homogeneous ocean in the subarctic Pacific make OWS P one of the best available candidates for modeling the chemistry of the upper ocean in one dimension. The chlorophyll and ocean currents composite for 1998 illustrates this premise. The chlorophyll concentration map was derived from SeaWiFS data and the currents are from an OGCM simulation (from R. Murtugudde).

Shallow Water Reverberation Measurement and Prediction

DTIC Science & Technology

1994-06-01

tool . The temporal signal processing consisted of a short-time Fourier transform spectral estimation method applied to data from a single hydrophone...The three-dimensional Hamiltonian Acoustic Ray-tracing Program for the Ocean (HARPO) was used as the primary propagation modeling tool . The temporal...summarizes the work completed and discusses lessons learned . Advice regarding future work to refine the present study will be provided. 6 our poiut source

Three-dimensional temperature fields of the North Patagonian Sea recorded by Magellanic penguins as biological sampling platforms

NASA Astrophysics Data System (ADS)

Sala, Juan E.; Pisoni, Juan P.; Quintana, Flavio

2017-04-01

Temperature is a primary determinant of biogeographic patterns and ecosystem processes. Standard techniques to study the ocean temperature in situ are, however, particularly limited by their time and spatial coverage, problems which might be partially mitigated by using marine top predators as biological platforms for oceanographic sampling. We used small archival tags deployed on 33 Magellanic penguins (Spheniscus magellanicus), and obtained 21,070 geo-localized profiles of water temperature, during late spring of 2008, 2011, 2012 and 2013; in a region of the North Patagonian Sea with limited oceanographic records in situ. We compared our in situ data of sea surface temperature (SST) with those available from satellite remote sensing; to describe the three-dimensional temperature fields around the area of influence of two important tidal frontal systems; and to study the inter-annual variation in the three-dimensional temperature fields. There was a strong positive relationship between satellite- and animal-derived SST data although there was an overestimation by remote-sensing by a maximum difference of +2 °C. Little inter-annual variability in the 3-dimensional temperature fields was found, with the exception of 2012 (and to a lesser extent in 2013) where the SST was significantly higher. In 2013, we found weak stratification in a region which was unexpected. In addition, during the same year, a warm small-scale vortex is indicated by the animal-derived temperature data. This allowed us to describe and better understand the dynamics of the water masses, which, so far, have been mainly studied by remote sensors and numerical models. Our results highlight again the potential of using marine top predators as biological platforms to collect oceanographic data, which will enhance and accelerate studies on the Southwest Atlantic Ocean. In a changing world, threatened by climate change, it is urgent to fill information gaps on the coupled ocean-atmosphere system allowing to link the hydrothermal process to the at-sea distribution of top predators.

Mesoscale Eddies in the Northwestern Pacific Ocean: Three-Dimensional Eddy Structures and Heat/Salt Transports

NASA Astrophysics Data System (ADS)

Dong, Di; Brandt, Peter; Chang, Ping; Schütte, Florian; Yang, Xiaofeng; Yan, Jinhui; Zeng, Jisheng

2017-12-01

The region encompassing the Kuroshio Extension (KE) in the Northwestern Pacific Ocean (25°N-45°N and 130°E-180°E) is one of the most eddy-energetic regions of the global ocean. The three-dimensional structures and transports of mesoscale eddies in this region are comprehensively investigated by combined use of satellite data and Argo profiles. With the allocation of Argo profiles inside detected eddies, the spatial variations of structures of eddy temperature and salinity anomalies are analyzed. The results show that eddies predominantly have subsurface (near-surface) intensified temperature and salinity anomalies south (north) of the KE jet, which is related to different background stratifications between these regions. A new method based on eddy trajectories and the inferred three-dimensional eddy structures is proposed to estimate heat and salt transports by eddy movements in a Lagrangian framework. Spatial distributions of eddy transports are presented over the vicinity of the KE for the first time. The magnitude of eddy-induced meridional heat (freshwater volume) transport is on the order of 0.01 PW (103 m3/s). The eddy heat transport divergence results in an oceanic heat loss south and heat gain north of the KE, thereby reinforcing and counteracting the oceanic heat loss from air-sea fluxes south and north of the KE jet, respectively. It also suggests a poleward heat transport across the KE jet due to eddy propagation.

Complementary constraints from carbon (13C) and nitrogen (15N) isotopes on the glacial ocean's soft-tissue biological pump

NASA Astrophysics Data System (ADS)

Schmittner, A.; Somes, C. J.

2016-06-01

A three-dimensional, process-based model of the ocean's carbon and nitrogen cycles, including 13C and 15N isotopes, is used to explore effects of idealized changes in the soft-tissue biological pump. Results are presented from one preindustrial control run (piCtrl) and six simulations of the Last Glacial Maximum (LGM) with increasing values of the spatially constant maximum phytoplankton growth rate μmax, which accelerates biological nutrient utilization mimicking iron fertilization. The default LGM simulation, without increasing μmax and with a shallower and weaker Atlantic Meridional Overturning Circulation and increased sea ice cover, leads to 280 Pg more respired organic carbon (Corg) storage in the deep ocean with respect to piCtrl. Dissolved oxygen concentrations in the colder glacial thermocline increase, which reduces water column denitrification and, with delay, nitrogen fixation, thus increasing the ocean's fixed nitrogen inventory and decreasing δ15NNO3 almost everywhere. This simulation already fits sediment reconstructions of carbon and nitrogen isotopes relatively well, but it overestimates deep ocean δ13CDIC and underestimates δ15NNO3 at high latitudes. Increasing μmax enhances Corg and lowers deep ocean δ13CDIC, improving the agreement with sediment data. In the model's Antarctic and North Pacific Oceans modest increases in μmax result in higher δ15NNO3 due to enhanced local nutrient utilization, improving the agreement with reconstructions there. Models with moderately increased μmax fit both isotope data best, whereas large increases in nutrient utilization are inconsistent with nitrogen isotopes although they still fit the carbon isotopes reasonably well. The best fitting models reproduce major features of the glacial δ13CDIC, δ15N, and oxygen reconstructions while simulating increased Corg by 510-670 Pg compared with the preindustrial ocean. These results are consistent with the idea that the soft-tissue pump was more efficient during the LGM. Both circulation and biological nutrient utilization could contribute. However, these conclusions are preliminary given our idealized experiments, which do not consider changes in benthic denitrification and spatially inhomogenous changes in aeolian iron fluxes. The analysis illustrates interactions between the carbon and nitrogen cycles as well as the complementary constraints provided by their isotopes. Whereas carbon isotopes are sensitive to circulation changes and indicate well the three-dimensional Corg distribution, nitrogen isotopes are more sensitive to biological nutrient utilization.

Estimating the Ocean Flow Field from Combined Sea Surface Temperature and Sea Surface Height Data

NASA Technical Reports Server (NTRS)

Stammer, Detlef; Lindstrom, Eric (Technical Monitor)

2002-01-01

This project was part of a previous grant at MIT that was moved over to the Scripps Institution of Oceanography (SIO) together with the principal investigator. The final report provided here is concerned only with the work performed at SIO since January 2000. The primary focus of this project was the study of the three-dimensional, absolute and time-evolving general circulation of the global ocean from a combined analysis of remotely sensed fields of sea surface temperature (SST) and sea surface height (SSH). The synthesis of those two fields was performed with other relevant physical data, and appropriate dynamical ocean models with emphasis on constraining ocean general circulation models by a combination of both SST and SSH data. The central goal of the project was to improve our understanding and modeling of the relationship between the SST and its variability to internal ocean dynamics, and the overlying atmosphere, and to explore the relative roles of air-sea fluxes and internal ocean dynamics in establishing anomalies in SST on annual and longer time scales. An understanding of those problems will feed into the general discussion on how SST anomalies vary with time and the extend to which they interact with the atmosphere.

Multidisciplinary research in the space sciences

NASA Technical Reports Server (NTRS)

Broecker, W. S.; Flynn, G. W.

1983-01-01

Research activities were carried out in the following areas during this reporting period: (1) astrophysics; (2) climate and atmospheric modeling; and (3) climate applications of earth observations & geological studies. An ultra-low-noise 115 GHz receiver based upon a superconducting tunnel diode mixer has been designed and constructed. The first laboratory tests have yielded spectacular results: a single-sideband noise temperature of 75 K considerably more sensitive than any other receiver at this frequency. The receiver will replace that currently in use on the Columbia-GISS CO Sky Survey telescope. The 1.2 meter millimeter-wave telescope at Columbia University has been used to complete two large-scale surveys of molecular matter in the part of the inner galaxy which is visible from the Northern hemisphere (the first galactic quadrant); one of the distant galaxy and one of the solar neighborhood. The research conducted during the past year in the climate and atmospheric modeling programs has been focused on the development of appropriate atmospheric and upper ocean models, and preliminary applications of these models. Principal models are a one-dimensional radiative-convective model, a three-dimensional global climate model, and an upper ocean model. During the past year this project has focused on development of 2-channel satellite analysis methods and radiative transfer studies in support of multichannel analysis techniques.

Three Dimensional Modeling Analysis of the Transpacific Transport of Aerosols During PACDEX

NASA Astrophysics Data System (ADS)

Carmichael, G. R.; Adhikary, B.; Hatch, C.; Kulkarni, S.; Moen, J.; Mena, M.

2007-12-01

Mineral dust and aerosols emitted from Asia are known to traverse long distances across the Pacific Ocean and can reach North America within a few days. A pilot field study, the PACific Dust Experiment (PACDEX), was carried out in April and May of 2007, during the peak East Asian dust emission season. The NSF/NCAR-HIAPER (High Performance Instrumented Airborne Platform for Environmental Research) platform allowed for sampling the evolution of mineral aerosol/pollution plumes and their physical and chemical characteristics as they traverse the Pacific Ocean and interact with the Pacific cloud systems en route to North America in both the upper and lower troposphere. A comprehensive 3-dimensional regional-scale model developed at The University of Iowa (Sulfur Transport dEposition Model, STEM) has been used for the analysis of aerosol interactions to help define key measurement strategies during the mission and to help interpret observations from the HIAPER platform. In this study we will present model aerosol distribution inter-comparison with cloud fields and aircraft observations. Model analysis provides further insight into cloud/pollution/dust interactions as East Asian emissions transit the Pacific Ocean en route to North America. Trajectory analysis and emission markers are used to help understand the air mass history and aerosol aging processes of the aerosols sampled by the HIAPER platform. Estimates of the fluxes of aerosol dust, BC and sulfate due to transpacific transport will also be presented.

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.

Three dimensional lithospheric magnetization structures beneath Australia derived by inverse modeling of CHAMP satellite magnetic field model

NASA Astrophysics Data System (ADS)

Du, Jinsong; Chen, Chao; Lesur, Vincent; Li, Yaoguo; Lane, Richard; Liang, Qing; Wang, Haoran

2014-05-01

We present an inversion algorithm for magnetic anomaly data in spherical coordinates to image the three dimensional (3-D) susceptibility distributions in the lithosphere. The method assumes that remanent magnetization is absent and that the magnetic anomalies are solely the result of lateral variations in magnetic susceptibility. To take into account the curvature of the Earth, the 3-D model is comprised of a set of spherical prisms (referred to as tesseroids), each of which has a constant isotropic susceptibility. The inversion method is formulated with a specifically designed model objective function and radial weighting function in spherical coordinates. A Tikhonov regularization technique is used to obtain an optimal solution with data misfit consistent with the estimated error level. Results for regional synthetic models with different magnetized inclinations and declinations are presented to demonstrate the capability of the method to recover large scale lithospheric magnetic structures. We have applied the algoithm to study the lithospheric susceptibility structures in the Australia region using magnetic anomaly data from the GRIMM_L120v0.0 model, which is based on ten years of CHAMP satellite data. As a self-constrained inversion, the maximum depths variation of magnetization layer is estimated first and then incorporated to the three dimensional (3-D) inversion. Results showed that the susceptibility anomalies concentrate in the depth range from 25 km to 45 km, i.e., focused in the lower crust. In addition, the results showed that the susceptibilities in continental lithosphere are higher than those in oceanic lithosphere. The inverted 3-D susceptibility distribution in the region of Australia reveals significant features related to tectonics, surface heat-flux, crustal thickness and Curie isotherm depths. In general, the higher susceptibility anomalies correlate with Precambrian rocks, and the lower susceptibility anomalies correlate with younger orogenic belts, suture zones and modern uplifts. In details, the inverted susceptibility distribution shows differences in the magnetic structures between the eastern and western parts of the Yilgarn Craton, and three lower susceptibility belts from north to south in the Eromanga Basin and the Gawler Craton with high susceptibility that extend to the ocean and then to the west.

On the relationship between satellite-estimated bio-optical and thermal properties in the Gulf of Mexico

NASA Astrophysics Data System (ADS)

Jolliff, Jason K.; Kindle, John C.; Penta, Bradley; Helber, Robert; Lee, Zhongping; Shulman, Igor; Arnone, Robert; Rowley, Clark D.

2008-03-01

Three years of Sea-viewing Wide Field-of-view Sensor (SeaWiFS) ocean color data were combined with three-dimensional thermal fields generated by the U.S. Navy's Modular Ocean Data Assimilation System (MODAS) in order to examine the interdependencies between bio-optical fields and their relationship to seasonal and mesoscale changes in upper ocean thermal structure. The combined data set suggests that the oceanic boundary layer within the Gulf of Mexico may be broadly defined by two seasonally occurring bio-thermal periods. A winter mixing period, characterized by net heat losses to the atmosphere, deepening of the isothermal layer depth, and annual maxima of satellite-estimated colored detrital matter (CDM) absorption coefficients and surface pigment concentration, was followed by a thermally stratified period characterized by net surface ocean heating, reduced isothermal layer depths, and annual minima in surface bio-optical fields. Variability in the interdependencies of ocean color products was used to diagnose an attendant shift in the size-structure of surface phytoplankton communities as well as identify CDM as the constituent responsible for the majority of blue-light absorption in Gulf of Mexico surface waters. The mesoscale circulation, as resolved by MODAS thermal fields into cold and warm-core eddies, appears to significantly modulate the seasonal bio-optical cycle of CDM absorption and surface pigment concentration. An empirical model was developed to describe CDM absorption as a function of upper ocean thermal energy. The model accounted for nearly half the variance in the satellite-estimate of this bio-optical variable. Large mismatches between the model and satellite data implied episodes of shelf water export to the deep Gulf of Mexico.

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.