Integrated approach to estimate the ocean's time variable dynamic topography including its covariance matrix

NASA Astrophysics Data System (ADS)

Müller, Silvia; Brockmann, Jan Martin; Schuh, Wolf-Dieter

2015-04-01

The ocean's dynamic topography as the difference between the sea surface and the geoid reflects many characteristics of the general ocean circulation. Consequently, it provides valuable information for evaluating or tuning ocean circulation models. The sea surface is directly observed by satellite radar altimetry while the geoid cannot be observed directly. The satellite-based gravity field determination requires different measurement principles (satellite-to-satellite tracking (e.g. GRACE), satellite-gravity-gradiometry (GOCE)). In addition, hydrographic measurements (salinity, temperature and pressure; near-surface velocities) provide information on the dynamic topography. The observation types have different representations and spatial as well as temporal resolutions. Therefore, the determination of the dynamic topography is not straightforward. Furthermore, the integration of the dynamic topography into ocean circulation models requires not only the dynamic topography itself but also its inverse covariance matrix on the ocean model grid. We developed a rigorous combination method in which the dynamic topography is parameterized in space as well as in time. The altimetric sea surface heights are expressed as a sum of geoid heights represented in terms of spherical harmonics and the dynamic topography parameterized by a finite element method which can be directly related to the particular ocean model grid. Besides the difficult task of combining altimetry data with a gravity field model, a major aspect is the consistent combination of satellite data and in-situ observations. The particular characteristics and the signal content of the different observations must be adequately considered requiring the introduction of auxiliary parameters. Within our model the individual observation groups are combined in terms of normal equations considering their full covariance information; i.e. a rigorous variance/covariance propagation from the original measurements to the final product is accomplished. In conclusion, the developed integrated approach allows for estimating the dynamic topography and its inverse covariance matrix on arbitrary grids in space and time. The inverse covariance matrix contains the appropriate weights for model-data misfits in least-squares ocean model inversions. The focus of this study is on the North Atlantic Ocean. We will present the conceptual design and dynamic topography estimates based on time variable data from seven satellite altimeter missions (Jason-1, Jason-2, Topex/Poseidon, Envisat, ERS-2, GFO, Cryosat2) in combination with the latest GOCE gravity field model and in-situ data from the Argo floats and near-surface drifting buoys.

Adaptation of a general circulation model to ocean dynamics

NASA Technical Reports Server (NTRS)

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

1976-01-01

A primitive-variable general circulation model of the ocean was formulated in which fast external gravity waves are suppressed with rigid-lid surface constraint pressires which also provide a means for simulating the effects of large-scale free-surface topography. The surface pressure method is simpler to apply than the conventional stream function models, and the resulting model can be applied to both global ocean and limited region situations. Strengths and weaknesses of the model are also presented.

Altimeter measurements for the determination of the Earth's gravity field

NASA Technical Reports Server (NTRS)

Tapley, B. D.; Schutz, B. E.; Shum, C. K.

1987-01-01

The ability of satellite-borne radar altimeter data to measure the global ocean surface with high precision and dense spatial coverage provides a unique tool for the mapping of the Earth's gravity field and its geoid. The altimeter crossover measurements, created by differencing direct altimeter measurements at the subsatellite points where the orbit ground tracks intersect, have the distinct advantage of eliminating geoid error and other nontemporal or long period oceanographic features. In the 1990's, the joint U.S./French TOPEX/POSEIDON mission and the European Space Agency's ERS-1 mission will carry radar altimeter instruments capable of global ocean mapping with high precision. This investigation aims at the development and application of dynamically consistent direct altimeter and altimeter crossover measurement models to the simultaneous mapping of the Earth's gravity field and its geoid, the ocean tides and the quasi-stationary component of the dynamic sea surface topography. Altimeter data collected by SEASAT, GEOS-3, and GEOSAT are used for the investigation.

Tracing Acoustic-Gravity Waves from the Ocean into the Ionosphere

NASA Astrophysics Data System (ADS)

Zabotin, N. A.; Godin, O. A.; Bullett, T. W.; Negrea, C.

2013-12-01

Ionospheric manifestations of tsunamis provide dramatic evidence of a connection between wave processes in the ocean and in the atmosphere. But tsunamis are only a transient feature of a more general phenomenon, infragravity waves (IGWs). IGWs are permanently present surface gravity waves in the ocean with periods longer than the longest periods (~30 s) of wind-generated waves. IGWs propagate transoceanic distances and, because of their long wavelengths (from ~1 km to hundreds of km), provide a mechanism for coupling wave processes in the ocean, atmosphere, and the solid Earth. The notion that tsunamis may generate waves in the upper atmosphere has existed for a long time but no quantitative coupling theory for the background waves has been proposed. We provide a strict physical justification for the influence of the background IGWs on the upper atmosphere. Taking into account both fluid compressibility and the gravity in a coupled atmosphere-ocean system, we show that there exist two distinct regimes of IGW penetration into the atmosphere. At higher frequencies, one has evanescent waves in the atmosphere propagating horizontally along the ocean surface. At lower frequencies, IGWs continuously radiate their energy into the upper atmosphere in the form of acoustic gravity waves (AGWs). The transition frequency depends on the ocean depth; it varies slowly near 3 mHz for typical depth values and drops to zero sharply only for extremely large depths. Using semi-empirical model of the IGW power spectrum, we derive an estimate of the flux of the mechanical energy and mechanical momentum from the deep ocean into the atmosphere due to background IGWs and predict specific forcing on the atmosphere in coastal regions. We compare spectra of wave processes in the ionosphere measured using Dynasonde technique over Wallops Island, VA and San Juan, PR and interpret the differences in terms of the oceanic effects. We conclude that AGWs of oceanic origin may have an observable impact on the upper atmosphere and describe techniques for experimental verification of this finding.

The Indian Ocean gravity low - Evidence for an isostatically uncompensated depression in the upper mantle

NASA Technical Reports Server (NTRS)

Ihnen, S. M.; Whitcomb, J. H.

1983-01-01

The broad gravity low in the equatorial Indian Ocean south of Sri Lanka is the largest and most striking feature in the gravitational field of the earth. The most negative long-wavelength free-air gravity anomalies are found there and the sea surface (geoid) lies more than 100 meters below the best fitting ellipsoid. A model of the lithosphere and upper mantle is proposed which accurately predicts the observed free-air gravity and geoid elevation. This model is consistent with bathymetry and sediment thickness data and suggests that the crust south of India currently floats as much as 600 meters lower than would be expected if the region were isostatically compensated. This residual depression of the crust is apparently confirmed by observations of ocean depth. An uncompensated depression is consistent with the presence of a mechanical wake left in the upper mantle behind India as it traveled toward Asia.

Global variations in gravity-derived oceanic crustal thickness: Implications on oceanic crustal accretion and hotspot-lithosphere interactions

NASA Astrophysics Data System (ADS)

Lin, J.; Zhu, J.

2012-12-01

We present a new global model of oceanic crustal thickness based on inversion of global oceanic gravity anomaly with constrains from seismic crustal thickness profiles. We first removed from the observed marine free-air gravity anomaly all gravitational effects that can be estimated and removed using independent constraints, including the effects of seafloor topography, marine sediment thickness, and the age-dependent thermal structure of the oceanic lithosphere. We then calculated models of gravity-derived crustal thickness through inversion of the residual mantle Bouguer anomaly using best-fitting gravity-modeling parameters obtained from comparison with seismically determined crustal thickness profiles. Modeling results show that about 5% of the global crustal volume (or 9% of the global oceanic surface area) is associated with model crustal thickness <5.2 km (designated as "thin" crust), while 56% of the crustal volume (or 65% of the surface area) is associated with crustal thickness of 5.2-8.6 km thick (designated as "normal" crust). The remaining 39% of the crustal volume (or 26% of the surface area) is associated with crustal thickness >8.6 km and is interpreted to have been affected by excess magmatism. The percentage of oceanic crustal volume that is associated with thick crustal thickness (>8.6 km) varies greatly among tectonic plates: Pacific (33%), Africa (50%), Antarctic (33%), Australia (30%), South America (34%), Nazca (23%), North America (47%), India (74%), Eurasia (68%), Cocos (20%), Philippine (26%), Scotia (41%), Caribbean (89%), Arabian (82%), and Juan de Fuca (21%). We also found that distribution of thickened oceanic crust (>8.6 km) seems to depend on spreading rate and lithospheric age: (1) On ocean basins younger than 5 Ma, regions of thickened crust are predominantly associated with slow and ultraslow spreading ridges. The relatively strong lithospheric plate at slow and ultraslow ridges might facilitate the loading of large magmatic emplacements on the plate. (2) In contrast, crustal thickness near fast and intermediately fast spreading ridges typically does not exceed 7-8 km. The relatively weak lithosphere at fast and intermediately fast ridges might make it harder for excess magmatism to accrete. We further speculate that the relatively wide partial melting zones in the upper mantle beneath the fast and intermediately fast ridges might act as "buffer" zones, thus diluting the melt anomalies from the underlying hotspots or regions of mantle heterogeneities. (3) As the crustal age increases and the lithospheric plate thickens, regions of thickened crust start to develop on ocean basins that were originally created at fast and intermediately fast ridges. The integrated crustal volume for fast and intermediately fast ocean crust appears to reach peak values for certain geological periods, such as 40-50 Ma and 70-80 Ma. The newly constructed global models of gravity-derived crustal thickness, combining with geochemical and other constraints, can be used to investigate the processes of oceanic crustal accretion and hotspot-lithosphere interactions.

Altimeter Observations of Baroclinic Oceanic Inertia-Gravity Wave Turbulence

NASA Technical Reports Server (NTRS)

Glazman, R. E.; Cheng, B.

1996-01-01

For a wide range of nonlinear wave processes - from capillary to planetary waves - theory predicts the existence of Kolmogorov-type spectral cascades of energy and other conserved quantities occuring via nonlinear resonant wave-wave interactions. So far, observations of wave turbulence (WT) have been limited to small-scale processes such as surface gravity and capillary-gravity waves.

How does ice sheet loading affect ocean flow around Antarctica?

NASA Astrophysics Data System (ADS)

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

2012-12-01

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

Surfing surface gravity waves

NASA Astrophysics Data System (ADS)

Pizzo, Nick

2017-11-01

A simple criterion for water particles to surf an underlying surface gravity wave is presented. It is found that particles travelling near the phase speed of the wave, in a geometrically confined region on the forward face of the crest, increase in speed. The criterion is derived using the equation of John (Commun. Pure Appl. Maths, vol. 6, 1953, pp. 497-503) for the motion of a zero-stress free surface under the action of gravity. As an example, a breaking water wave is theoretically and numerically examined. Implications for upper-ocean processes, for both shallow- and deep-water waves, are discussed.

On the generation of internal wave modes by surface waves

NASA Astrophysics Data System (ADS)

Harlander, Uwe; Kirschner, Ian; Maas, Christian; Zaussinger, Florian

2016-04-01

Internal gravity waves play an important role in the ocean since they transport energy and momentum and the can lead to mixing when they break. Surface waves and internal gravity waves can interact. On the one hand, long internal waves imply a slow varying shear current that modifies the propagation of surface waves. Surface waves generated by the atmosphere can, on the other hand, excite internal waves by nonlinear interaction. Thereby a surface wave packet consisting of two close frequencies can resonate with a low frequency internal wave (Phillips, 1966). From a theoretical point of view, the latter has been studied intensively by using a 2-layer model, i.e. a surface layer with a strong density contrast and an internal layer with a comparable weak density contrast (Ball, 1964; Craig et al., 2010). In the present work we analyse the wave coupling for a continuously stratified fluid using a fully non-linear 2D numerical model (OpenFoam) and compare this with laboratory experiments (see Lewis et al. 1974). Surface wave modes are used as initial condition and the time development of the dominant surface and internal waves are studied by spectral and harmonic analysis. For the simple geometry of a box, the results are compared with analytical spectra of surface and gravity waves. Ball, F.K. 1964: Energy transfer between external and internal gravity waves. J. Fluid Mech. 19, 465. Craig, W., Guyenne, P., Sulem, C. 2010: Coupling between internal and surface waves. Natural Hazards 57, 617-642. Lewis, J.E., Lake, B.M., Ko, D.R.S 1974: On the interaction of internal waves and surfacr gravity waves, J. Fluid Mech. 63, 773-800. Phillips, O.M. 1966: The dynamics of the upper ocean, Cambridge University Press, 336pp.

Application of precise altimetry to the study of precise leveling of the sea surface, the Earth's gravity field, and the rotation of the Earth

NASA Technical Reports Server (NTRS)

Segawa, J.; Ganeko, Y.; Sasaki, M.; Mori, T.; Ooe, M.; Nakagawa, I.; Ishii, H.; Hagiwara, Y.

1991-01-01

Our program includes five research items: (1) determination of a precision geoid and gravity anomaly field; (2) precise leveling and detection of tidal changes of the sea surface and study of the role of the tide in the global energy exchange; (3) oceanic effect on the Earth's rotation and polar motion; (4) geological and geophysical interpretation of the altimetry gravity field; and (5) evaluation of the effectiveness of local tracking of TOPEX/POSEIDON by use of a laser tracker.

Internal gravity wave contributions to global sea surface variability

NASA Astrophysics Data System (ADS)

Savage, A.; Arbic, B. K.; Richman, J. G.; Shriver, J. F.; Buijsman, M. C.; Zamudio, L.; Wallcraft, A. J.; Sharma, H.

2016-02-01

High-resolution (1/12th and 1/25th degree) 41-layer simulations of the HYbrid Coordinate Ocean Model (HYCOM), forced by both atmospheric fields and the astronomical tidal potential, are used to construct global maps of sea-surface height (SSH). The HYCOM output has been separated into steric, non-steric, and total sea-surface height and the maps display variance in subtidal, tidal, and supertidal bands. Two of the global maps are of particular interest in planning for the upcoming Surface Water and Ocean Topography (SWOT) wide-swath satellite altimeter mission; (1) a map of the nonstationary tidal signal (estimated after removing the stationary tidal signal via harmonic analysis), and (2) a map of the steric supertidal contributions, which are dominated by the internal gravity wave continuum. Both of these maps display signals of order 1 cm2, the target accuracy for the SWOT mission. Therefore, both non-stationary internal tides and non-tidal internal gravity waves are likely to be important sources of "noise" that must be accurately removed before examination of lower-frequency phenomena can take place.

Two-frequency /Delta k/ microwave scatterometer measurements of ocean wave spectra from an aircraft

NASA Technical Reports Server (NTRS)

Johnson, J. W.; Jones, W. L.; Weissman, D. E.

1981-01-01

A technique for remotely sensing the large-scale gravity wave spectrum on the ocean surface using a two frequency (Delta k) microwave scatterometer has been demonstrated from stationary platforms and proposed from moving platforms. This measurement takes advantage of Bragg type resonance matching between the electromagnetic wavelength at the difference frequency and the length of the large-scale surface waves. A prominent resonance appears in the cross product power spectral density (PSD) of the two backscattered signals. Ku-Band aircraft scatterometer measurements were conducted by NASA in the North Sea during the 1979 Maritime Remote Sensing (MARSEN) experiment. Typical examples of cross product PSD's computed from the MARSEN data are presented. They demonstrate strong resonances whose frequency and bandwidth agree with the surface characteristics and the theory. Directional modulation spectra of the surface reflectivity are compared to the gravity wave spectrum derived from surface truth measurements.

Detailed gravimetric geoid confirmation of short wavelength features of sea surface topography detected by the Skylab S-193 altimeter in the Atlantic Ocean

NASA Technical Reports Server (NTRS)

Marsh, J. G.; Vincent, S.; Mcclinton, A. T.; Chang, E. S.

1975-01-01

A detailed gravimetric geoid was computed for the Northwest Atlantic Ocean and Caribbean Sea area in support of the calibration and evaluation of the GEOS-C altimeter. This geoid, computed on a 15 ft. x 15 ft. grid was based upon a combination of surface gravity data with the GSFC GEM-6 satellite derived gravity data. A comparison of this gravimetric geoid with 10 passes of SKYLAB altimeter data is presented. The agreement of the two data types is quite good with the differences generally less than 2 meters. Sea surface manifestations of numerous short wavelength (approximately 100 km) oceanographic features are now indicated in the gravimetric geoid and are also confirmed by the altimetry data.

A numerical model for ocean ultra-low frequency noise: wave-generated acoustic-gravity and Rayleigh modes.

PubMed

Ardhuin, Fabrice; Lavanant, Thibaut; Obrebski, Mathias; Marié, Louis; Royer, Jean-Yves; d'Eu, Jean-François; Howe, Bruce M; Lukas, Roger; Aucan, Jerome

2013-10-01

The generation of ultra-low frequency acoustic noise (0.1 to 1 Hz) by the nonlinear interaction of ocean surface gravity waves is well established. More controversial are the quantitative theories that attempt to predict the recorded noise levels and their variability. Here a single theoretical framework is used to predict the noise level associated with propagating pseudo-Rayleigh modes and evanescent acoustic-gravity modes. The latter are dominant only within 200 m from the sea surface, in shallow or deep water. At depths larger than 500 m, the comparison of a numerical noise model with hydrophone records from two open-ocean sites near Hawaii and the Kerguelen islands reveal: (a) Deep ocean acoustic noise at frequencies 0.1 to 1 Hz is consistent with the Rayleigh wave theory, in which the presence of the ocean bottom amplifies the noise by 10 to 20 dB; (b) in agreement with previous results, the local maxima in the noise spectrum support the theoretical prediction for the vertical structure of acoustic modes; and (c) noise level and variability are well predicted for frequencies up to 0.4 Hz. Above 0.6 Hz, the model results are less accurate, probably due to the poor estimation of the directional properties of wind-waves with frequencies higher than 0.3 Hz.

3D Dynamics of the Near-Surface Layer of the Ocean in the Presence of Freshwater Influx

NASA Astrophysics Data System (ADS)

Dean, C.; Soloviev, A.

2015-12-01

Freshwater inflow due to convective rains or river runoff produces lenses of freshened water in the near surface layer of the ocean. These lenses are localized in space and typically involve both salinity and temperature anomalies. Due to significant density anomalies, strong pressure gradients develop, which result in lateral spreading of freshwater lenses in a form resembling gravity currents. Gravity currents inherently involve three-dimensional dynamics. The gravity current head can include the Kelvin-Helmholtz billows with vertical density inversions. In this work, we have conducted a series of numerical experiments using computational fluid dynamics tools. These numerical simulations were designed to elucidate the relationship between vertical mixing and horizontal advection of salinity under various environmental conditions and potential impact on the pollution transport including oil spills. The near-surface data from the field experiments in the Gulf of Mexico during the SCOPE experiment were available for validation of numerical simulations. In particular, we observed a freshwater layer within a few-meter depth range and, in some cases, a density inversion at the edge of the freshwater lens, which is consistent with the results of numerical simulations. In conclusion, we discuss applicability of these results to the interpretation of Aquarius and SMOS sea surface salinity satellite measurements. The results of this study indicate that 3D dynamics of the near-surface layer of the ocean are essential in the presence of freshwater inflow.

Long wavelength gravity and topography anomalies

NASA Technical Reports Server (NTRS)

Watts, A. B.; Daly, S. F.

1981-01-01

It is shown that gravity and topography anomalies on the earth's surface may provide new information about deep processes occurring in the earth, such as those associated with mantle convection. Two main reasons are cited for this. The first is the steady improvement that has occurred in the resolution of the long wavelength gravity field, particularly in the wavelength range of a few hundred to a few thousand km, mainly due to increased coverage of terrestrial gravity measurements and the development of radar altimeters in orbiting satellites. The second reason is the large number of numerical and laboratory experiments of convection in the earth, including some with deformable upper and lower boundaries and temperature-dependent viscosity. The oceans are thought to hold the most promise for determining long wavelength gravity and topography anomalies, since their evolution has been relatively simple in comparison with that of the continents. It is also shown that good correlation between long wavelength gravity and topography anomalies exists over some portions of the ocean floor

Near Field Ocean Surface Waves Acoustic Radiation Observation and Modeling

NASA Astrophysics Data System (ADS)

Ardhuin, F.; Peureux, C.; Royer, J. Y.

2016-12-01

The acoustic noise generation by nonlinearly interacting surface gravity waves has been studied for a long time both theoretically and experimentally [Longuet-Higgins 1951]. The associated far field noise is continuously measured by a vast network of seismometers at the ocean bottom and on the continents. It can especially be used to infer the time variability of short ocean waves statistics [Peureux and Ardhuin 2016]. However, better quantitative estimates of the latter are made difficult due to a poor knowledge of the Earth's crust characteristics, whose coupling with acoustic modes can affect large uncertainties to the frequency response at the bottom of the ocean.The pressure field at depths less than an acoustic wave length to the surface is made of evanescent modes which vanish away from their sources (near field) [Cox and Jacobs 1989]. For this reason, they are less affected by the ocean bottom composition. This near field is recorded and analyzed in the frequency range 0.1 to 0.5 Hz approximately, at two locations : at a shallow site in the North-East Atlantic continental shelf and a deep water site in the Southern Indian ocean, where pressure measurements are performed at the ocean bottom (ca. 100 m) and at 300 m water depth respectively. Evanescent and propagating Rayleigh modes are compared against theoretical predictions. Comparisons against surface waves hindcast based on WAVEWATCH(R) III modeling framework help assessing its performances and can be used to help future model improvements.References Longuet-Higgins, M. S., A Theory of the Origin of Microseisms, Philos. Trans. Royal Soc. A, 1950, 243, 1-3. Peureux, C. and Ardhuin, F., Ocean bottom pressure records from the Cascadia array and short surface gravity waves, J. Geophys. Res. Oceans, 2016, 121, 2862-2873. Cox, C. S. & Jacobs, D. C., Cartesian diver observations of double frequency pressure fluctuations in the upper levels of the ocean, Geophys. Res. Lett., 1989, 16, 807-810.

Numerical modeling of marine Gravity data for tsunami hazard zone mapping

NASA Astrophysics Data System (ADS)

Porwal, Nipun

2012-07-01

Tsunami is a series of ocean wave with very high wavelengths ranges from 10 to 500 km. Therefore tsunamis act as shallow water waves and hard to predict from various methods. Bottom Pressure Recorders of Poseidon class considered as a preeminent method to detect tsunami waves but Acoustic Modem in Ocean Bottom Pressure (OBP) sensors placed in the vicinity of trenches having depth of more than 6000m fails to propel OBP data to Surface Buoys. Therefore this paper is developed for numerical modeling of Gravity field coefficients from Bureau Gravimetric International (BGI) which do not play a central role in the study of geodesy, satellite orbit computation, & geophysics but by mathematical transformation of gravity field coefficients using Normalized Legendre Polynomial high resolution ocean bottom pressure (OBP) data is generated. Real time sea level monitored OBP data of 0.3° by 1° spatial resolution using Kalman filter (kf080) for past 10 years by Estimating the Circulation and Climate of the Ocean (ECCO) has been correlated with OBP data from gravity field coefficients which attribute a feasible study on future tsunami detection system from space and in identification of most suitable sites to place OBP sensors near deep trenches. The Levitus Climatological temperature and salinity are assimilated into the version of the MITGCM using the ad-joint method to obtain the sea height segment. Then TOPEX/Poseidon satellite altimeter, surface momentum, heat, and freshwater fluxes from NCEP reanalysis product and the dynamic ocean topography DOT_DNSCMSS08_EGM08 is used to interpret sea-bottom elevation. Then all datasets are associated under raster calculator in ArcGIS 9.3 using Boolean Intersection Algebra Method and proximal analysis tools with high resolution sea floor topographic map. Afterward tsunami prone area and suitable sites for set up of BPR as analyzed in this research is authenticated by using Passive microwave radiometry system for Tsunami Hazard Zone Mapping by network of seismometers. Thus using such methodology for early Tsunami Hazard Zone Mapping also increase accuracy and reduce time period for tsunami predictions. KEYWORDS:, Tsunami, Gravity Field Coefficients, Ocean Bottom Pressure, ECCO, BGI, Sea Bottom Temperature, Sea Floor Topography.

On the Direct Assimilation of Along-track Sea Surface Height Observations into a Free-surface Ocean Model Using a Weak Constraints Four Dimensional Variational (4dvar) Method

NASA Astrophysics Data System (ADS)

Ngodock, H.; Carrier, M.; Smith, S. R.; Souopgui, I.; Martin, P.; Jacobs, G. A.

2016-02-01

The representer method is adopted for solving a weak constraints 4dvar problem for the assimilation of ocean observations including along-track SSH, using a free surface ocean model. Direct 4dvar assimilation of SSH observations along the satellite tracks requires that the adjoint model be integrated with Dirac impulses on the right hand side of the adjoint equations for the surface elevation equation. The solution of this adjoint model will inevitably include surface gravity waves, and it constitutes the forcing for the tangent linear model (TLM) according to the representer method. This yields an analysis that is contaminated by gravity waves. A method for avoiding the generation of the surface gravity waves in the analysis is proposed in this study; it consists of removing the adjoint of the free surface from the right hand side (rhs) of the free surface mode in the TLM. The information from the SSH observations will still propagate to all other variables via the adjoint of the balance relationship between the barotropic and baroclinic modes, resulting in the correction to the surface elevation. Two assimilation experiments are carried out in the Gulf of Mexico: one with adjoint forcing included on the rhs of the TLM free surface equation, and the other without. Both analyses are evaluated against the assimilated SSH observations, SSH maps from Aviso and independent surface drifters, showing that the analysis that did not include adjoint forcing in the free surface is more accurate. This study shows that when a weak constraint 4dvar approach is considered for the assimilation of along-track SSH observations using a free surface model, with the aim of correcting the mesoscale circulation, an independent model error should not be assigned to the free surface.

Triad Resonance in the Gravity-Acoustic Family

NASA Astrophysics Data System (ADS)

Kadri, U.

2015-12-01

Resonance interactions of waves play a prominent role in energy share among the different wave types involved. Such interactions may significantly contribute, among others, to the evolution of the ocean energy spectrum by exchanging energy between surface-gravity waves; surface and internal gravity waves; or even surface and compression-type waves, that can transfer energy from the upper ocean through the whole water column reaching down to the seafloor. A resonant triad occurs among a triplet of waves, usually involving interaction of nonlinear terms of second order perturbed equations. Until recently, it has been believed that in a homogeneous fluid a resonant triad is possible only when tension forces are included, or at the limit of a shallow water, and that when the compressibility of water is considered, no resonant triads can occur within the family of gravity-acoustic waves. However, more recently it has been proved that, under some circumstances, resonant triads comprising two opposing surface-gravity waves of similar periods (though not identical) and a much longer acoustic-gravity wave, of almost double the frequency, exist [Kadri and Stiassnie 2013, J. Fluid Mech.735 R6]. Here, I report on a new resonant triad involving a gravity wave and two acoustic waves of almost double the length. Interestingly, the two acoustic waves propagate in the same direction with similar wavelengths, that are almost double of that of the gravity wave. The evolution of the wave triad amplitudes is periodic and it is derived analytically, in terms of Jacobian elliptic functions and elliptic integrals. The physical importance of this type of triad interactions is the modulation of pertinent acoustic signals, leading to inaccurate signal perceptions. Enclosed figure: presents an example spatio-temporal evolution of the wave triad amplitudes. The gravity wave (top) remains almost unaltered, while the envelope slowly displaces to the left. However, the prescribed acoustic envelope (middle) travels relatively fast to the right minimising the interaction time. Consequently, the resultant acoustic wave envelope (bottom) might be significantly smaller. As the two acoustic beams concurrently move away from the gravity wave, with disparate group velocities, the resonant interaction gradually vanishes.

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.

Basic research and data analysis for the earth and ocean physics applications program and for the National Geodetic Satellite Program

NASA Technical Reports Server (NTRS)

1976-01-01

Data analysis and supporting research in connection with the following objectives are discussed: (1) provide a precise and accurate geometric description of the earth's surface, (2) provide a precise and accurate mathematical description of the earth's gravitational field, and (3) determine time variations of the geometry of the ocean surface, the solid earth, the gravity field and other geophysical parameters.

Effect of Gravity Waves from Small Islands in the Southern Ocean on the Southern Hemisphere Atmospheric Circulation

NASA Technical Reports Server (NTRS)

Garfinkel, C. I.; Oman, L. D.

2018-01-01

The effect of small islands in the Southern Ocean on the atmospheric circulation in the Southern Hemisphere is considered with a series of simulations using the NASA Goddard Earth Observing System Chemistry-Climate Model in which the gravity wave stress generated by these islands is increased to resemble observed values. The enhanced gravity wave drag leads to a 2 K warming of the springtime polar stratosphere, partially ameliorating biases in this region. Resolved wave drag declines in the stratospheric region in which the added orographic gravity waves deposit their momentum, such that changes in gravity waves are partially compensated by changes in resolved waves, though resolved wave drag increases further poleward. The orographic drag from these islands has impacts for surface climate, as biases in tropospheric jet position are also partially ameliorated. These results suggest that these small islands are likely contributing to the missing drag near 60 degrees S in the upper stratosphere evident in many data assimilation products.

Gravity anomaly detection: Apollo/Soyuz

NASA Technical Reports Server (NTRS)

Vonbun, F. O.; Kahn, W. D.; Bryan, J. W.; Schmid, P. E.; Wells, W. T.; Conrad, D. T.

1976-01-01

The Goddard Apollo-Soyuz Geodynamics Experiment is described. It was performed to demonstrate the feasibility of tracking and recovering high frequency components of the earth's gravity field by utilizing a synchronous orbiting tracking station such as ATS-6. Gravity anomalies of 5 MGLS or larger having wavelengths of 300 to 1000 kilometers on the earth's surface are important for geologic studies of the upper layers of the earth's crust. Short wavelength Earth's gravity anomalies were detected from space. Two prime areas of data collection were selected for the experiment: (1) the center of the African continent and (2) the Indian Ocean Depression centered at 5% north latitude and 75% east longitude. Preliminary results show that the detectability objective of the experiment was met in both areas as well as at several additional anomalous areas around the globe. Gravity anomalies of the Karakoram and Himalayan mountain ranges, ocean trenches, as well as the Diamantina Depth, can be seen. Maps outlining the anomalies discovered are shown.

CryoSat-2 altimetry derived Arctic bathymetry map: first results and validation

NASA Astrophysics Data System (ADS)

Andersen, O. B.; Abulaitijiang, A.; Cancet, M.; Knudsen, P.

2017-12-01

The Technical University of Denmark (DTU), DTU Space has been developing high quality high resolution gravity fields including the new highly accurate CryoSat-2 radar altimetry satellite data which extends the global coverage of altimetry data up to latitude 88°. With its exceptional Synthetic Aperture Radar (SAR) mode being operating throughout the Arctic Ocean, leads, i.e., the ocean surface heights, is used to retrieve the sea surface height with centimeter-level range precision. Combined with the long repeat cycle ( 369 days), i.e., dense cross-track coverage, the high-resolution Arctic marine gravity can be modelled using the CryoSat-2 altimetry. Further, the polar gap can be filled by the available ArcGP product, thus yielding the complete map of the Arctic bathymetry map. In this presentation, we will make use of the most recent DTU17 marine gravity, to derive the arctic bathymetry map using inversion based on best available hydrographic maps. Through the support of ESA a recent evaluation of existing hydrographic models of the Arctic Ocean Bathymetry models (RTOPO, GEBCO, IBCAO etc) and various inconsistencies have been identified and means to rectify these inconsistencies have been taken prior to perform the inversion using altimetry. Simultaneously DTU Space has been placing great effort on the Arctic data screening, filtering, and de-noising using various altimetry retracking solutions and classifications. All the pre-processing contributed to the fine modelling of Actic gravity map. Thereafter, the arctic marine gravity grids will eventually be translated (downward continuation operation) to a new altimetry enhanced Arctic bathymetry map using appropriate band-pass filtering.

Surface topography due to convection in a variable viscosity fluid - Application to short wavelength gravity anomalies in the central Pacific Ocean

NASA Technical Reports Server (NTRS)

Lin, J.; Parmentier, E. M.

1985-01-01

Finite difference calculations of thermal convection in a fluid layer with a viscosity exponentially decreasing with temperature are performed in the context of examining the topography and gravity anomalies due to mantle convection. The surface topography and gravity anomalies are shown to be positive over regions of ascending flow and negative over regions of descending flow; at large Rayleigh numbers the amplitude of surface topography is inferred to depend on Rayleigh number to the power of 7/9. Compositional stratifications of the mantle is proposed as a mechanism for confining small-scale convection to a thin layer. A comparative analysis of the results with other available models is included.

The observation of ocean surface phenomena using imagery from the SEASAT synthetic aperture radar: An assessment

NASA Astrophysics Data System (ADS)

Vesecky, John F.; Stewart, Robert H.

1982-04-01

Over the period July 4 to October 10, 1978, the SEASAT synthetic aperture radar (SAR) gathered 23 cm wavelength radar images of some 108 km2 of the earth's surface, mainly of ocean areas, at 25-40 m resolution. Our assessment is in terms of oceanographic and ocean monitoring objectives and is directed toward discovering the proper role of SAR imagery in these areas of interest. In general, SAR appears to have two major and somewhat overlapping roles: first, quantitative measurement of ocean phenomena, like long gravity waves and wind fields, as well as measurement of ships; second, exploratory observations of large-scale ocean phenomena, such as the Gulf Stream and its eddies, internal waves, and ocean fronts. These roles are greatly enhanced by the ability of 23 cm SAR to operate day or night and through clouds. To begin we review some basics of synthetic aperture radar and its implementation on the SEASAT spacecraft. SEASAT SAR imagery of the ocean is fundamentally a map of the radar scattering characteristics of ˜30 cm wavelength ocean waves, distorted in some cases by ocean surface motion. We discuss how wind stress, surface currents, long gravity waves, and surface films modulate the scattering properties of these resonant waves with particular emphasis on the mechanisms that could produce images of long gravity waves. Doppler effects by ocean motion are also briefly described. Measurements of long (wavelength ≳100 m) gravity waves, using SEASAT SAR imagery, are compared with surface measurements during several experiments. Combining these results we find that dominant wavelength and direction are measured by SEASAT SAR within ±12% and ±15°, respectively. However, we note that ocean waves are not always visible in SAR images and discuss detection criteria in terms of wave height, length, and direction. SAR estimates of omnidirectional wave height spectra made by assuming that SAR image intensity is proportional to surface height fluctuations are more similar to corresponding surface measurements of wave height spectra than to wave slope spectra. Because SEASAT SAR images show the radar cross section σ° of ˜30 cm waves (neglecting doppler effects), and because these waves are raised by wind stress on the ocean surface, wind measurements are possible. Comparison between wind speeds estimated from SEASAT SAR imagery and from the SEASAT satellite scatterometer (SASS) agreed to within ±0.7 m s- over a 350-km comparison track and for wind speeds from 2 to 15 m s-. The great potential of SAR wind measurements lies in studying the spatial structure of the wind field over a range of spatial scales of from ≲1 km to ≳100 km. At present, the spatial and temporal structure of ocean wind fields is largely unknown. Because SAR responds to short waves whose energy density is a function of wind stress at the surface rather than wind speed at some distance above the surface, variations in image intensity may also reflect changes in air-sea temperature difference (thus complicating wind measurements by SAR). Because SAR images show the effects of surface current shear, air-sea temperature difference, and surface films through their modulation of the ˜30 cm waves, SEASAT images can be used to locate and study the Gulf Stream and related warm water rings, tidal flows at inlets, internal waves, and slicks resulting from surface films. In many of these applications, SAR provides a remote sensing capability that is complementary to infrared imagery because the two techniques sense largely different properties, namely, surface roughness and temperature. Both stationary ships and moving ships with their attendant wakes are often seen in SAR images. Ship images can be used to estimate ship size, heading, and speed. However, ships known to be in areas imaged by SAR are not always detectable. Clearly, a variety of factors, such as image resolution, ship size, sea state, and winds could affect ship detection. Overall, the role of SAR imagery in oceanography is definitely evolving at this time, but its ultimate role is unclear. We have assessed the ability of SEASAT SAR to measure a variety of ocean phenomena and have commented briefly on applications. In the end, oceanographers and others will have to judge from these capabilities the proper place for SAR in oceanography and remote sensing of the ocean.

Observation of `third sound' in superfluid 3He

NASA Astrophysics Data System (ADS)

Schechter, A. M. R.; Simmonds, R. W.; Packard, R. E.; Davis, J. C.

1998-12-01

Waves on the surface of a fluid provide a powerful tool for studying the fluid itself and the surrounding physical environment. For example, the wave speed is determined by the force per unit mass at the surface, and by the depth of the fluid: the decreasing speed of ocean waves as they approach the shore reveals the changing depth of the sea and the strength of gravity. Other examples include propagating waves in neutron-star oceans and on the surface of levitating liquid drops. Although gravity is a common restoring force, others exist, including the electrostatic force which causes a thin liquid film to adhere to a solid. Usually surface waves cannot occur on such thin films because viscosity inhibits their motion. However, in the special case of thin films of superfluid 4He, surface waves do exist and are called `third sound'. Here we report the detection of similar surface waves in thin films of superfluid 3He. We describe studies of the speed of these waves, the properties of the surface force, and the film's superfluid density.

Impact location of objects hitting the water surface

NASA Astrophysics Data System (ADS)

Kadri, Usama

2017-04-01

Analysis of data, recorded on March 8th 2014 at the Comprehensive Test ban Treaty Organisation's hydroacoustic station off Cape Leeuwin Western Australia, reveal pressure signatures of objects impacting at the sea surface which could be associated with falling meteorites as well as the missing Malaysian MH370 airplane. The location of the sources are identified analytically by an inverse solution based on acoustic-gravity wave theory (e.g. see references below) which have been developed and validated experimentally. Apart from the direct contribution to the search efforts after the missing airplane, the method we describe here is very efficient for identifying the location of sources that result in a sudden change in the water pressure in general. References 1. T.Yamamoto,1982.Gravity waves and acoustic waves generated by submarine earthquakes, Soil Dyn. Earthquake Eng., 1, 75-82. 2. M. Stiassnie, 2010. Tsunamis and acoustic-gravity waves from underwater earthquakes, J. Eng. Math., 67, 23-32, doi:10.1007/s10665-009-9323-x. 3. U. Kadri and M. Staissnie, 2012. Acoustic-gravity waves interacting with the shelf break. J. Geophys. Res., 117, C03035, doi: 10.1029/2011JC007674. 4. E. Eyov, A. Klar, U. Kadri and M. Stiassnie, 2013. Progressive waves in a compressible ocean with elastic bottom, Wave Motion 50, 929-939. doi: 10.1016/j.wavemoti.2013.03.003 5. G. Hendin and M. Stiassnie, 2013. Tsunami and acoustic-gravity waves in water of constant depth, Phys. Fluids 25, 086103, doi: 10.1063/1.481799. 6. U. Kadri, 2016. Acoustic-gravity waves from an oscillating ice-block in arctic zones. Advances in Acoustics and Vibration, 8076108, http://dx.doi.org/10.1155/2016/8076108 7. T.C.A. Oliveira, U. Kadri, 2016. Acoustic-gravity waves from the 2004 Indian Ocean earthquake and tsunami. Journal of Geophysical Research: Oceans. doi: 10.1002/2016JC011742

New Estimates of Hydrological and Oceanic Excitations of Variations of Earth's Rotation, Geocenter and Gravitational Field

NASA Technical Reports Server (NTRS)

Chao, Benjamin F.; Chen, J. L.; Johnson, T.; Au, A. Y.

1998-01-01

Hydrological mass transport in the geophysical fluids of the atmosphere-hydrosphere-solid Earth surface system can excite Earth's rotational variations in both length-of-day and polar motion. These effects can be computed in terms of the hydrological angular momentum by proper integration of global meteorological data. We do so using the 40-year NCEP data and the 18-year NASA GEOS-1 data, where the precipitation and evapotranspiration budgets are computed via the water mass balance of the atmosphere based on Oki et al.'s (1995) algorithm. This hydrological mass redistribution will also cause geocenter motion and changes in Earth's gravitational field, which are similarly computed using the same data sets. Corresponding geodynamic effects due to the oceanic mass transports (i.e. oceanic angular momentum and ocean-induced geocenter/gravity changes) have also been computed in a similar manner. We here compare two independent sets of the result from: (1) non-steric ocean surface topography observations based on Topex/Poseidon, and (2) the model output of the mass field by the Parallel Ocean Climate Model. Finally, the hydrological and the oceanic time series are combined in an effort to better explain the observed non-atmospheric effects. The latter are obtained by subtracting the atmospheric angular momentum from Earth rotation observations, and the atmosphere- induced geocenter/gravity effects from corresponding geodetic observations, both using the above-mentioned atmospheric data sets.

Ocean dynamics studies. [of current-wave interactions

NASA Technical Reports Server (NTRS)

1974-01-01

Both the theoretical and experimental investigations into current-wave interactions are discussed. The following three problems were studied: (1) the dispersive relation of a random gravity-capillary wave field; (2) the changes of the statistical properties of surface waves under the influence of currents; and (3) the interaction of capillary-gravity with the nonuniform currents. Wave current interaction was measured and the feasibility of using such measurements for remote sensing of surface currents was considered. A laser probe was developed to measure the surface statistics, and the possibility of using current-wave interaction as a means of current measurement was demonstrated.

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.

Characterization of the surface wave variability in the California Current region from satellite altimetry.

NASA Astrophysics Data System (ADS)

Villas Boas, A. B.; Gille, S. T.; Mazloff, M. R.

2016-02-01

Surface gravity waves play a crucial role in upper-ocean dynamics, and they are an important mechanism by which the ocean exchanges energy with the overlying atmosphere. Surface waves are largely wind forced and can also be modulated by ocean currents via nonlinear wave-current interactions, leading to either an amplification or attenuation of the wave amplitude. Even though individual waves cannot be detected by present satellite altimeters, surface waves have the potential to produce a sea-state bias in altimeter measurements and can impact the sea-surface-height spectrum at high wavenumbers or frequencies. Knowing the wave climatology is relevant for the success of future altimeter missions, such as the Surface Water and Ocean Topography (SWOT). We analyse the seasonal, intra-annual and interannual variability of significant wave heights retrieved from over two decades of satellite altimeter data and assess the extent to which the variability of the surface wave field in the California Current region is modulated by the local wind and current fields.

Comparison of undulation difference accuracies using gravity anomalies and gravity disturbances. [for ocean geoid

NASA Technical Reports Server (NTRS)

Jekeli, C.

1980-01-01

Errors in the outer zone contribution to oceanic undulation differences computed from a finite set of potential coefficients based on satellite measurements of gravity anomalies and gravity disturbances are analyzed. Equations are derived for the truncation errors resulting from the lack of high-degree coefficients and the commission errors arising from errors in the available lower-degree coefficients, and it is assumed that the inner zone (spherical cap) is sufficiently covered by surface gravity measurements in conjunction with altimetry or by gravity anomaly data. Numerical computations of error for various observational conditions reveal undulation difference errors ranging from 13 to 15 cm and from 6 to 36 cm in the cases of gravity anomaly and gravity disturbance data, respectively for a cap radius of 10 deg and mean anomalies accurate to 10 mgal, with a reduction of errors in both cases to less than 10 cm as mean anomaly accuracy is increased to 1 mgal. In the absence of a spherical cap, both cases yield error estimates of 68 cm for an accuracy of 1 mgal and between 93 and 160 cm for the lesser accuracy, which can be reduced to about 110 cm by the introduction of a perfect 30-deg reference field.

Feeling Gravity's Pull: Gravity Modeling. The Gravity Field of Mars

NASA Technical Reports Server (NTRS)

Lemoine, Frank; Smith, David; Rowlands, David; Zuber, Maria; Neumann, G.; Chinn, Douglas; Pavlis, D.

2000-01-01

Most people take the constant presence of gravitys pull for granted. However, the Earth's gravitational strength actually varies from location to location. This variation occurs because mass, which influences an object's gravitational pull, is not evenly distributed within the planet. Changes in topography, such as glacial movement, an earthquake, or a rise in the ocean level, can subtly affect the gravity field. An accurate measurement of the Earth's gravity field helps us understand the distribution of mass beneath the surface. This insight can assist us in locating petroleum, mineral deposits, ground water, and other valuable substances. Gravity mapping can also help notice or verify changes in sea surface height and other ocean characteristics. Such changes may indicate climate change from polar ice melting and other phenomena. In addition, gravity mapping can indicate how land moves under the surface after earthquakes and other plate tectonic processes. Finally, changes in the Earth's gravity field might indicate a shift in water distribution that could affect agriculture, water supplies for population centers, and long-term weather prediction. Scientists can map out the Earth's gravity field by watching satellite orbits. When a satellite shifts in vertical position, it might be passing over an area where gravity changes in strength. Gravity is only one factor that may shape a satellite's orbital path. To derive a gravity measurement from satellite movement, scientists must remove other factors that might affect a satellite's position: 1. Drag from atmospheric friction. 2. Pressure from solar radiation as it heads toward Earth and. as it is reflected off the surface of the Earth 3. Gravitational pull from the Sun, the Moon, and other planets in the Solar System. 4. The effect of tides. 5. Relativistic effects. Scientists must also correct for the satellite tracking process. For example, the tracking signal must be corrected for refraction through the atmosphere of the Earth. Supercomputers can calculate the effect of gravity for specific locations in space following a mathematical process known as spherical harmonics, which quantifies the gravity field of a planetary body. The process is based on Laplace's fundamental differential equation of gravity. The accuracy of a spherical harmonic solution is rated by its degree and order. Minute variations in gravity are measured against the geoid, a surface of constant gravity acceleration at mean sea level. The geoid reference gravity model strength includes the central body gravitational attraction (9.8 m/sq s) and a geopotential variation in latitude partially caused by the rotation of the Earth. The rotational effect modifies the shape of the geoid to be more like an ellipsoid, rather than a perfect, circle. Variations of gravity strength from the ellipsoidal reference model are measured in units called milli-Galileos (mGals). One mGal equals 10(exp -5) m/sq s. Research projects have also measured the gravity fields of other planetary bodies, as noted in the user profile that follows. From this information, we may make inferences about our own planet's internal structure and evolution. Moreover, mapping the gravity fields of other planets can help scientists plot the most fuel-efficient course for spacecraft expeditions to those planets.

Basic research and data analysis for the National Geodetic Satellite Program and for the Earth and Ocean Physics Application Program

NASA Technical Reports Server (NTRS)

1974-01-01

Accomplishments in the continuing programs are reported. The data were obtained in support of the following broad objectives: (1) to provide a precise and accurate geometric description of the earth's surface; (2) to provide a precise and accurate mathematical description of the earth's gravitational field; and (3) to determine time variations of the geometry of the ocean surface, the solid earth, the gravity field, and other geophysical parameters.

An Assessment of State-of-the-Art Mean Sea Surface and Geoid Models of the Arctic Ocean: Implications for Sea Ice Freeboard Retrieval

NASA Astrophysics Data System (ADS)

Skourup, Henriette; Farrell, Sinéad Louise; Hendricks, Stefan; Ricker, Robert; Armitage, Thomas W. K.; Ridout, Andy; Andersen, Ole Baltazar; Haas, Christian; Baker, Steven

2017-11-01

State-of-the-art Arctic Ocean mean sea surface (MSS) models and global geoid models (GGMs) are used to support sea ice freeboard estimation from satellite altimeters, as well as in oceanographic studies such as mapping sea level anomalies and mean dynamic ocean topography. However, errors in a given model in the high-frequency domain, primarily due to unresolved gravity features, can result in errors in the estimated along-track freeboard. These errors are exacerbated in areas with a sparse lead distribution in consolidated ice pack conditions. Additionally model errors can impact ocean geostrophic currents, derived from satellite altimeter data, while remaining biases in these models may impact longer-term, multisensor oceanographic time series of sea level change in the Arctic. This study focuses on an assessment of five state-of-the-art Arctic MSS models (UCL13/04 and DTU15/13/10) and a commonly used GGM (EGM2008). We describe errors due to unresolved gravity features, intersatellite biases, and remaining satellite orbit errors, and their impact on the derivation of sea ice freeboard. The latest MSS models, incorporating CryoSat-2 sea surface height measurements, show improved definition of gravity features, such as the Gakkel Ridge. The standard deviation between models ranges 0.03-0.25 m. The impact of remaining MSS/GGM errors on freeboard retrieval can reach several decimeters in parts of the Arctic. While the maximum observed freeboard difference found in the central Arctic was 0.59 m (UCL13 MSS minus EGM2008 GGM), the standard deviation in freeboard differences is 0.03-0.06 m.

Tsunami generation and associated waves in the water column and seabed due to an asymmetric earthquake motion within an anisotropic substratum

NASA Astrophysics Data System (ADS)

Bagheri, Amirhossein; Greenhalgh, Stewart; Khojasteh, Ali; Rahimian, Mohammad; Attarnejad, Reza

2016-10-01

In this paper, closed-form integral expressions are derived to describe how surface gravity waves (tsunamis) are generated when general asymmetric ground displacement (due to earthquake rupturing), involving both horizontal and vertical components of motion, occurs at arbitrary depth within the interior of an anisotropic subsea solid beneath the ocean. In addition, we compute the resultant hydrodynamic pressure within the seawater and the elastic wavefield within the seabed at any position. The method of potential functions and an integral transform approach, accompanied by a special contour integration scheme, are adopted to handle the equations of motion and produce the numerical results. The formulation accounts for any number of possible acoustic-gravity modes and is valid for both shallow and deep water situations as well as for any focal depth of the earthquake source. Phase and group velocity dispersion curves are developed for surface gravity (tsunami mode), acoustic-gravity, Rayleigh, and Scholte waves. Several asymptotic cases which arise from the general analysis are discussed and compared to existing solutions. The role of effective parameters such as hypocenter location and frequency of excitation is examined and illustrated through several figures which show the propagation pattern in the vertical and horizontal directions. Attention is directed to the unexpected contribution from the horizontal ground motion. The results have important application in several fields such as tsunami hazard prediction, marine seismology, and offshore and coastal engineering. In a companion paper, we examine the effect of ocean stratification on the appearance and character of internal and surface gravity waves.

Global accuracy estimates of point and mean undulation differences obtained from gravity disturbances, gravity anomalies and potential coefficients

NASA Technical Reports Server (NTRS)

Jekeli, C.

1979-01-01

Through the method of truncation functions, the oceanic geoid undulation is divided into two constituents: an inner zone contribution expressed as an integral of surface gravity disturbances over a spherical cap; and an outer zone contribution derived from a finite set of potential harmonic coefficients. Global, average error estimates are formulated for undulation differences, thereby providing accuracies for a relative geoid. The error analysis focuses on the outer zone contribution for which the potential coefficient errors are modeled. The method of computing undulations based on gravity disturbance data for the inner zone is compared to the similar, conventional method which presupposes gravity anomaly data within this zone.

Research activities of the Geodynamics Branch

NASA Technical Reports Server (NTRS)

Kahn, W. D. (Editor); Cohen, S. C. (Editor)

1984-01-01

A broad spectrum of geoscience disciplines including space geodesy, geopotential field modeling, tectonophysics, and dynamic oceanography are discussed. The NASA programs, include the Geodynamics and Ocean Programs, the Crustal Dynamics Project, the proposed Ocean Topography Experiment (TOPEX), and the Geopotential Research Mission (GRM). The papers are grouped into chapters on Crustal Movements, Global Earth Dynamics, Gravity Field Model Development, Sea Surface Topography, and Advanced Studies.

The relationship between surface topography, gravity anomalies, and temperature structure of convection

NASA Technical Reports Server (NTRS)

Parsons, B.; Daly, S.

1983-01-01

Consideration is given to the relationship between the temperature structure of mantle convection and the resulting surface topography and gravity anomalies, which are used in its investigation. Integral expressions relating the three variables as a function of wavelength are obtained with the use of Green's function solutions to the equations of motion for the case of constant-viscosity convection in a plane layer subject to a uniform gravitational field. The influence of the boundary conditions, particularly at large wavelengths, is pointed out, and surface topographies and gravity produced by convection are illustrated for a number of simple temperature distributions. It is shown that the upper thermal boundary layer plays an important role in determining the surface observables, while temperatures near the bottom of the layer affect mainly that boundary. This result is consistent with an explanation of geoid anomalies over mid-ocean swells in terms of convection beneath the lithosphere.

Regional Geoid Modeling Compared to Ocean Surface Observations

NASA Astrophysics Data System (ADS)

Roman, D. R.; Saleh, J.; Wang, Y. M.

2007-05-01

Aerogravity over a limited coastal region of the northern Gulf of Mexico enhanced and rectified the local gravity field signal. In turn, these data improved the derived geoid height model based on comparison with dynamic ocean topography (DOT) and tide gage information at eleven stations. Additionally, lidar observations were analyzed along nearly 50 profiles to estimate the reliability of these models into the offshore region. The overall comparison shows dm-level agreement between the various geoid and DOT models and ocean surface observations. An approximate 30 cm bias must still be explained; however, the results of this study point to the potential for further cooperative studies between oceanographers and geodesists.

Air-sea interaction with SSM/I and altimeter

NASA Technical Reports Server (NTRS)

1985-01-01

A number of important developments in satellite remote sensing techniques have occurred recently which offer the possibility of studying over vast areas of the ocean the temporally evolving energy exchange between the ocean and the atmosphere. Commencing in spring of 1985, passive and active microwave sensors that can provide valuable data for scientific utilization will start to become operational on Department of Defense (DOD) missions. The passive microwave radiometer can be used to estimate surface wind speed, total air column humidity, and rain rate. The active radar, or altimeter, senses surface gravity wave height and surface wind speed.

Spectral decomposition of internal gravity wave sea surface height in global models

NASA Astrophysics Data System (ADS)

Savage, Anna C.; Arbic, Brian K.; Alford, Matthew H.; Ansong, Joseph K.; Farrar, J. Thomas; Menemenlis, Dimitris; O'Rourke, Amanda K.; Richman, James G.; Shriver, Jay F.; Voet, Gunnar; Wallcraft, Alan J.; Zamudio, Luis

2017-10-01

Two global ocean models ranging in horizontal resolution from 1/12° to 1/48° are used to study the space and time scales of sea surface height (SSH) signals associated with internal gravity waves (IGWs). Frequency-horizontal wavenumber SSH spectral densities are computed over seven regions of the world ocean from two simulations of the HYbrid Coordinate Ocean Model (HYCOM) and three simulations of the Massachusetts Institute of Technology general circulation model (MITgcm). High wavenumber, high-frequency SSH variance follows the predicted IGW linear dispersion curves. The realism of high-frequency motions (>0.87 cpd) in the models is tested through comparison of the frequency spectral density of dynamic height variance computed from the highest-resolution runs of each model (1/25° HYCOM and 1/48° MITgcm) with dynamic height variance frequency spectral density computed from nine in situ profiling instruments. These high-frequency motions are of particular interest because of their contributions to the small-scale SSH variability that will be observed on a global scale in the upcoming Surface Water and Ocean Topography (SWOT) satellite altimetry mission. The variance at supertidal frequencies can be comparable to the tidal and low-frequency variance for high wavenumbers (length scales smaller than ˜50 km), especially in the higher-resolution simulations. In the highest-resolution simulations, the high-frequency variance can be greater than the low-frequency variance at these scales.

Displacements of the earth's surface due to atmospheric loading - Effects of gravity and baseline measurements

NASA Technical Reports Server (NTRS)

Van Dam, T. M.; Wahr, J. M.

1987-01-01

Atmospheric mass loads and deforms the earth's crust. By performing a convolution sum between daily, global barometric pressure data and mass loading Green's functions, the time dependent effects of atmospheric loading, including those associated with short-term synoptic storms, on surface point positioning measurements and surface gravity observations are estimated. The response for both an oceanless earth and an earth with an inverted barometer ocean is calculated. Load responses for near-coastal stations are significantly affected by the inclusion of an inverted barometer ocean. Peak-to-peak vertical displacements are frequently 15-20 mm with accompanying gravity perturbations of 3-6 micro Gal. Baseline changes can be as large as 20 mm or more. The perturbations are largest at higher latitudes and during winter months. These amplitudes are consistent with the results of Rabbel and Zschau (1985), who modeled synoptic pressure disturbances as Gaussian functions of radius around a central point. Deformation can be adequately computed using real pressure data from points within about 1000 km of the station. Knowledge of local pressure, alone, is not sufficient. Rabbel and Zschau's hypothesized corrections for these displacements, which use local pressure and the regionally averaged pressure, prove accurate at points well inland but are, in general, inadequate within a few hundred kilometers of the coast.

Crustal Thickness Mapping of the Rifted Margin Ocean-Continent Transition using Satellite Gravity Inversion Incorporating a Lithosphere Thermal Correction

NASA Astrophysics Data System (ADS)

Hurst, N. W.; Kusznir, N. J.

2005-05-01

A new method of inverting satellite gravity at rifted continental margins to give crustal thickness, incorporating a lithosphere thermal correction, has been developed which does not use a priori information about the location of the ocean-continent transition (OCT) and provides an independent prediction of OCT location. Satellite derived gravity anomaly data (Sandwell and Smith 1997) and bathymetry data (Gebco 2003) are used to derive the mantle residual gravity anomaly which is inverted in 3D in the spectral domain to give Moho depth. Oceanic lithosphere and stretched continental margin lithosphere produce a large negative residual thermal gravity anomaly (up to -380 mgal), which must be corrected for in order to determine Moho depth. This thermal gravity correction may be determined for oceanic lithosphere using oceanic isochron data, and for the thinned continental margin lithosphere using margin rift age and beta stretching estimates iteratively derived from crustal basement thickness determined from the gravity inversion. The gravity inversion using the thermal gravity correction predicts oceanic crustal thicknesses consistent with seismic observations, while that without the thermal correction predicts much too great oceanic crustal thicknesses. Predicted Moho depth and crustal thinning across the Hatton and Faroes rifted margins, using the gravity inversion with embedded thermal correction, compare well with those produced by wide-angle seismology. A new gravity inversion method has been developed in which no isochrons are used to define the thermal gravity correction. The new method assumes all lithosphere to be initially continental and a uniform lithosphere stretching age is used corresponding to the time of continental breakup. The thinning factor produced by the gravity inversion is used to predict the thickness of oceanic crust. This new modified form of gravity inversion with embedded thermal correction provides an improved estimate of rifted continental margin crustal thinning and an improved (and isochron independent) prediction of OCT location. The new method uses an empirical relationship to predict the thickness of oceanic crust as a function of lithosphere thinning factor controlled by two input parameters: a critical thinning factor for the start of ocean crust production and the maximum oceanic crustal thickness produced when the thinning factor = 1, corresponding to infinite lithosphere stretching. The disadvantage of using a uniform stretching age corresponding to the age of continental breakup is that the inversion fails to predict increasing thermal gravity correction towards the ocean ridge and incorrectly predicts thickening of oceanic crust with decreasing oceanic age. The new gravity inversion method has been applied to N. Atlantic rifted margins. This work forms part of the NERC Margins iSIMM project. iSIMM investigators are from Liverpool and Cambridge Universities, Badley Geoscience & Schlumberger Cambridge Research supported by the NERC, the DTI, Agip UK, BP, Amerada Hess Ltd, Anadarko, ConocoPhillips, Shell, Statoil and WesternGeco. The iSIMM team comprises NJ Kusznir, RS White, AM Roberts, PAF Christie, A Chappell, J Eccles, R Fletcher, D Healy, N Hurst, ZC Lunnon, CJ Parkin, AW Roberts, LK Smith, V Tymms & R Spitzer.

Continental and oceanic magnetic anomalies: Enhancement through GRM

NASA Technical Reports Server (NTRS)

Vonfrese, R. R. B.; Hinze, W. J.

1985-01-01

In contrast to the POGO and MAGSAT satellites, the Geopotential Research Mission (GRM) satellite system will orbit at a minimum elevation to provide significantly better resolved lithospheric magnetic anomalies for more detailed and improved geologic analysis. In addition, GRM will measure corresponding gravity anomalies to enhance our understanding of the gravity field for vast regions of the Earth which are largely inaccessible to more conventional surface mapping. Crustal studies will greatly benefit from the dual data sets as modeling has shown that lithospheric sources of long wavelength magnetic anomalies frequently involve density variations which may produce detectable gravity anomalies at satellite elevations. Furthermore, GRM will provide an important replication of lithospheric magnetic anomalies as an aid to identifying and extracting these anomalies from satellite magnetic measurements. The potential benefits to the study of the origin and characterization of the continents and oceans, that may result from the increased GRM resolution are examined.

Earth's gravity field to the eighteenth degree and geocentric coordinates for 104 stations from satellite and terrestrial data

NASA Technical Reports Server (NTRS)

Gaposchkin, E. M.

1973-01-01

Geodetic parameters describing the earth's gravity field and the positions of satellite-tracking stations in a geocentric reference frame were computed. These parameters were estimated by means of a combination of five different types of data: routine and simultaneous satellite observations, observations of deep-space probes, measurements of terrestrial gravity, and surface-triangulation data. The combination gives better parameters than does any subset of data types. The dynamic solution used precision-reduced Baker-Nunn observations and laser range data of 25 satellites. Data from the 49-station National Oceanic and Atmospheric Administration BC-4 network, the 19-station Smithsonian Astrophysical Observatory Baker-Nunn network, and independent camera stations were employed in the geometrical solution. Data from the tracking of deep-space probes were converted to relative longitudes and distances to the earth's axis of rotation of the tracking stations. Surface-gravity data in the form of 550-km squares were derived from 19,328 1 deg X 1 deg mean gravity anomalies.

Annual, Seasonal, and Secular Changes in Time-Variable Gravity from GRACE

NASA Astrophysics Data System (ADS)

Lemoine, F. G.; Luthcke, S. B.; Klosko, S. M.; Rowlands, D. D.; Chinn, D. S.; McCarthy, J. J.; Ray, R. D.; Boy, J.

2007-12-01

The NASA/DLR GRACE mission, launched in 2002, has now operated for more than five years, producing monthly and ten-day snapshots of the variations of the gravity field of the Earth. The available solutions, either from spherical harmonics or from mascons, allow us new insights into the variations of surface gravity on the Earth at annual, inter-annual, and secular time scales. Our baseline time series, based on GGM02C, NCEP Atmospheric Gravity with IB, and GOT00 tides now is extended to July 2007, spanning four+ years, and we analyze both mascon and spherical harmonic solutions from this time series with respect to global hydrology variations. Our 4degx4deg mascon solutions are extended to cover all continental regions of the globe. Comparisons with hydrology (land-surface) models can offer insights into how these models might be improved. We compare our baseline time series, with new time series that include an updated Goddard Ocean Tide (GOT) model, ECMWF- 3hr atmosphere de-aliasing data, and the MOG-2D ocean dealiasing product. Finally, we intercompare the spherical harmonic solutions at low degree from GRACE from the various product centers (e.g., GFZ, CSR, GRGS), and look for secular signals in both the GSFC mascon and spherical harmonic solutions, taking care to compare the results for secular gravity field change with independent solutions developed over 25 years of independent tracking to geodetic satellites by Satellite Laser Ranging (SLR) and DORIS.

Detailed geoid computations for GEOS-C altimeter experiment areas. [gravimetric geoid for Atlantic and northeast Pacific Ocean areas

NASA Technical Reports Server (NTRS)

Marsh, J. G.; Vincent, S.

1973-01-01

The GEOS-C spacecraft is scheduled to carry onboard a radar altimeter for the purpose of measuring the geoid undulations in oceanic areas. An independently derived geoid map will provide a valuable complement to these experiments. A detailed gravimetric geoid is presented for the Atlantic and northeast Pacific Ocean areas based upon a combination of the Goddard Space Flight Center GEM-6 earth model and surface 1 deg x 1 deg gravity data. As part of this work a number of satellite derived gravity models were evaluated to establish the model which best represented the long wave length features of the geoid in the above mentioned area. Comparisons of the detailed geoid with the astrogeodetic data provided by the National Ocean Survey and dynamically derived tracking station heights indicate that the accuracy of this combined geoid is on the order of 2 meters or better where data was dense and 5 to 7 meters where data was less dense.

A Century of Stability of Avannarleq and Kujalleq Glaciers, West Greenland, Explained Using High-Resolution Airborne Gravity and Other Data

NASA Astrophysics Data System (ADS)

An, L.; Rignot, E.; Mouginot, J.; Millan, R.

2018-04-01

The evolution of Greenland glaciers in a warming climate depends on their depth below sea level, flow speed, surface melt, and ocean-induced undercutting at the calving front. We present an innovative mapping of bed topography in the frontal regions of Sermeq Avannarleq and Kujalleq, two major glaciers flowing into the ice-choked Torssukatak Fjord, central west Greenland. The mapping combines a mass conservation algorithm inland, multibeam echo sounding data in the fjord, and high-resolution airborne gravity data at the ice-ocean transition where other approaches have traditionally failed. We obtain a reliable, precision (±40 m) solution for bed topography across the ice-ocean boundary. The results reveal a 700 m deep fjord that abruptly ends on a 100-300 m deep sill along the calving fronts. The shallow sills explain the presence of stranded icebergs, the resilience of the glaciers to ocean-induced undercutting by warm Atlantic water, and their remarkable stability over the past century.

Airborne Sea-Surface Topography in an Absolute Reference Frame

NASA Astrophysics Data System (ADS)

Brozena, J. M.; Childers, V. A.; Jacobs, G.; Blaha, J.

2003-12-01

Highly dynamic coastal ocean processes occur at temporal and spatial scales that cannot be captured by the present generation of satellite altimeters. Space-borne gravity missions such as GRACE also provide time-varying gravity and a geoidal msl reference surface at resolution that is too coarse for many coastal applications. The Naval Research Laboratory and the Naval Oceanographic Office have been testing the application of airborne measurement techniques, gravity and altimetry, to determine sea-surface height and height anomaly at the short scales required for littoral regions. We have developed a precise local gravimetric geoid over a test region in the northern Gulf of Mexico from historical gravity data and recent airborne gravity surveys. The local geoid provides a msl reference surface with a resolution of about 10-15 km and provides a means to connect airborne, satellite and tide-gage observations in an absolute (WGS-84) framework. A series of altimetry reflights over the region with time scales of 1 day to 1 year reveal a highly dynamic environment with coherent and rapidly varying sea-surface height anomalies. AXBT data collected at the same time show apparent correlation with wave-like temperature anomalies propagating up the continental slope of the Desoto Canyon. We present animations of the temporal evolution of the surface topography and water column temperature structure down to the 800 m depth of the AXBT sensors.

GOCE, Satellite Gravimetry and Antarctic Mass Transports

NASA Astrophysics Data System (ADS)

Rummel, Reiner; Horwath, Martin; Yi, Weiyong; Albertella, Alberta; Bosch, Wolfgang; Haagmans, Roger

2011-09-01

In 2009 the European Space Agency satellite mission GOCE (Gravity Field and Steady-State Ocean Circulation Explorer) was launched. Its objectives are the precise and detailed determination of the Earth's gravity field and geoid. Its core instrument, a three axis gravitational gradiometer, measures the gravity gradient components V xx , V yy , V zz and V xz (second-order derivatives of the gravity potential V) with high precision and V xy , V yz with low precision, all in the instrument reference frame. The long wavelength gravity field is recovered from the orbit, measured by GPS (Global Positioning System). Characteristic elements of the mission are precise star tracking, a Sun-synchronous and very low (260 km) orbit, angular control by magnetic torquing and an extremely stiff and thermally stable instrument environment. GOCE is complementary to GRACE (Gravity Recovery and Climate Experiment), another satellite gravity mission, launched in 2002. While GRACE is designed to measure temporal gravity variations, albeit with limited spatial resolution, GOCE is aiming at maximum spatial resolution, at the expense of accuracy at large spatial scales. Thus, GOCE will not provide temporal variations but is tailored to the recovery of the fine scales of the stationary field. GRACE is very successful in delivering time series of large-scale mass changes of the Antarctic ice sheet, among other things. Currently, emphasis of respective GRACE analyses is on regional refinement and on changes of temporal trends. One of the challenges is the separation of ice mass changes from glacial isostatic adjustment. Already from a few months of GOCE data, detailed gravity gradients can be recovered. They are presented here for the area of Antarctica. As one application, GOCE gravity gradients are an important addition to the sparse gravity data of Antarctica. They will help studies of the crustal and lithospheric field. A second area of application is ocean circulation. The geoid surface from the gravity field model GOCO01S allows us now to generate rather detailed maps of the mean dynamic ocean topography and of geostrophic flow velocities in the region of the Antarctic Circumpolar Current.

The absolute dynamic ocean topography (ADOT)

NASA Astrophysics Data System (ADS)

Bosch, Wolfgang; Savcenko, Roman

The sea surface slopes relative to the geoid (an equipotential surface) basically carry the in-formation on the absolute velocity field of the surface circulation. Pure oceanographic models may remain unspecific with respect to the absolute level of the ocean topography. In contrast, the geodetic approach to estimate the ocean topography as difference between sea level and the geoid gives by definition an absolute dynamic ocean topography (ADOT). This approach requires, however, a consistent treatment of geoid and sea surface heights, the first being usually derived from a band limited spherical harmonic series of the Earth gravity field and the second observed with much higher spectral resolution by satellite altimetry. The present contribution shows a procedure for estimating the ADOT along the altimeter profiles, preserving as much sea surface height details as the consistency w.r.t. the geoid heights will allow. The consistent treatment at data gaps and the coast is particular demanding and solved by a filter correction. The ADOT profiles are inspected for their innocent properties towards the coast and compared to external estimates of the ocean topography or the velocity field of the surface circulation as derived, for example, by ARGO floats.

3D Gravimetric Modeling of the Spreading System North and Southeast of the Rodriguez Triple Junction (Indian Ocean)

NASA Astrophysics Data System (ADS)

Heyde, I.; Girolami, C.; Barckhausen, U.; Freitag, R.

2017-12-01

Hydrothermal vent fields along mid-ocean ridges can be metal-rich and thus of great importance for the industries in the future. By order of the German Federal Ministry of Economics and in coordination with the International Seabed Authority (ISA), BGR explores potential areas of the active spreading system in the Indian Ocean. A main goal is the identification of inactive seafloor massive sulfides (SMS) with the aid of modern exploration techniques. Important contributions could be expected from bathymetric, magnetic, and gravity datasets, which can be acquired simultaneously time from the sea surface within relatively short ship time. The area of interest is located between 21°S and 28°S and includes the southern Central Indian Ridge (CIR) and the northern Southeast Indian Ridge (SEIR). In this study we analyzed the marine gravity and bathymetric data acquired during six research cruises. The profiles running perpendicular to the ridge axis have a mean length of 60 km. Magnetic studies reveal that the parts of the ridges covered are geologically very young with the oldest crust dating back to about 1 Ma. To extend the area outside the ridges, the shipboard data were complemented with data derived from satellite radar altimeter measurements. We analyzed the gravity anomalies along sections which cross particular geologic features (uplifted areas, accommodation zones, hydrothermal fields, and areas with hints for extensional processes e.g. oceanic core complexes) to establish a correlation between the gravity anomalies and the surface geology. Subsequently, for both ridge segments 3D density models were developed. We started with simple horizontally layered models, which, however, do not explain the measured anomalies satisfyingly. The density values of the crust and the upper mantle in the ridge areas had to be reduced. Finally, the models show the lateral heterogeneity and the variations in the thickness of the oceanic crust. There are areas characterized by crustal thickening related to magmatic accretion and areas of crustal thinning related to depleted accretion and exposure of OCCs.

Mapping Antarctic Crustal Thickness using Gravity Inversion and Comparison with Seismic Estimates

NASA Astrophysics Data System (ADS)

Kusznir, Nick; Ferraccioli, Fausto; Jordan, Tom

2017-04-01

Using gravity anomaly inversion, we produce comprehensive regional maps of crustal thickness and oceanic lithosphere distribution for Antarctica and the Southern Ocean. Crustal thicknesses derived from gravity inversion are consistent with seismic estimates. We determine Moho depth, crustal basement thickness, continental lithosphere thinning (1-1/β) and ocean-continent transition location using a 3D spectral domain gravity inversion method, which incorporates a lithosphere thermal gravity anomaly correction (Chappell & Kusznir 2008). The gravity anomaly contribution from ice thickness is included in the gravity inversion, as is the contribution from sediments which assumes a compaction controlled sediment density increase with depth. Data used in the gravity inversion are elevation and bathymetry, free-air gravity anomaly, the Bedmap 2 ice thickness and bedrock topography compilation south of 60 degrees south and relatively sparse constraints on sediment thickness. Ocean isochrons are used to define the cooling age of oceanic lithosphere. Crustal thicknesses from gravity inversion are compared with independent seismic estimates, which are still relatively sparse over Antarctica. Our gravity inversion study predicts thick crust (> 45 km) under interior East Antarctica, which is penetrated by narrow continental rifts featuring relatively thinner crust. The largest crustal thicknesses predicted from gravity inversion lie in the region of the Gamburtsev Subglacial Mountains, and are consistent with seismic estimates. The East Antarctic Rift System (EARS), a major Permian to Cretaceous age rift system, is imaged by our inversion and appears to extend from the continental margin at the Lambert Rift to the South Pole region, a distance of 2500 km. Offshore an extensive region of either thick oceanic crust or highly thinned continental crust lies adjacent to Oates Land and north Victoria Land, and also off West Antarctica around the Amundsen Ridges. Thin crust is predicted under the Ross Sea and beneath the West Antarctic Ice Sheet and delineates the regional extent of the broad West Antarctic Rift System (WARS). Substantial regional uplift is required under Marie Byrd Land to reconcile gravity and seismic estimates. A mantle dynamic uplift origin of the uplift is preferred to a thermal anomaly from a very young rift. The new maps produced by this study support the hypothesis that one branch of the WARS links through to the De Gerlache sea-mounts and Peter I Island in the Bellingshausen Sea region, while another branch may link to the George V Sound Rift in the Antarctic Peninsula region. Crustal thickness and lithosphere thinning derived from gravity inversion also allows the determination of circum-Antarctic ocean-continent transition structure and the mapping of continent-ocean boundary location. Superposition of illuminated satellite gravity data onto crustal thickness maps from gravity inversion provides improved determination of Southern Ocean rift orientation, pre-breakup rifted margin conjugacy and continental breakup trajectory. The continental lithosphere thinning distribution, used to define the initial thermal model temperature perturbation, is derived from the gravity inversion and uses no a priori isochron information; as a consequence the gravity inversion method provides a prediction of ocean-continent transition location, which is independent of ocean isochron information.

Can we observe the fronts of the Antarctic Circumpolar Current using GRACE OBP?

NASA Astrophysics Data System (ADS)

Makowski, J.; Chambers, D. P.; Bonin, J. A.

2014-12-01

The Antarctic Circumpolar Current (ACC) and the Southern Ocean remains one of the most undersampled regions of the world's oceans. The ACC is comprised of four major fronts: the Sub-Tropical Front (STF), the Polar Front (PF), the Sub-Antarctic Front (SAF), and the Southern ACC Front (SACCF). These were initially observed individually from repeat hydrographic sections and their approximate locations globally have been quantified using all available temperature data from the World Ocean and Climate Experiment (WOCE). More recent studies based on satellite altimetry have found that the front positions are more dynamic and have shifted south by up to 1° on average since 1993. Using ocean bottom pressure (OBP) data from the current Gravity Recovery and Climate Experiment (GRACE) we have measured integrated transport variability of the ACC south of Australia. However, differentiation of variability of specific fronts has been impossible due to the necessary smoothing required to reduce noise and correlated errors in the measurements. The future GRACE Follow-on (GFO) mission and the post 2020 GRACE-II mission are expected to produce higher resolution gravity fields with a monthly temporal resolution. Here, we study the resolution and error characteristics of GRACE gravity data that would be required to resolve variations in the front locations and transport. To do this, we utilize output from a high-resolution model of the Southern Ocean, hydrology models, and ice sheet surface mass balance models; add various amounts of random and correlated errors that may be expected from GFO and GRACE-II; and quantify requirements needed for future satellite gravity missions to resolve variations along the ACC fronts.

Anomalous Subsidence at Rifted Continental Margins: Distinguishing Mantle Dynamic Topography from Anomalous Oceanic Crustal Thickness

NASA Astrophysics Data System (ADS)

Cowie, L.; Kusznir, N. J.

2012-12-01

It has been proposed that some continental rifted margins have anomalous subsidence histories and that at breakup they were elevated at shallower bathymetries than the isostatic response of classical rift models (McKenzie 1978) would predict. The existence of anomalous syn or post breakup subsidence of this form would have important implications for our understanding of the geodynamics of continental breakup and rifted continental margin formation, margin subsidence history and the evolution of syn and post breakup depositional systems. We have investigated three rifted continental margins; the Gulf of Aden, Galicia Bank and the Gulf of Lions, to determine whether the oceanic crust in the ocean-continent transition of these margins has present day anomalous subsidence and if so, whether it is caused by mantle dynamic topography or anomalous oceanic crustal thickness. Residual depth anomalies (RDA) corrected for sediment loading, using flexural backstripping and decompaction, have been calculated by comparing observed and age predicted oceanic bathymetries in order to identify anomalous oceanic bathymetry and subsidence at these margins. Age predicted bathymetric anomalies have been calculated using the thermal plate model predictions from Crosby & McKenzie (2009). Non-zero sediment corrected RDAs may result from anomalous oceanic crustal thickness with respect to the global average, or from mantle dynamic uplift. Positive RDAs may result from thicker than average oceanic crust or mantle dynamic uplift; negative RDAs may result from thinner than average oceanic crust or mantle dynamic subsidence. Gravity inversion incorporating a lithosphere thermal gravity anomaly correction and sediment thickness from 2D seismic data has been used to determine Moho depth and oceanic crustal basement thickness. The reference Moho depths used in the gravity inversion have been calibrated against seismic refraction Moho depths. The gravity inversion crustal basement thicknesses together with Airy isostasy have been used to predict a "synthetic" gravity derived RDA. Sediment corrected RDA for oceanic crust in the Gulf of Aden are positive (+750m) indicating anomalous uplift with respect to normal subsidence. Gravity inversion predicts normal thickness oceanic crust and a zero "synthetic" gravity derived RDA in the oceanic domain. The difference between the positive sediment corrected RDA and the zero "synthetic" gravity derived RDA, implies that the anomalous subsidence reported in the Gulf of Aden is the result of mantle dynamic uplift. For the oceanic crust outboard of Galicia Bank both the sediment corrected RDA and the "synthetic" gravity derived RDA are negative (-800m) and of similar magnitude, indicating anomalous subsidence, which is the result of anomalously thin oceanic crust, not mantle dynamic topography. We conclude that there is negligible mantle dynamic topography influencing the Galicia Bank region. In the Gulf of Lions, gravity inversion predicts thinner than average oceanic crust. Both sediment corrected RDA (-1km) and "synthetic" gravity derived RDA (-500m) are negative. The more negative sediment corrected RDA compared with the "synthetic" gravity derived RDA implies that the anomalous subsidence in the Gulf of Lions is the result of mantle dynamic subsidence as well as thinner than average oceanic crust.

Measuring ocean waves from space; Proceedings of the Symposium, Johns Hopkins University, Laurel, MD, Apr. 15-17, 1986

NASA Technical Reports Server (NTRS)

Beal, Robert C. (Editor)

1987-01-01

Papers are presented on ocean-wave prediction; the quasi-universal form of the spectra of wind-generated gravity waves at different stages of their development; the limitations of the spectral measurements and observations of the group structure of surface waves; the effect of swell on the growth of wind wave; operational wave forecasting; ocean-wave models, and seakeeping using directional wave spectra. Consideration is given to microwave measurements of the ocean-wave directional spectra; SIR research; estimating wave energy spectra from SAR imagery, with the radar ocean-wave spectrometer, and SIR-B; the wave-measurement capabilities of the surface contour radar and the airborne oceanographic lidar; and SIR-B ocean-wave enhancement with fast-Fourier transform techniques. Topics discussed include wave-current interaction; the design and applicability of Spectrasat; the need for a global wave monitoring system; the age and source of ocean swell observed in Hurricane Josephine; and the use of satellite technology for insulin treatment.

Gravity field of the Western Weddell Sea: Comparison of airborne gravity and Geosat derived gravity

NASA Technical Reports Server (NTRS)

Bell, R. E.; Brozena, J. M.; Haxby, W. F.; Labrecque, J. L.

1989-01-01

Marine gravity surveying in polar regions was typically difficult and costly, requiring expensive long range research vessels and ice-breakers. Satellite altimetry can recover the gravity field in these regions where it is feasible to survey with a surface vessel. Unfortunately, the data collected by the first global altimetry mission, Seasat, was collected only during the austral winter, producing a very poor quality gravitational filed for the southern oceans, particularly in the circum-Antarctic regions. The advent of high quality airborne gravity (Brozena, 1984; Brozena and Peters, 1988; Bell, 1988) and the availability of satellite altimetry data during the austral summer (Sandwell and McAdoo, 1988) has allowed the recovery of a free air gravity field for most of the Weddell Sea. The derivation of the gravity field from both aircraft and satellite measurements are briefly reviewed, before presenting along track comparisons and shaded relief maps of the Weddell Sea gravity field based on these two data sets.

Evidence for Thin Oceanic Crust on the Extinct Aegir Ridge, Norwegian Basin, N.E. Atlantic Derived from Satellite Gravity Inversion

NASA Astrophysics Data System (ADS)

Greenhalgh, E. E.; Kusznir, N. J.

2006-12-01

Satellite gravity inversion incorporating a lithosphere thermal gravity correction has been used to map crustal thickness and lithosphere thinning factor for the N.E. Atlantic. The inversion of gravity data to determine crustal thickness incorporates a lithosphere thermal gravity anomaly correction for both oceanic and continental margin lithosphere. Predicted crustal thicknesses in the Norwegian Basin are between 7 and 4 km on the extinct Aegir oceanic ridge which ceased sea-floor spreading in the Oligocene. Crustal thickness estimates do not include a correction for sediment thickness and are upper bounds. Crustal thicknesses determined by gravity inversion for the Aegir Ridge are consistent with recent estimates derived using refraction seismology by Breivik et al. (2006). Failure to incorporate a lithosphere thermal gravity anomaly correction produces an over-estimate of crustal thickness. Oceanic crustal thicknesses within the Norwegian Basin are predicted by the gravity inversion to increase to 9-10 km eastwards towards the Norwegian (Moere) and westwards towards the Jan Mayen micro-continent, consistent with volcanic margin continental breakup at the end of the Palaeocene. The observation (from gravity inversion and seismic refraction studies) of thin oceanic crust produced by the Aegir ocean ridge in the Oligocene has implications for the temporal evolution of asthenosphere temperature under the N.E. Atlantic during the Tertiary. Thin Oligocene oceanic crust may imply cool (normal) asthenosphere temperatures during the Oligocene in contrast to elevated asthenosphere temperatures in the Palaeocene and Miocene-Recent as indicated by volcanic margin formation and the formation of Iceland respectively. Gravity inversion also predicts a region of thin oceanic crust to the west of the northern part of the Jan Mayen micro-continent and to the east of the thicker oceanic crust currently being formed at the Kolbeinsey Ridge. Thicker crust (c.f. ocean basins) is predicted for the Jan Mayen micro- continent south of Jan Mayen Island, with crust of the order of 20 km thickness extending southwards to connect with both the Faroes-Iceland Ridge and N.E. Iceland. Predicted crustal thicknesses under the Faroes- Iceland Ridge are approximately 25 km. The lithosphere thermal model used to predict the lithosphere thermal gravity anomaly correction may be conditioned using magnetic isochron data to provide the age of oceanic lithosphere. The resulting crustal thickness determination and the location of ocean-continent transition (OCT) are however sensitive to errors in the magnetic isochron data. An alternative method of inverting satellite gravity to give crustal thickness, incorporating a lithosphere thermal correction, has been used which does not use magnetic isochron data and provides an independent prediction of crustal thickness and OCT location. The crustal thickness estimates and OCT locations detailed above are robust to these sensitivity tests.

Basic research and data analysis for the earth and ocean physics applications program and for the National Geodetic Satellite program

NASA Technical Reports Server (NTRS)

1975-01-01

Data acquisition using single image and seven image data processing is used to provide a precise and accurate geometric description of the earth's surface. Transformation parameters and network distortions are determined, Sea slope along the continental boundaries of the U.S. and earth rotation are examined, along with close grid geodynamic satellite system. Data are derived for a mathematical description of the earth's gravitational field; time variations are determined for geometry of the ocean surface, the solid earth, gravity field, and other geophysical parameters.

Dynamic sea surface topography, gravity and improved orbit accuracies from the direct evaluation of SEASAT altimeter data

NASA Technical Reports Server (NTRS)

Marsh, J. G.; Lerch, F.; Koblinsky, C. J.; Klosko, S. M.; Robbins, J. W.; Williamson, R. G.; Patel, G. B.

1989-01-01

A method for the simultaneous solution of dynamic ocean topography, gravity and orbits using satellite altimeter data is described. A GEM-T1 based gravitational model called PGS-3337 that incorporates Seasat altimetry, surface gravimetry and satellite tracking data has been determined complete to degree and order 50. The altimeter data is utilized as a dynamic observation of the satellite's height above the sea surface with a degree 10 model of dynamic topography being recovered simultaneously with the orbit parameters, gravity and tidal terms in this model. PGS-3337 has a geoid uncertainty of 60 cm root-mean-square (RMS) globally, with the uncertainty over the altimeter tracked ocean being in the 25 cm range. Doppler determined orbits for Seasat, show large improvements, with the sub-30 cm radial accuracies being achieved. When altimeter data is used in orbit determination, radial orbital accuracies of 20 cm are achieved. The RMS of fit to the altimeter data directly gives 30 cm fits for Seasat when using PGS-3337 and its geoid and dynamic topography model. This performance level is two to three times better than that achieved with earlier Goddard earth models (GEM) using the dynamic topography from long-term oceanographic averages. The recovered dynamic topography reveals the global long wavelength circulation of the oceans with a resolution of 1500 km. The power in the dynamic topography recovery is now found to be closer to that of oceanographic studies than for previous satellite solutions. This is attributed primarily to the improved modeling of the geoid which has occurred. Study of the altimeter residuals reveals regions where tidal models are poor and sea state effects are major limitations.

The significance of ultra-refracted surface gravity waves on sheltered coasts, with application to San Francisco Bay

USGS Publications Warehouse

Hanes, D.M.; Erikson, L.H.

2013-01-01

Ocean surface gravity waves propagating over shallow bathymetry undergo spatial modification of propagation direction and energy density, commonly due to refraction and shoaling. If the bathymetric variations are significant the waves can undergo changes in their direction of propagation (relative to deepwater) greater than 90° over relatively short spatial scales. We refer to this phenomenon as ultra-refraction. Ultra-refracted swell waves can have a powerful influence on coastal areas that otherwise appear to be sheltered from ocean waves. Through a numerical modeling investigation it is shown that San Francisco Bay, one of the earth's largest and most protected natural harbors, is vulnerable to ultra-refracted ocean waves, particularly southwest incident swell. The flux of wave energy into San Francisco Bay results from wave transformation due to the bathymetry and orientation of the large ebb tidal delta, and deep, narrow channel through the Golden Gate. For example, ultra-refracted swell waves play a critical role in the intermittent closure of the entrance to Crissy Field Marsh, a small restored tidal wetland located on the sheltered north-facing coast approximately 1.5 km east of the Golden Gate Bridge.

Current state of art of satellite altimetry

NASA Astrophysics Data System (ADS)

Łyszkowicz, Adam Bolesław; Bernatowicz, Anna

2017-12-01

One of the fundamental problems of modern geodesy is precise defi nition of the gravitational fi eld and its changes in time. This is essential in positioning and navigation, geophysics, geodynamics, oceanography and other sciences related to the climate and Earth's environment. One of the major sources of gravity data is satellite altimetry that provides gravity data with almost 75% surface of the Earth. Satellite altimetry also provides data to study local, regional and global geophysical processes, the geoid model in the areas of oceans and seas. This technique can be successfully used to study the ocean mean dynamic topography. The results of the investigations and possible products of altimetry will provide a good material for the GGOS (Global Geodetic Observing System) and institutions of IAS (International Altimetry Service). This paper presents the achievements in satellite altimetry in all the above disciplines obtained in the last years. First very shorly basic concept of satellite altimetry is given. In order to obtain the highest accuracy on range measurements over the ocean improved of altimetry waveforms performed on the ground is described. Next, signifi cant improvements of sea and ocean gravity anomalies models developed presently is shown. Study of sea level and its extremes examined, around European and Australian coasts using tide gauges data and satellite altimetry measurements were described. Then investigations of the phenomenon of the ocean tides, calibration of altimeters, studies of rivers and ice-sheets in the last years are given.

A coupling modulation model of capillary waves from gravity waves: Theoretical analysis and experimental validation

NASA Astrophysics Data System (ADS)

Chen, Pengzhen; Wang, Xiaoqing; Liu, Li; Chong, Jinsong

2016-06-01

According to Bragg theory, capillary waves are the predominant scatterers of high-frequency band (such as Ka-band) microwave radiation from the surface of the ocean. Therefore, understanding the modulation mechanism of capillary waves is an important foundation for interpreting high-frequency microwave remote sensing images of the surface of the sea. In our experiments, we discovered that modulations of capillary waves are significantly larger than the values predicted by the classical theory. Further, analysis shows that the difference in restoring force results in an inflection point while the phase velocity changes from gravity waves region to capillary waves region, and this results in the capillary waves being able to resonate with gravity waves when the phase velocity of the gravity waves is equal to the group velocity of the capillary waves. Consequently, we propose a coupling modulation model in which the current modulates the capillary wave indirectly by modulating the resonant gravity waves, and the modulation of the former is approximated by that of the latter. This model very effectively explains the results discovered in our experiments. Further, based on Bragg scattering theory and this coupling modulation model, we simulate the modulation of normalized radar cross section (NRCS) of typical internal waves and show that the high-frequency bands are superior to the low-frequency bands because of their greater modulation of NRCS and better radiometric resolution. This result provides new support for choice of radar band for observation of wave-current modulation oceanic phenomena such as internal waves, fronts, and shears.

An improved dual-frequency technique for the remote sensing of ocean currents and wave spectra

NASA Technical Reports Server (NTRS)

Schuler, D. L.; Eng, W. P.

1984-01-01

A two frequency microwave radar technique for the remote sensing of directional ocean wave spectra and surface currents is investigated. This technique is conceptually attractive because its operational physical principle involves a spatial electromagnetic scattering resonance with a single, but selectable, long gravity wave. Multiplexing of signals having different spacing of the two transmitted frequencies allows measurements of the entire long wave ocean spectrum to be carried out. A new scatterometer is developed and experimentally tested which is capable of making measurements having much larger signal/background values than previously possible. This instrument couples the resonance technique with coherent, frequency agility radar capabilities. This scatterometer is presently configured for supporting a program of surface current measurements.

Development and Utilization of Regional Oceanic Modeling System (ROMS). Delicacy, Imprecision, and Uncertainty of Oceanic Simulations: An Investigation with the Regional Oceanic Modeling System (ROMS). Submesoscale Flows and Mixing in the Ocean Surface Layer Using the Regional Oceanic Modeling System (ROMS). Eddy Effects in General Circulation, Spanning Mean Currents, Mesoscale Eddies, and Topographic Generation, including Submesoscale Nests

DTIC Science & Technology

2012-09-30

unbalanced motions is likely to occur. Due to an rapidly expanding set of investigation on oceanic flows at submesoscales, it is increasingly clear...Uchiyama, E. M. Lane, J. M. Restrepo, & J. C. McWilliams, 2011: A vortex force analysis of the interaction of rip currents and gravity waves. J. Geophys...particular topographic features, the torque is pervasively positive (cyclonic) along the Stream, in opposition to the anticyclonic wind curl in the

Arctic and N Atlantic Crustal Thickness and Oceanic Lithosphere Distribution from Gravity Inversion

NASA Astrophysics Data System (ADS)

Kusznir, Nick; Alvey, Andy

2014-05-01

The ocean basins of the Arctic and N. Atlantic formed during the Mesozoic and Cenozoic as a series of distinct ocean basins, both small and large, leading to a complex distribution of oceanic crust, thinned continental crust and rifted continental margins. The plate tectonic framework of this region was demonstrated by the pioneering work of Peter Ziegler in AAPG Memoir 43 " Evolution of the Arctic-North Atlantic and the Western Tethys" published in 1988. The spatial evolution of Arctic Ocean and N Atlantic ocean basin geometry and bathymetry are critical not only for hydrocarbon exploration but also for understanding regional palaeo-oceanography and ocean gateway connectivity, and its influence on global climate. Mapping crustal thickness and oceanic lithosphere distribution represents a substantial challenge for the Polar Regions. Using gravity anomaly inversion we have produced comprehensive maps of crustal thickness and oceanic lithosphere distribution for the Arctic and N Atlantic region, We determine Moho depth, crustal basement thickness, continental lithosphere thinning and ocean-continent transition location using a 3D spectral domain gravity inversion method, which incorporates a lithosphere thermal gravity anomaly correction (Chappell & Kusznir 2008). Gravity anomaly and bathymetry data used in the gravity inversion are from the NGA (U) Arctic Gravity Project and IBCAO respectively; sediment thickness is from a new regional compilation. The resulting maps of crustal thickness and continental lithosphere thinning factor are used to determine continent-ocean boundary location and the distribution of oceanic lithosphere. Crustal cross-sections using Moho depth from the gravity inversion allow continent-ocean transition structure to be determined and magmatic type (magma poor, "normal" or magma rich). Our gravity inversion predicts thin crust and high continental lithosphere thinning factors in the Eurasia, Canada, Makarov, Podvodnikov and Baffin Basins consistent with these basins being oceanic. Larger crustal thicknesses, in the range 20 - 30 km, are predicted for the Lomonosov, Alpha and Mendeleev Ridges. Crustal basement thicknesses of 10-15 km are predicted under the Laptev Sea which is interpreted as highly thinned continental crust formed at the eastward continuation of Eurasia Basin sea-floor spreading. Thin continental or oceanic crust of thickness 7 km or less is predicted under the North Chukchi Basin and has major implications for understanding the Mesozoic and Cenozoic plate tectonic history of the Siberian and Chukchi Amerasia Basin margins. Restoration of crustal thickness and continent-ocean boundary location from gravity inversion may be used to test and refine plate tectonic reconstructions. Using crustal thickness and continental lithosphere thinning factor maps with superimposed shaded-relief free-air gravity anomaly, we improve the determination of pre-breakup rifted margin conjugacy and sea-floor spreading trajectory within the Arctic and N Atlantic basins. By restoring crustal thickness & continental lithosphere thinning maps of the Eurasia Basin & NE Atlantic to their initial post-breakup configuration we show the geometry and segmentation of the rifted continental margins at their time of breakup, together with the location of highly-stretched failed breakup basins and rifted micro-continents. We interpret gravity inversion crustal thicknesses underneath Morris Jessop Rise & Yermak Plateau as continental crust which provided a barrier to the tectonic and palaeo-oceanic linkage between the Arctic & North Atlantic until the Oligocene. Before this time, we link the seafloor spreading within the Eurasia Basin to that in Baffin Bay.

The non-Gaussian joint probability density function of slope and elevation for a nonlinear gravity wave field. [in ocean surface

NASA Technical Reports Server (NTRS)

Huang, N. E.; Long, S. R.; Bliven, L. F.; Tung, C.-C.

1984-01-01

On the basis of the mapping method developed by Huang et al. (1983), an analytic expression for the non-Gaussian joint probability density function of slope and elevation for nonlinear gravity waves is derived. Various conditional and marginal density functions are also obtained through the joint density function. The analytic results are compared with a series of carefully controlled laboratory observations, and good agreement is noted. Furthermore, the laboratory wind wave field observations indicate that the capillary or capillary-gravity waves may not be the dominant components in determining the total roughness of the wave field. Thus, the analytic results, though derived specifically for the gravity waves, may have more general applications.

The Formation and Fate of Internal Waves in the South China Sea

DTIC Science & Technology

2015-11-05

FOf’miiiiiiM and Fate at Internal Waves In the South •C:hln;~t Sea --- --------· . _.,.. --- -------Author(s) Name{s) (Firsi,MI,La$t), Code, Atfi(iation...Tswen-Yung (David) Tang7 Internal gravity waves , the subsurface analogue of the familiar surface gravity waves that break on beaches, are ubiquitous in...for man-made structures in the ocean4. Generated primarily by the wind and the tides, internal waves can travel thousands of kilometres from their

Mars - Crustal structure inferred from Bouguer gravity anomalies.

NASA Technical Reports Server (NTRS)

Phillips, R. J.; Saunders, R. S.; Conel, J. E.

1973-01-01

Bouguer gravity has been computed for the equatorial region of Mars by differencing free air gravity and the gravity predicted from topographic variations. The free air gravity was generated from an eighth-order set of spherical harmonic coefficients. The gravity from topographic variations was generated by integrating a two-dimensional Green's function over each contour level. The Bouguer gravity indicates crustal inhomogeneities on Mars that are postulated to be variations in crustal thickness. The Tharsis ridge is a region of thick continental type crust. The gravity data, structural patterns, topography, and surface geology of this region lead to the interpretation of the Tharsis topographic high as a broad crustal upwarp possibly associated with local formation of lower-density crustal material and subsequent rise of a thicker crust. The Amazonis region is one of several basins of relatively thin crust, analogous to terrestrial ocean basins. The Libya and Hellas basins, which are probable impact features, are also underlain by thin crust and are possible regions of mantle upwelling.

Detection and Characterization of Deep Water Wave Breaking Using Moderate Incidence Angle Microwave Backscatter from the Sea Surface

DTIC Science & Technology

1990-06-01

interaction and wave breaking. The ocean surface can be modelled as a two-scale or composite surface - 21 - made up of short wind-generated ripples... composite or two-scale rough surface (Barrick and Peake, 1968). For radar wavelengths on the order of a few centimeters, the resonant scatterers are...short wind ripples which ride on top of long gravity waves, and a - 46 - composite model is used to describe the two-scale nature of the sea surface

Waves in Radial Gravity Using Magnetic Fluid

NASA Technical Reports Server (NTRS)

Ohlsen, D. R.; Hart, J. E.; Weidman, P. D.

1999-01-01

Terrestrial laboratory experiments studying various fluid dynamical processes are constrained, by being in an Earth laboratory, to have a gravitational body force which is uniform and unidirectional. Therefore fluid free-surfaces are horizontal and flat. Such free surfaces must have a vertical solid boundary to keep the fluid from spreading horizontally along a gravitational potential surface. In atmospheric, oceanic, or stellar fluid flows that have a horizontal scale of about one-tenth the body radius or larger, sphericity is important in the dynamics. Further, fluids in spherical geometry can cover an entire domain without any sidewall effects, i.e. have truly periodic boundary conditions. We describe spherical body-force laboratory experiments using ferrofluid. Ferrofluids are dilute suspensions of magnetic dipoles, for example magnetite particles of order 10 nm diameter, suspended in a carrier fluid. Ferrofluids are subject to an additional body force in the presence of an applied magnetic field gradient. We use this body force to conduct laboratory experiments in spherical geometry. The present study is a laboratory technique improvement. The apparatus is cylindrically axisymmetric. A cylindrical ceramic magnet is embedded in a smooth, solid, spherical PVC ball. The geopotential field and its gradient, the body force, were made nearly spherical by careful choice of magnet height-to-diameter ratio and magnet size relative to the PVC ball size. Terrestrial gravity is eliminated from the dynamics by immersing the "planet" and its ferrofluid "ocean" in an immiscible silicone oil/freon mixture of the same density. Thus the earth gravity is removed from the dynamics of the ferrofluid/oil interface and the only dynamically active force there is the radial magnetic gravity. The entire apparatus can rotate, and waves are forced on the ferrofluid surface by exterior magnets. The biggest improvement in technique is in the wave visualization. Fluorescing dye is added to the oil/freon mixture and an argon ion laser generates a horizontal light that can be scanned vertically. Viewed from above, the experiment is a black circle with wave deformations surrounded by a light background. A contour of the image intensity at any light sheet position gives the surface of the ferrofluid "ocean" at that "latitude". Radial displacements of the waves as a function of longitude are obtained by subtracting the contour line positions from a no-motion contour at that laser sheet latitude. The experiments are run by traversing the forcing magnet with the laser sheet height fixed and images are frame grabbed to obtain a time-series at one latitude. The experiment is then re-run with another laser-sheet height to generate a full picture of the three-dimensional wave structure in the upper hemisphere of the ball as a function of time. We concentrate here on results of laboratory studies of waves that are important in Earth's atmosphere and especially the ocean. To get oceanic scaling in the laboratory, the experiment must rotate rapidly (4-second rotation period) so that the wave speed is slow compared to the planetary rotation speed as in the ocean. In the Pacific Ocean, eastward propagating Kelvin waves eventually run into the South American coast. Theory predicts that some of the wave energy should scatter into coastal-trapped Kelvin waves that propagate north and south along the coast. Some of this coastal wave energy might then scatter into mid-latitude Rossby waves that propagate back westward. Satellite observations of the Pacific Ocean sea-surface temperature and height seem to show signatures of westward propagating mid-latitude Rossby waves, 5 to 10 years after the 1982-83 El Nino. The observational data is difficult to interpret unambiguously owing to the large range of motions that fill the ocean at shorter timescales. This series of reflections giving eastward, north- ward, and then westward traveling waves is observed cleanly in the laboratory experiments, confirming the theoretical expectations.

That gravity thing

NASA Astrophysics Data System (ADS)

Jewess, Mike

2009-05-01

Your news article "New probe plots Earth's gravity field" (March p11) reports on the European Space Agency's Gravity Field and Steady-State Ocean Circulation Explorer (GOCE) - a satellite that will measure the Earth's gravitational field. It describes the way that g, the acceleration of free fall at the Earth's surface, varies with latitude; this variation is great enough to require adjustment of pendulum clocks between latitudes and also the recalibration of all balances that do not directly compare one mass with a reference mass. The article also notes that the spin of the (effectively fluid) Earth causes it to bulge at the equator, a realization that goes back to Newton's Principia.

Reprocessing the GRACE-derived gravity field time series based on data-driven method for ocean tide alias error mitigation

NASA Astrophysics Data System (ADS)

Liu, Wei; Sneeuw, Nico; Jiang, Weiping

2017-04-01

GRACE mission has contributed greatly to the temporal gravity field monitoring in the past few years. However, ocean tides cause notable alias errors for single-pair spaceborne gravimetry missions like GRACE in two ways. First, undersampling from satellite orbit induces the aliasing of high-frequency tidal signals into the gravity signal. Second, ocean tide models used for de-aliasing in the gravity field retrieval carry errors, which will directly alias into the recovered gravity field. GRACE satellites are in non-repeat orbit, disabling the alias error spectral estimation based on the repeat period. Moreover, the gravity field recovery is conducted in non-strictly monthly interval and has occasional gaps, which result in an unevenly sampled time series. In view of the two aspects above, we investigate the data-driven method to mitigate the ocean tide alias error in a post-processing mode.

Surface Gravity Data Contribution to the Puerto Rico and U.S. Virgin Islands Geoid Model

NASA Astrophysics Data System (ADS)

Li, X.; Gerhards, C.; Holmes, S. A.; Saleh, J.; Shaw, B.

2015-12-01

The Gravity for the Redefinition of the American Vertical Datum (GRAV-D) project provides updated local gravity field information for the XGEOID15 models. In particular, its airborne gravity data in the area of Puerto Rico and U.S. Virgin Islands (PRVI) made substantial improvements (~60%) on the precision of the geoid models at the local GNSS/Leveling bench marks in the target area. Fortunately, PRVI is free of the huge systematic error in the North American Vertical Datum of 1988 (NAVD88). Thus, the airborne contribution was evaluated more realistically. In addition, the airborne data picked up more detailed gravity field information in the medium wavelength band (spherical harmonic degree 200 to 600) that are largely beyond the resolution of the current satellite missions, especially along the nearby ocean trench areas. Under this circumstance (significant airborne contributions in the medium band), local surface gravity data need to be examined more carefully than before during merging with the satellite and airborne information for local geoid improvement, especially considering the well-known systematic problems in the NGS historical gravity holdings (Saleh et al 2013 JoG). Initial tests showed that it is very important to maintain high consistency between the surface data sets and the airborne enhanced reference model. In addition, a new aggregation method (Gerhards 2014, Inverse Problems) will also be tested to optimally combine the local surface data with the reference model. The data cleaning and combining procedures in the target area will be summarized here as reference for future applications.

Langmuir cells and mixing in the upper ocean

NASA Astrophysics Data System (ADS)

Carniel, S.; Sclavo, M.; Kantha, L. H.; Clayson, C. A.

2005-01-01

The presence of surface gravity waves at the ocean surface has two important effects on turbulence in the oceanic mixed layer (ML): the wave breaking and the Langmuir cells (LC). Both these effects act as additional sources of turbulent kinetic energy (TKE) in the oceanic ML, and hence are important to mixing in the upper ocean. The breaking of high wave-number components of the wind wave spectrum provides an intense but sporadic source of turbulence in the upper surface; turbulence thus injected diffuses downward, while decaying rapidly, modifying oceanic near-surface properties which in turn could affect the air-sea transfer of heat and dissolved gases. LC provide another source of additional turbulence in the water column; they are counter-rotating cells inside the ML, with their axes roughly aligned in the direction of the wind (Langmuir I., Science871938119). These structures are usually made evident by the presence of debris and foam in the convergence area of the cells, and are generated by the interaction of the wave-field-induced Stokes drift with the wind-induced shear stress. LC have long been thought to have a substantial influence on mixing in the upper ocean, but the difficulty in their parameterization have made ML modelers consistently ignore them in the past. However, recent Large Eddy Simulations (LES) studies suggest that it is possible to include their effect on mixing by simply adding additional production terms in the turbulence equations, thus enabling even 1D models to incorporate LC-driven turbulence. Since LC also modify the Coriolis terms in the mean momentum equations by the addition of a term involving the Stokes drift, their effect on the velocity structure in the ML is also quite significant and could have a major impact on the drift of objects and spilled oil in the upper ocean. In this paper we examine the effect of surface gravity waves on mixing in the upper ocean, focusing on Langmuir circulations, which is by far the dominant part of the surface wave contribution to mixing. Oceanic ML models incorporating these effects are applied to an observation station in the Northern Adriatic Sea to see what the extent of these effects might be. It is shown that the surface wave effects can indeed be significant; in particular, the modification of the velocity profile due to LC-generated turbulence can be large under certain conditions. However, the surface wave effects on the bulk properties of the ML, such as the associated temperature, while significant, are generally speaking well within the errors introduced by uncertainties in the external forcing of the models. This seems to be the reason why ML models, though pretty much ignoring surface wave effects until recently, have been reasonably successful in depicting the evolution of the mixed layer temperature (MLT) at various timescales.

Geodynamics Branch research report, 1982

NASA Technical Reports Server (NTRS)

Kahn, W. D. (Editor); Cohen, S. C. (Editor)

1983-01-01

The research program of the Geodynamics Branch is summarized. The research activities cover a broad spectrum of geoscience disciplines including space geodesy, geopotential field modeling, tectonophysics, and dynamic oceanography. The NASA programs which are supported by the work described include the Geodynamics and Ocean Programs, the Crustal Dynamics Project, the proposed Ocean Topography Experiment (TOPEX) and Geopotential Research Mission. The individual papers are grouped into chapters on Crustal Movements, Global Earth Dynamics, Gravity Field Model Development, Sea Surface Topography, and Advanced Studies.

Satellite Laser Ranging in the 1990s: Report of the 1994 Belmont Workshop

NASA Technical Reports Server (NTRS)

Degnan, John J. (Editor)

1994-01-01

An international network of 43 stations in 30 countries routinely collects satellite ranging data which is used to study the solid Earth and its interactions with the oceans, atmosphere, and Moon. Data products include centimeter accuracy site positions on a global scale, tectonic plate motions, regional crustal deformation, long wavelength gravity field and geoid, polar motion, and variations in the Earth's spin rate. By calibrating and providing precise orbits for spaceborne microwave altimeters, satellite laser ranging also enables global measurement of sea and ice surface topography, mean sea level, global ocean circulation, and short wavelength gravity fields and marine geoids. It provides tests of general relativity and a means or subnanosecond time transfer. This workshop was convened to define future roles and directions in satellite laser ranging.

A comparison of annual vertical crustal displacements from GPS and Gravity Recovery and Climate Experiment (GRACE) over Europe

NASA Astrophysics Data System (ADS)

van Dam, T.; Wahr, J.; LavalléE, David

2007-03-01

We compare approximately 3 years of GPS height residuals (with respect to the International Terrestrial Reference Frame) with predictions of vertical surface displacements derived from the Gravity Recovery and Climate Experiment (GRACE) gravity fields for stations in Europe. An annual signal fit to the residual monthly heights, corrected for atmospheric pressure and barotropic ocean loading effects, should primarily represent surface displacements due to long-wavelength variations in water storage. A comparison of the annual height signal from GPS and GRACE over Europe indicates that at most sites, the annual signals do not agree in amplitude or phase. We find that unlike the annual signal predicted from GRACE, the annual signal in the GPS heights is not coherent over the region, displaying significant variability from site to site. Confidence in the GRACE data and the unlikely possibility of large-amplitude small-scale features in the load field not captured by the GRACE data leads us to conclude that some of the discrepancy between the GPS and GRACE observations is due to technique errors in the GPS data processing. This is evidenced by the fact that the disagreement between GPS and GRACE is largest at coastal sites, where mismodeling of the semidiurnal ocean tidal loading signal can result in spurious annual signals.

Geologic implications of topographic, gravity, and aeromagnetic data in the northern Yukon-Koyukuk Province and its borderlands, Alaska

NASA Astrophysics Data System (ADS)

Cady, John W.

1989-11-01

The northern Yukon-Koyukuk province is characterized by low elevation and high Bouguer gravity and aeromagnetic anomalies in contrast to the adjacent Brooks Range and Ruby geanticline. Using newly compiled digital topographic, gravity, and aeromagnetic maps, I have divided the province into three geophysical domains. The Koyukuk domain, which is nearly equivalent to the Koyukuk lithotectonic terrane, is a horseshoe-shaped area, open to the south, of low topography, high gravity, and high-amplitude magnetic anomalies caused by an intraoceanic magmatic arc. The Angayucham and Kanuti domains are geophysical subdivisions of the Angayucham lithotectonic terrane that occur along the northern and southeastern margins of the Yukon-Koyukuk province, where oceanic rocks have been thrust over continental rocks of the Brooks Range and Ruby geanticline. Basalt of the Angayucham domain causes strong gravity highs and weak magnetic highs. The Kanuti domain is distinguished from the Angayucham domain by intense magnetic highs caused by cumulus mafic and ultramafic plutonic rocks, abundant ultramafic mantle tectonites, and magnetic syenite and monzonite. Long-wavelength, low-intensity magnetic highs and undulating gravity anomalies indicate an undulating basement surface of varied lithology beneath the Kobuk-Koyukuk and Lower Yukon basins. Modeling of gravity and magnetic anomalies shows that oceanic rocks of the Angayucham and Kanuti domains dip inward beneath the Kobuk-Koyukuk basin. The modeling supports, but does not prove, the hypothesis that the crust of the Kobuk-Koyukuk basin is 32-35 km thick, consisting of a tectonically thickened section of Cretaceous volcanic and sedimentary rocks and older oceanic crust. Plutons of the Brooks Range and the southern Ruby geanticline are nonmagnetic, ilmenite series, S-type granites that cause magnetic lows. Plutons of the northern Ruby geanticline are variable in their magnetic properties and cause both highs and lows. Plutons of both the eastern and western Yukon-Koyukuk province are variable in their magnetic expression but commonly cause magnetic lows in contrast to andesite.

Observations with the ROWS instrument during the Grand Banks calibration/validation experiments

NASA Technical Reports Server (NTRS)

Vandemark, D.; Chapron, B.

1994-01-01

As part of a global program to validate the ocean surface sensors on board ERS-1, a joint experiment on the Grand Banks of Newfoundland was carried out in Nov. 1991. The principal objective was to provide a field validation of ERS-1 Synthetic Aperture Radar (SAR) measurement of ocean surface structure. The NASA-P3 aircraft measurements made during this experiment provide independent measurements of the ocean surface along the validation swath. The Radar Ocean Wave Spectrometer (ROWS) is a radar sensor designed to measure direction of the long wave components using spectral analysis of the tilt induced radar backscatter modulation. This technique greatly differs from SAR and thus, provides a unique set of measurements for use in evaluating SAR performance. Also, an altimeter channel in the ROWS gives simultaneous information on the surface wave height and radar mean square slope parameter. The sets of geophysical parameters (wind speed, significant wave height, directional spectrum) are used to study the SAR's ability to accurately measure ocean gravity waves. The known distortion imposed on the true directional spectrum by the SAR imaging mechanism is discussed in light of the direct comparisons between ERS-1 SAR, airborne Canadian Center for Remote Sensing (CCRS) SAR, and ROWS spectra and the use of the nonlinear ocean SAR transform.

Mean Dynamic Topography of the Arctic Ocean

NASA Technical Reports Server (NTRS)

Farrell, Sinead Louise; Mcadoo, David C.; Laxon, Seymour W.; Zwally, H. Jay; Yi, Donghui; Ridout, Andy; Giles, Katherine

2012-01-01

ICESat and Envisat altimetry data provide measurements of the instantaneous sea surface height (SSH) across the Arctic Ocean, using lead and open water elevation within the sea ice pack. First, these data were used to derive two independent mean sea surface (MSS) models by stacking and averaging along-track SSH profiles gathered between 2003 and 2009. The ICESat and Envisat MSS data were combined to construct the high-resolution ICEn MSS. Second, we estimate the 5.5-year mean dynamic topography (MDT) of the Arctic Ocean by differencing the ICEn MSS with the new GOCO02S geoid model, derived from GRACE and GOCE gravity. Using these satellite-only data we map the major features of Arctic Ocean dynamical height that are consistent with in situ observations, including the topographical highs and lows of the Beaufort and Greenland Gyres, respectively. Smaller-scale MDT structures remain largely unresolved due to uncertainties in the geoid at short wavelengths.

Evidence for infragravity wave-tide resonance in deep oceans.

PubMed

Sugioka, Hiroko; Fukao, Yoshio; Kanazawa, Toshihiko

2010-10-05

Ocean tides are the oscillatory motions of seawater forced by the gravitational attraction of the Moon and Sun with periods of a half to a day and wavelengths of the semi-Pacific to Pacific scale. Ocean infragravity (IG) waves are sea-surface gravity waves with periods of several minutes and wavelengths of several dozen kilometres. Here we report the first evidence of the resonance between these two ubiquitous phenomena, mutually very different in period and wavelength, in deep oceans. The evidence comes from long-term, large-scale observations with arrays of broadband ocean-bottom seismometers located at depths of more than 4,000 m in the Pacific Ocean. This observational evidence is substantiated by a theoretical argument that IG waves and the tide can resonantly couple and that such coupling occurs over unexpectedly wide areas of the Pacific Ocean. Through this resonant coupling, some of ocean tidal energy is transferred in deep oceans to IG wave energy.

Imaging Small-scale Seafloor and Sub-seafloor Tectonic Fabric Using Satellite Altimetry

NASA Astrophysics Data System (ADS)

Sandwell, D. T.; Müller, D.; Matthews, K. J.; Smith, W. H. F.

2017-12-01

Marine gravity anomalies derived from satellite radar altimetry now provide an unprecedented resolution of about 7 km for mapping small-scale seafloor and sub-seafloor tectonic fabric. These gravity maps are improving rapidly because three satellite altimeters are currently collecting data with dense track coverage: (1) CryoSat-2 has routinely collected altimetry data over ice, land, and ocean since July 2010. The satellite has a long 369-day repeat cycle resulting in an average ground track spacing of 3.5 km at the equator. To date it has completed more than 7 geodetic mappings of the ocean surface. (2) The SARAL AltiKa altimeter began a non-repeat orbit phase in July 2016. AltiKa has a new Ka-band instrument with a factor of 2 better range precision than all previous altimeters. (3) Jason-2 was placed in a geodetic orbit starting July 2017. It has lower inclination coverage to provide improved gravity recovery for N-S trending anomalies. These data combined with sparse soundings will provide a dramatic improvement in predicted bathymetry and thus help guide future deep ocean surveys. The most recent global marine gravity anomaly map based on these geodetic mission data with 2-pass retracking for optimal range precision has an accuracy that is 2-4 times better than the maps derived from Geosat and ERS-1. The new data reveal the detailed fabric of fracture zones, previously unmapped, now extinct oceanic microplates in the central Pacific, and fault networks buried beneath thick sediments along continental margins. By combining satellite altimetry with marine magnetic anomalies and seafloor age dates from rock samples we are able to pinpoint the geometry and age of major plate reorganizations, particularly the enigmatic 100 Ma event, which occurred during the Cretaceous Magnetic Superchron.

Complete synthetic seismograms based on a spherical self-gravitating Earth model with an atmosphere-ocean-mantle-core structure

NASA Astrophysics Data System (ADS)

Wang, Rongjiang; Heimann, Sebastian; Zhang, Yong; Wang, Hansheng; Dahm, Torsten

2017-04-01

A hybrid method is proposed to calculate complete synthetic seismograms based on a spherically symmetric and self-gravitating Earth with a multi-layered structure of atmosphere, ocean, mantle, liquid core and solid core. For large wavelengths, a numerical scheme is used to solve the geodynamic boundary-value problem without any approximation on the deformation and gravity coupling. With the decreasing wavelength, the gravity effect on the deformation becomes negligible and the analytical propagator scheme can be used. Many useful approaches are used to overcome the numerical problems that may arise in both analytical and numerical schemes. Some of these approaches have been established in the seismological community and the others are developed for the first time. Based on the stable and efficient hybrid algorithm, an all-in-one code QSSP is implemented to cover the complete spectrum of seismological interests. The performance of the code is demonstrated by various tests including the curvature effect on teleseismic body and surface waves, the appearance of multiple reflected, teleseismic core phases, the gravity effect on long period surface waves and free oscillations, the simulation of near-field displacement seismograms with the static offset, the coupling of tsunami and infrasound waves, and free oscillations of the solid Earth, the atmosphere and the ocean. QSSP is open source software that can be used as a stand-alone FORTRAN code or may be applied in combination with a Python toolbox to calculate and handle Green's function databases for efficient coding of source inversion problems.

Complete synthetic seismograms based on a spherical self-gravitating Earth model with an atmosphere-ocean-mantle-core structure

NASA Astrophysics Data System (ADS)

Wang, Rongjiang; Heimann, Sebastian; Zhang, Yong; Wang, Hansheng; Dahm, Torsten

2017-09-01

A hybrid method is proposed to calculate complete synthetic seismograms based on a spherically symmetric and self-gravitating Earth with a multilayered structure of atmosphere, ocean, mantle, liquid core and solid core. For large wavelengths, a numerical scheme is used to solve the geodynamic boundary-value problem without any approximation on the deformation and gravity coupling. With decreasing wavelength, the gravity effect on the deformation becomes negligible and the analytical propagator scheme can be used. Many useful approaches are used to overcome the numerical problems that may arise in both analytical and numerical schemes. Some of these approaches have been established in the seismological community and the others are developed for the first time. Based on the stable and efficient hybrid algorithm, an all-in-one code QSSP is implemented to cover the complete spectrum of seismological interests. The performance of the code is demonstrated by various tests including the curvature effect on teleseismic body and surface waves, the appearance of multiple reflected, teleseismic core phases, the gravity effect on long period surface waves and free oscillations, the simulation of near-field displacement seismograms with the static offset, the coupling of tsunami and infrasound waves, and free oscillations of the solid Earth, the atmosphere and the ocean. QSSP is open source software that can be used as a stand-alone FORTRAN code or may be applied in combination with a Python toolbox to calculate and handle Green's function databases for efficient coding of source inversion problems.

Assessment of radar altimetry correction slopes for marine gravity recovery: A case study of Jason-1 GM data

NASA Astrophysics Data System (ADS)

Zhang, Shengjun; Li, Jiancheng; Jin, Taoyong; Che, Defu

2018-04-01

Marine gravity anomaly derived from satellite altimetry can be computed using either sea surface height or sea surface slope measurements. Here we consider the slope method and evaluate the errors in the slope of the corrections supplied with the Jason-1 geodetic mission data. The slope corrections are divided into three groups based on whether they are small, comparable, or large with respect to the 1 microradian error in the current sea surface slope models. (1) The small and thus negligible corrections include dry tropospheric correction, inverted barometer correction, solid earth tide and geocentric pole tide. (2) The moderately important corrections include wet tropospheric correction, dual-frequency ionospheric correction and sea state bias. The radiometer measurements are more preferred than model values in the geophysical data records for constraining wet tropospheric effect owing to the highly variable water-vapor structure in atmosphere. The items of dual-frequency ionospheric correction and sea state bias should better not be directly added to range observations for obtaining sea surface slopes since their inherent errors may cause abnormal sea surface slopes and along-track smoothing with uniform distribution weight in certain width is an effective strategy for avoiding introducing extra noises. The slopes calculated from radiometer wet tropospheric corrections, and along-track smoothed dual-frequency ionospheric corrections, sea state bias are generally within ±0.5 microradians and no larger than 1 microradians. (3) Ocean tide has the largest influence on obtaining sea surface slopes while most of ocean tide slopes distribute within ±3 microradians. Larger ocean tide slopes mostly occur over marginal and island-surrounding seas, and extra tidal models with better precision or with extending process (e.g. Got-e) are strongly recommended for updating corrections in geophysical data records.

Large scale mass redistribution and surface displacement from GRACE and SLR

NASA Astrophysics Data System (ADS)

Cheng, M.; Ries, J. C.; Tapley, B. D.

2012-12-01

Mass transport between the atmosphere, ocean and solid earth results in the temporal variations in the Earth gravity field and loading induced deformation of the Earth. Recent space-borne observations, such as GRACE mission, are providing extremely high precision temporal variations of gravity field. The results from 10-yr GRACE data has shown a significant annual variations of large scale vertical and horizontal displacements occurring over the Amazon, Himalayan region and South Asia, African, and Russian with a few mm amplitude. Improving understanding from monitoring and modeling of the large scale mass redistribution and the Earth's response are a critical for all studies in the geosciences, in particular for determination of Terrestrial Reference System (TRS), including geocenter motion. This paper will report results for the observed seasonal variations in the 3-dimentional surface displacements of SLR and GPS tracking stations and compare with the prediction from time series of GRACE monthly gravity solution.

Anisotropic models of the upper mantle

NASA Technical Reports Server (NTRS)

Regan, J.; Anderson, D. L.

1983-01-01

Long period Rayleigh wave and Love wave dispersion data, particularly for oceanic areas, were not simultaneously satisfied by an isotropic structure. Available phase and group velocity data are inverted by a procedure which includes the effects of transverse anisotropy, an elastic dispersion, sphericity, and gravity. The resulting models, for the average Earth, average ocean and oceanic regions divided according to the age of the ocean floor, are quite different from previous results which ignore the above effects. The models show a low velocity zone with age dependent anisotropy and velocities higher than derived in previous surface wave studies. The correspondence between the anisotropy variation with age and a physical model based on flow aligned olivine is suggested.

Satellite Laser Ranging operations

NASA Technical Reports Server (NTRS)

Pearlman, Michael R.

1994-01-01

Satellite Laser Ranging (SLR) is currently providing precision orbit determination for measurements of: 1) Ocean surface topography from satellite borne radar altimetry, 2) Spatial and temporal variations of the gravity field, 3) Earth and ocean tides, 4) Plate tectonic and regional deformation, 5) Post-glacial uplift and subsidence, 6) Variations in the Earth's center-of-mass, and 7) Variations in Earth rotation. SLR also supports specialized programs in time transfer and classical geodetic positioning, and will soon provide precision ranging to support experiments in relativity.

From Germany to Antarctica: Airborne geodesy and geophysics and the utilization of the research aircraft HALO (Invited)

NASA Astrophysics Data System (ADS)

Scheinert, M.; Barthelmes, F.; Foerste, C.; Heyde, I.

2013-12-01

The geoid as an equipotential surface of the gravity potential plays a crucial role for the realiziation of the Global Geodetic Observation System (GGOS) of IAG (International Association of Geodesy). It is the major reference surface for physical height systems. The gravity potential is needed to precisely predict the orbits of artificial satellites of the earth. A precise static solution enters analyses of temporal changes of the gravity field due to mass transport processes between the different subsystems of the earth. However, also in neighbouring disciplines the geoid is applied. In oceanography, for example, the geoid serves as a reference surface for the determination of the mean sea-surface topography (MSST). In glaciology, it enters analyses of the thickness of ice bodies floating in polar waters, based on freeboard heights and the equilibrium supposition. To come up with high resolution global gravity field models, satellite observations - preferably of the dedicated satellite gravity missions - have to be combined with surface gravity data. Although the majority of the continental surface is captured by ground-based or near-surface gravity measurements - and gravity over the oceans is determined by satellite altimetry - still large gaps in surface gravity data exist. In this respect it is the Antarctic continent which suffers large data gaps, not only in surface gravity but also due to the polar gap of GOCE satellite gravimetry. Chairing the IAG Subcommission 2.4f 'Gravity and Geoid in Antarctica' (AntGG) the author will discuss the current status of gravity surveys in Antarctica. Especially airborne gravimetry has been and is being widely applied as the only reasonable method to survey large areas in this vast and hostile environment. As a novel application the German research aircraft HALO was utilized for a geodetic-geophysical flight mission. Measurements were realized to acquire data of the gravity and magnetic fields, of GNSS remote sensing and of laser altimetry over Italy and adjacent (Tyrrhenian, Adriatic and Ionian) seas. This so-called GEOHALO flight mission was carried out in the time period from June 2 to 12, 2012. The flights comprised seven parallel profiles directing from north-west to south-east, in a height of about 3,500 m, with a length of about 1,000 km each and a line spacing of about 40 km. These long profiles were complemented by four crossing profiles and a profile at an altitude of approx. 10 km along the same track as the center long profile. Special focus will be given to the results of airborne gravimetry and laser altimetry to further investigate the gravity field and the sea-surface topography in the Mediterranean. Furthermore, the status of HALO and future plans to utilize HALO for an Antarctic flight mission will be discussed. Applications of airborne gravimetry to investigate geodetic problems in Antarctica shall be shortly discussed, together with an outlook of AntGG.

Ocean wave-radar modulation transfer functions from the West Coast experiment

NASA Technical Reports Server (NTRS)

Wright, J. W.; Plant, W. J.; Keller, W. C.; Jones, W. L.

1980-01-01

Short gravity-capillary waves, the equilibrium, or the steady state excitations of the ocean surface are modulated by longer ocean waves. These short waves are the predominant microwave scatterers on the ocean surface under many viewing conditions so that the modulation is readily measured with CW Doppler radar used as a two-scale wave probe. Modulation transfer functions (the ratio of the cross spectrum of the line-of-sight orbital speed and backscattered microwave power to the autospectrum of the line-of-sight orbital speed) were measured at 9.375 and 1.5 GHz (Bragg wavelengths of 2.3 and 13 cm) for winds up to 10 m/s and ocean wave periods from 2-18 s. The measurements were compared with the relaxation-time model; the principal result is that a source of modulation other than straining by the horizontal component of orbital speed, possibly the wave-induced airflow, is responsible for most of the modulation by waves of typical ocean wave period (10 s). The modulations are large; for unit coherence, spectra of radar images of deep-water waves should be proportional to the quotient of the slope spectra of the ocean waves by the ocean wave frequency.

Determining OCT structure and COB Location of the Omani Gulf of Aden Continental Margin from Gravity Inversion, Residual Depth Anomaly and Subsidence Analysis.

NASA Astrophysics Data System (ADS)

Cowie, Leanne; Kusznir, Nick; Leroy, Sylvie; Manatshal, Gianreto

2013-04-01

Knowledge and understanding of the ocean-continent transition (OCT) structure and continent-ocean boundary (COB) location, the distribution of thinned continental crust and lithosphere, its distal extent and the start of unequivocal oceanic crust are of critical importance in evaluating rifted continental margin formation and evolution. In order to determine the OCT structure and COB location for the eastern Gulf of Aden, along the Oman margin, we use a combination of gravity inversion, subsidence analysis and residual depth anomaly (RDA) analysis. Gravity inversion has been used to determine Moho depth, crustal basement thickness and continental lithosphere thinning; subsidence analysis has been used to determine the distribution of continental lithosphere thinning; and RDAs have been used to investigate the OCT bathymetric anomalies with respect to expected oceanic bathymetries at rifted margins. The gravity inversion method, which is carried out in the 3D spectral domain, incorporates a lithosphere thermal gravity anomaly and includes a correction for volcanic addition due to decompression melting. Reference Moho depths used in the gravity inversion have been calibrated against seismic refraction Moho depths. RDAs have been calculated by comparing observed and age predicted oceanic bathymetries, using the thermal plate model predictions from Crosby and McKenzie (2009). RDAs have been computed along profiles and have been corrected for sediment loading using flexural back-stripping and decompaction. In addition, gravity inversion crustal basement thicknesses together with Airy isostasy have been used to predict a synthetic RDA. The RDA results show a change in RDA signature and may be used to estimate the distal extent of thinned continental crust and where oceanic crust begins. Continental lithosphere thinning has been determined using flexural back-stripping and subsidence analysis assuming the classical rift model of McKenzie (1978) with a correction for volcanic addition due to decompression melting based on White & McKenzie (1989). Gravity inversion and the "synthetic" gravity derived RDA both show generally normal thickness oceanic crust, with some localised thin oceanic crust. Continental lithosphere thinning factors determined from gravity inversion and subsidence analysis are in good agreement and have been used to constrain COB location along the profile lines. These techniques show that the OCT in the eastern Gulf of Aden, is relatively narrow, with the distance between the COB and the margin hinge measuring less than 100km.

Arctic Ocean Gravity Field Derived From ERS-1 Satellite Altimetry.

PubMed

Laxon, S; McAdoo, D

1994-07-29

The derivation of a marine gravity field from satellite altimetry over permanently ice-covered regions of the Arctic Ocean provides much new geophysical information about the structure and development of the Arctic sea floor. The Arctic Ocean, because of its remote location and perpetual ice cover, remains from a tectonic point of view the most poorly understood ocean basin on Earth. A gravity field has been derived with data from the ERS-1 radar altimeter, including permanently ice-covered regions. The gravity field described here clearly delineates sections of the Arctic Basin margin along with the tips of the Lomonosov and Arctic mid-ocean ridges. Several important tectonic features of the Amerasia Basin are clearly expressed in this gravity field. These include the Mendeleev Ridge; the Northwind Ridge; details of the Chukchi Borderland; and a north-south trending, linear feature in the middle of the Canada Basin that apparently represents an extinct spreading center that "died" in the Mesozoic. Some tectonic models of the Canada Basin have proposed such a failed spreading center, but its actual existence and location were heretofore unknown.

Application of Satellite Gravimetry for Water Resource Vulnerability Assessment

NASA Technical Reports Server (NTRS)

Rodell, Matthew

2012-01-01

The force of Earth's gravity field varies in proportion to the amount of mass near the surface. Spatial and temporal variations in the gravity field can be measured via their effects on the orbits of satellites. The Gravity Recovery and Climate Experiment (GRACE) is the first satellite mission dedicated to monitoring temporal variations in the gravity field. The monthly gravity anomaly maps that have been delivered by GRACE since 2002 are being used to infer changes in terrestrial water storage (the sum of groundwater, soil moisture, surface waters, and snow and ice), which are the primary source of gravity variability on monthly to decadal timescales after atmospheric and oceanic circulation effects have been removed. Other remote sensing techniques are unable to detect water below the first few centimeters of the land surface. Conventional ground based techniques can be used to monitor terrestrial water storage, but groundwater, soil moisture, and snow observation networks are sparse in most of the world, and the countries that do collect such data rarely are willing to share them. Thus GRACE is unique in its ability to provide global data on variations in the availability of fresh water, which is both vital to life on land and vulnerable to climate variability and mismanagement. This chapter describes the unique and challenging aspects of GRACE terrestrial water storage data, examples of how the data have been used for research and applications related to fresh water vulnerability and change, and prospects for continued contributions of satellite gravimetry to water resources science and policy.

Intercomparison of Satellite Derived Gravity Time Series with Inferred Gravity Time Series from TOPEX/POSEIDON Sea Surface Heights and Climatological Model Output

NASA Technical Reports Server (NTRS)

Cox, C.; Au, A.; Klosko, S.; Chao, B.; Smith, David E. (Technical Monitor)

2001-01-01

The upcoming GRACE mission promises to open a window on details of the global mass budget that will have remarkable clarity, but it will not directly answer the question of what the state of the Earth's mass budget is over the critical last quarter of the 20th century. To address that problem we must draw upon existing technologies such as SLR, DORIS, and GPS, and climate modeling runs in order to improve our understanding. Analysis of long-period geopotential changes based on SLR and DORIS tracking has shown that addition of post 1996 satellite tracking data has a significant impact on the recovered zonal rates and long-period tides. Interannual effects such as those causing the post 1996 anomalies must be better characterized before refined estimates of the decadal period changes in the geopotential can be derived from the historical database of satellite tracking. A possible cause of this anomaly is variations in ocean mass distribution, perhaps associated with the recent large El Nino/La Nina. In this study, a low-degree spherical harmonic gravity time series derived from satellite tracking is compared with a TOPEX/POSEIDON-derived sea surface height time series. Corrections for atmospheric mass effects, continental hydrology, snowfall accumulation, and ocean steric model predictions will be considered.

Goce and Its Role in Combined Global High Resolution Gravity Field Determination

NASA Astrophysics Data System (ADS)

Fecher, T.; Pail, R.; Gruber, T.

2013-12-01

Combined high-resolution gravity field models serve as a mandatory basis to describe static and dynamic processes in system Earth. Ocean dynamics can be modeled referring to a high-accurate geoid as reference surface, solid earth processes are initiated by the gravity field. Also geodetic disciplines such as height system determination depend on high-precise gravity field information. To fulfill the various requirements concerning resolution and accuracy, any kind of gravity field information, that means satellite as well as terrestrial and altimetric gravity field observations have to be included in one combination process. A key role is here reserved for GOCE observations, which contribute with its optimal signal content in the long to medium wavelength part and enable a more accurate gravity field determination than ever before especially in areas, where no high-accurate terrestrial gravity field observations are available, such as South America, Asia or Africa. For our contribution we prepare a combined high-resolution gravity field model up to d/o 720 based on full normal equation including recent GOCE, GRACE and terrestrial / altimetric data. For all data sets, normal equations are set up separately, relative weighted to each other in the combination step and solved. This procedure is computationally challenging and can only be performed using super computers. We put special emphasis on the combination process, for which we modified especially our procedure to include GOCE data optimally in the combination. Furthermore we modified our terrestrial/altimetric data sets, what should result in an improved outcome. With our model, in which we included the newest GOCE TIM4 gradiometry results, we can show how GOCE contributes to a combined gravity field solution especially in areas of poor terrestrial data coverage. The model is validated by independent GPS leveling data in selected regions as well as computation of the mean dynamic topography over the oceans. Further, we analyze the statistical error estimates derived from full covariance propagation and compare them with the absolute validation with independent data sets.

Coseismic Compression/Dilatation and Viscoelastic Uplift/Subsidence Following the 2012 Indian Ocean Earthquakes Quantified from Satellite Gravity Observations

NASA Technical Reports Server (NTRS)

Han, Shin-Chan; Sauber, Jeanne; Pollitz, Fred

2015-01-01

The 2012 Indian Ocean earthquake sequence (M(sub w) 8.6, 8.2) is a rare example of great strike slip earthquakes in an intra-oceanic setting. With over a decade of GRACE data, we were able to measure and model the unanticipated large co-, and post-seismic gravity changes of these events. Using the approach of normal mode decomposition and spatial localization, we computed the gravity changes corresponding to five moment tensor components. Our analysis revealed that the gravity changes are produced predominantly by coseismic compression and dilatation within the oceanic crust and upper mantle and by post-seismic vertical motion. Our results suggest that the post-seismic positive gravity and the post-seismic uplift measured with GPS within the coseismic compressional quadrant are best fit by ongoing uplift associated with viscoelastic mantle relaxation. Our study demonstrates that the GRACE data are suitable for analyzing strike-slip earthquakes as small as M(sub w) 8.2 with the noise characteristics of this region.

Frequency content of sea surface height variability from internal gravity waves to mesoscale eddies

NASA Astrophysics Data System (ADS)

Savage, Anna C.; Arbic, Brian K.; Richman, James G.; Shriver, Jay F.; Alford, Matthew H.; Buijsman, Maarten C.; Thomas Farrar, J.; Sharma, Hari; Voet, Gunnar; Wallcraft, Alan J.; Zamudio, Luis

2017-03-01

High horizontal-resolution (1/12.5° and 1/25°) 41-layer global simulations of the HYbrid Coordinate Ocean Model (HYCOM), forced by both atmospheric fields and the astronomical tidal potential, are used to construct global maps of sea surface height (SSH) variability. The HYCOM output is separated into steric and nonsteric and into subtidal, diurnal, semidiurnal, and supertidal frequency bands. The model SSH output is compared to two data sets that offer some geographical coverage and that also cover a wide range of frequencies—a set of 351 tide gauges that measure full SSH and a set of 14 in situ vertical profilers from which steric SSH can be calculated. Three of the global maps are of interest in planning for the upcoming Surface Water and Ocean Topography (SWOT) two-dimensional swath altimeter mission: (1) maps of the total and (2) nonstationary internal tidal signal (the latter calculated after removing the stationary internal tidal signal via harmonic analysis), with an average variance of 1.05 and 0.43 cm2, respectively, for the semidiurnal band, and (3) a map of the steric supertidal contributions, which are dominated by the internal gravity wave continuum, with an average variance of 0.15 cm2. Stationary internal tides (which are predictable), nonstationary internal tides (which will be harder to predict), and nontidal internal gravity waves (which will be very difficult to predict) may all be important sources of high-frequency "noise" that could mask lower frequency phenomena in SSH measurements made by the SWOT mission.

New observations of Yanai waves and equatorial inertia-gravity waves in the Pacific Ocean

NASA Astrophysics Data System (ADS)

Farrar, J. T.; Durland, T.

2011-12-01

In the 1970's and 1980's, there was a great deal of research activity on near-equatorial variability at periods of days to weeks associated with oceanic equatorial inertia-gravity waves and Yanai waves. At that time, the measurements available for studying these waves were much more limited than today: most of the available observations were from island tide gauges and a handful of short mooring records. We use more than a decade of the extensive modern data record from the TAO/TRITON mooring array in the Pacific Ocean to re-examine the internal-wave climate in the equatorial Pacific, with a focus on interpretation of the zonal-wavenumber/frequency spectrum of surface dynamic height relative to 500-m depth. Many equatorial-wave meridional modes can be identified, for both the first and second baroclinic mode. We also estimated zonal-wavenumber/frequency spectra for the zonal and meridional wind stress components. The location and extent of spectral peaks in dynamic height is readily rationalized using basic, linear theory of forced equatorial waves and the observed wind stress spectrum.

Determination of Mean Dynamic Topography (MDT) to Bridge Geoid and Mean Sea Surface Height (SSH) with a New Elliptic Equation

NASA Astrophysics Data System (ADS)

Chu, P. C.

2016-12-01

Mean dynamic topography (MDT, η) bridges the geoid and the mean sea surface (from satellite altimetry) and constrains large scale surface geostrophic circulations. It can be estimated from either satellite or underwater ocean temperature (T) and salinity (S) data. Satellite altimeter measures sea surface height (SSH) with high precision and unique resolution above a reference ellipsoid (not geoid). Two Gravity Recovery and Climate Experiment (GRACE) satellites launched in 2002, provide data to compute the marine geoid [called the GRACE Gravity Model (GGM)] (see website: http://www.csr.utexas.edu/grace/). The MDT is the difference of altimetry-derived mean SSH and the mean marine geoid (using GGM or pre-GRACE gravity model such as EGM96). A major difficulty arises that the spatial variations in mean SSH and marine geoid are approximately two orders of magnitude larger than the spatial variations in η.The second approach (using T, Sdata) is based on geostrophic balance, which is at the minimum energy state in the linear Boussinesq primitive equations with conservation of potential vorticity. In this paper, a new elliptic equation, -[∂x(gh/f2)∂xη+∂y(gh/f2)∂yη]+η = (g/f2)(∂C/∂x-∂B/∂y)is derived to determine MDT with H the water depth, g the gravitational acceleration, and coefficients (B, C) depend on 3D mean temperature (T) and salinity (S) data. Numerical approach transforms the elliptic equation into a set of well-posed linear algebraic equations of η at grid points. The solution for the North Atlantic Ocean (100oW-6oW, 7oN-72oN) on 1oX1ogrids with the coefficients (B, C) calculated from the three-dimensional (T, S) data of the NOAA National Centers for Environmental Information (NCEI) World Ocean Atlas 2013 version 2 (http://www.nodc.noaa.gov/OC5/woa13/woa13data.html) and H from the NOAA ETOPO5 (https://www.ngdc.noaa.gov/mgg/fliers/93mgg01.html), compares well with the difference (also considered as the MDT) between the time-averaged SSH and the geoid from the NASA/JPL (http://gracetellus.jpl.nasa.gov/data/dot/). Further application of this elliptic equation method on the high-precision altimetry measurements of SSH such as the Surface Water and Ocean Topography (SWOT) is also presented.

Assessing the performance of formulations for nonlinear feedback of surface gravity waves on ocean currents over coastal waters

NASA Astrophysics Data System (ADS)

Wang, Pengcheng; Sheng, Jinyu; Hannah, Charles

2017-08-01

This study presents applications of a two-way coupled wave-circulation modelling system over coastal waters, with a special emphasis of performance assessments of two different methods for nonlinear feedback of ocean surface gravity waves on three-dimensional (3D) ocean currents. These two methods are the vortex force (VF) formulation suggested by Bennis et al. (2011) and the latest version of radiation stress (RS) formulation suggested by Mellor (2015). The coupled modelling system is first applied to two idealized test cases of surf-zone scales to validate implementations of these two methods in the coupled wave-circulation system. Model results show that the latest version of RS has difficulties in producing the undertow over the surf zone. The coupled system is then applied to Lunenburg Bay (LB) of Nova Scotia during Hurricane Juan in 2003. The coupled system using both the VF and RS formulations generates much stronger and more realistic 3D circulation in the Bay during Hurricane Juan than the circulation-only model, demonstrating the importance of surface wave forces to the 3D ocean circulation over coastal waters. However, the RS formulation generates some weak unphysical currents outside the wave breaking zone due to a less reasonable representation for the vertical distribution of the RS gradients over a slopping bottom. These weak unphysical currents are significantly magnified in a two-way coupled system when interacting with large surface waves, degrading the model performance in simulating currents at one observation site. Our results demonstrate that the VF formulation with an appropriate parameterization of wave breaking effects is able to produce reasonable results for applications over coastal waters during extreme weather events. The RS formulation requires a complex wave theory rather than the linear wave theory for the approximation of a vertical RS term to improve its performance under both breaking and non-breaking wave conditions.

Inversion of gravity and bathymetry in oceanic regions for long-wavelength variations in upper mantle temperature and composition

NASA Technical Reports Server (NTRS)

Solomon, Sean C.; Jordan, Thomas H.

1993-01-01

Long-wavelength variations in geoid height, bathymetry, and SS-S travel times are all relatable to lateral variations in the characteristic temperature and bulk composition of the upper mantle. The temperature and composition are in turn relatable to mantle convection and the degree of melt extraction from the upper mantle residuum. Thus the combined inversion of the geoid or gravity field, residual bathymetry, and seismic velocity information offers the promise of resolving fundamental aspects of the pattern of mantle dynamics. The use of differential body wave travel times as a measure of seismic velocity information, in particular, permits resolution of lateral variations at scales not resolvable by conventional global or regional-scale seismic tomography with long-period surface waves. These intermediate scale lengths, well resolved in global gravity field models, are crucial for understanding the details of any chemical or physical layering in the mantle and of the characteristics of so-called 'small-scale' convection beneath oceanic lithosphere. In 1991 a three-year project to the NASA Geophysics Program was proposed to carry out a systematic inversion of long-wavelength geoid anomalies, residual bathymetric anomalies, and differential SS-S travel time delays for the lateral variation in characteristic temperature and bulk composition of the oceanic upper mantle. The project was funded as a three-year award, beginning on 1 Jan. 1992.

Estimating Gravity Biases with Wavelets in Support of a 1-cm Accurate Geoid Model

NASA Astrophysics Data System (ADS)

Ahlgren, K.; Li, X.

2017-12-01

Systematic errors that reside in surface gravity datasets are one of the major hurdles in constructing a high-accuracy geoid model at high resolutions. The National Oceanic and Atmospheric Administration's (NOAA) National Geodetic Survey (NGS) has an extensive historical surface gravity dataset consisting of approximately 10 million gravity points that are known to have systematic biases at the mGal level (Saleh et al. 2013). As most relevant metadata is absent, estimating and removing these errors to be consistent with a global geopotential model and airborne data in the corresponding wavelength is quite a difficult endeavor. However, this is crucial to support a 1-cm accurate geoid model for the United States. With recently available independent gravity information from GRACE/GOCE and airborne gravity from the NGS Gravity for the Redefinition of the American Vertical Datum (GRAV-D) project, several different methods of bias estimation are investigated which utilize radial basis functions and wavelet decomposition. We estimate a surface gravity value by incorporating a satellite gravity model, airborne gravity data, and forward-modeled topography at wavelet levels according to each dataset's spatial wavelength. Considering the estimated gravity values over an entire gravity survey, an estimate of the bias and/or correction for the entire survey can be found and applied. In order to assess the accuracy of each bias estimation method, two techniques are used. First, each bias estimation method is used to predict the bias for two high-quality (unbiased and high accuracy) geoid slope validation surveys (GSVS) (Smith et al. 2013 & Wang et al. 2017). Since these surveys are unbiased, the various bias estimation methods should reflect that and provide an absolute accuracy metric for each of the bias estimation methods. Secondly, the corrected gravity datasets from each of the bias estimation methods are used to build a geoid model. The accuracy of each geoid model provides an additional metric to assess the performance of each bias estimation method. The geoid model accuracies are assessed using the two GSVS lines and GPS-leveling data across the United States.

Three-Dimensional Acoustic Propagation Through Shallow Water Internal, Surface Gravity and Bottom Sediment Waves

DTIC Science & Technology

2011-09-01

energy never ends. I am also very pleased to have Dr. William M. Carey, Dr. Henrik Schmidt, Dr. Glen G. Gawarkiewicz and Dr. Pierre Lermusiaux on my...Internal Waves for Multi- Megameter Acoustic Propagation in the Ocean, J. Acoust. Soc. Amer., Vol. 100, P. 3607-3620, 1996. [6] J.R. Apel , M. Badiey

South China Sea crustal thickness and lithosphere thinning from satellite gravity inversion incorporating a lithospheric thermal gravity anomaly correction

NASA Astrophysics Data System (ADS)

Kusznir, Nick; Gozzard, Simon; Alvey, Andy

2016-04-01

The distribution of ocean crust and lithosphere within the South China Sea (SCS) are controversial. Sea-floor spreading re-orientation and ridge jumps during the Oligocene-Miocene formation of the South China Sea led to the present complex distribution of oceanic crust, thinned continental crust, micro-continents and volcanic ridges. We determine Moho depth, crustal thickness and continental lithosphere thinning (1- 1/beta) for the South China Sea using a gravity inversion method which incorporates a lithosphere thermal gravity anomaly correction (Chappell & Kusznir, 2008). The gravity inversion method provides a prediction of ocean-continent transition structure and continent-ocean boundary location which is independent of ocean isochron information. A correction is required for the lithosphere thermal gravity anomaly in order to determine Moho depth accurately from gravity inversion; the elevated lithosphere geotherm of the young oceanic and rifted continental margin lithosphere of the South China Sea produces a large lithosphere thermal gravity anomaly which in places exceeds -150 mGal. The gravity anomaly inversion is carried out in the 3D spectral domain (using Parker 1972) to determine 3D Moho geometry and invokes Smith's uniqueness theorem. The gravity anomaly contribution from sediments assumes a compaction controlled sediment density increase with depth. The gravity inversion includes a parameterization of the decompression melting model of White & McKenzie (1999) to predict volcanic addition generated during continental breakup lithosphere thinning and seafloor spreading. Public domain free air gravity anomaly, bathymetry and sediment thickness data are used in this gravity inversion. Using crustal thickness and continental lithosphere thinning factor maps with superimposed shaded-relief free-air gravity anomaly, we improve the determination of pre-breakup rifted margin conjugacy, rift orientation and sea-floor spreading trajectory. SCS conjugate margins are highly asymmetric and have several striking features such as the Macclesfield Bank, Xisha Trough, Reed Bank and Dangerous Grounds. Thin continental crust is predicted extending westwards from thin oceanic crust north of Macclesfield Bank into the Quiondongnan (QDN) basin and is interpreted as being generated ahead of westward propagating sea-floor spreading most in the Oligocene. Further south, highly thinned continental crust or possibly serpentinised exhumed mantle is predicted in the Phu Khanh Basin. Ahead of the failed propagating tip of seafloor spreading, offshore southern Vietnam, thinned continental crust is predicted for the Cuu Long and Nam Con Son Basins. Crustal thicknesses from gravity inversion confirms that the southern margin of the SCS consists of fragmented blocks of thinned continental crust separated by thinner regions of continental crust that have undergone higher degrees of stretching and thinning. The Reed Bank is predicted to have a crustal thickness of 20 to 25km, similar to that of Macclesfield Bank. The Dangerous Grounds, west of the Reed Bank, are also predicted to consist of continental crust. This region has been thinned to a higher degree than the Reed Bank, with continental crustal thickness ranging between 10 and 20km thick.

Gravity is the Key Experiment to Address the Habitability of the Ocean in Jupiter's Moon Europa

NASA Astrophysics Data System (ADS)

Sessa, A. M.; Dombard, A. J.

2013-12-01

Life requires three constituents: a liquid solvent (i.e., water), a chemical system that can form large molecules to record genetic information (e.g., carbon based) as well as chemical nutrients (e.g., nitrogen, phosphorous), and a chemical disequilibrium system that can provide metabolic energy. While it is believed that there is a saline water layer located between the rock and ice layers in Jupiter's moon Europa, which would satisfy the first requirement, it is unknown if the other conditions are currently met. The likelihood that Europa is a haven for life in our Solar System skyrockets, however, if there is currently active volcanism at the rock-water interface, much the same that volcanic processes enable the chemosynthetic life that forms the basis of deep sea-vent communities at the bottom of Earth's oceans. Exploring the volcanic activity on this interface is challenging, as direct observation via a submersible or high-resolution indirect observations via a dense global seismic network on the surface is at present technically (and fiscally!) untenable. Thus, gravity studies are the best way to explore currently the structure of this all-important interface. Though mostly a silicate body with only a relatively thin (~100 km) layer of water, Europa is different from the terrestrial planets in that this rock-water interface, and not the surface, represents the largest density contrast across the moon's near-surface layers, and thus topography on this interface could conceivably dominate the gravity. Here, we calculate the potential anomalies that arise from topography on the surface, the water-ice interface (at 20 km depth), and the rock-water interface, finding that the latter dominates the free-air gravity at the longest wavelengths (spherical harmonic degrees < 10) and the Bouguer gravity at intermediate wavelengths (degrees ~10-50), and only for the shortest wavelengths (degrees > 50) does the water-ice interface (and presumably mass-density anomalies within the ice shell) dominate the Bouguer gravity. Thus, gravity can be used to explore this interface. To test whether active volcanism can be detected, we scale gravity models for the terrestrial planets down to a body the size of Europa's silicate core and with a density contrast consistent with a rock-water interface. Here, Venus and Earth serve as proxies for volcanically active bodies, while the Moon and Mars are proxies for inactive bodies. Additionally, we create gravity from synthetic topography on the base of the ice shell. Maps of the Bouguer-gravity and geoid anomalies reveal that active volcanism is characterized by small amplitudes (a few mGal and a few meters). Large-scale topography on the base of the ice shell adds larger geoid anomalies (tens of meters) but still small gravity anomalies. The absence of volcanic activity on the rock-water interface is likely characterized by larger anomalies (tens of mGal and tens of meters), plausibly because the cooler thermal structure permits the rocky lithosphere to support larger mass-density anomalies. Thus, study of the gravity may illuminate the habitability of Europa, and gravity and topography experiments on any future mission (e.g., the Europa Clipper) should be given the highest scientific priority.

Global mean dynamic topography based on GOCE data and Wiener filters

NASA Astrophysics Data System (ADS)

Gilardoni, Maddalena; Reguzzoni, Mirko; Albertella, Alberta

2015-04-01

A mean dynamic ocean topography (MDT) has been computed by using a GOCE-only gravity model and a given mean sea surface (MSS) obtained from satellite altimetry. Since the used gravity model, i.e. the fifth release of the time-wise solution covering the full mission lifetime, is truncated at a maximum harmonic degree of 280, the obtained MDT has to be consistently filtered. This has been done globally by using the spherical harmonic representation and following a Wiener minimization principle. This global filtering approach is convenient from the computational point of view but requires to have MDT values all over the Earth surface and therefore to fill the continents with fictitious data. The main improvements with respect to the already presented results are in the MDT filling procedure (to guarantee that the global signal has the same covariance of the one over the oceans), in the error modelling of the input MSS and in the error estimation of the filtered MDT and of the corresponding geostrophic velocities. The impact of GOCE data in the ocean circulation global modelling has been assessed by comparing the pattern of the obtained geostrophic currents with those computed by using EGM2008. Comparisons with independent circulation data based on drifters and other MDT models have been also performed with the aim of evaluating the accuracy of the obtained results.

Recent results on modelling the spatial and temporal structure of the Earth's gravity field.

PubMed

Moore, P; Zhang, Q; Alothman, A

2006-04-15

The Earth's gravity field plays a central role in sea-level change. In the simplest application a precise gravity field will enable oceanographers to capitalize fully on the altimetric datasets collected over the past decade or more by providing a geoid from which absolute sea-level topography can be recovered. However, the concept of a static gravity field is now redundant as we can observe temporal variability in the geoid due to mass redistribution in or on the total Earth system. Temporal variability, associated with interactions between the land, oceans and atmosphere, can be investigated through mass redistributions with, for example, flow of water from the land being balanced by an increase in ocean mass. Furthermore, as ocean transport is an important contributor to the mass redistribution the time varying gravity field can also be used to validate Global Ocean Circulation models. This paper will review the recent history of static and temporal gravity field recovery, from the 1980s to the present day. In particular, mention will be made of the role of satellite laser ranging and other space tracking techniques, satellite altimetry and in situ gravity which formed the basis of gravity field determination until the last few years. With the launch of Challenging Microsatellite Payload and Gravity and Circulation Experiment (GRACE) our knowledge of the spatial distribution of the Earth's gravity field is taking a leap forward. Furthermore, GRACE is now providing insight into temporal variability through 'monthly' gravity field solutions. Prior to this data we relied on satellite tracking, Global Positioning System and geophysical models to give us insight into the temporal variability. We will consider results from these methodologies and compare them to preliminary results from the GRACE mission.

Numerical Investigations of Wave-Induced Mixing in Upper Ocean Layer

NASA Astrophysics Data System (ADS)

Guan, Changlong

2017-04-01

The upper ocean layer is playing an important role in ocean-atmosphere interaction. The typical characteristics depicting the upper ocean layer are the sea surface temperature (SST) and the mixed layer depth (MLD). So far, the existing ocean models tend to over-estimate SST and to under-estimate MLD, due to the inadequate mixing in the mixing layer, which is owing to that several processes related mixing in physics are ignored in these ocean models. The mixing induced by surface gravity wave is expected to be able to enhance the mixing in the upper ocean layer, and therefore the over-estimation of SST and the under-estimate of MLD could be improved by including wave-induced mixing. The wave-induced mixing could be accomplished by the physical mechanisms, such as wave breaking (WB), wave-induced Reynolds stress (WR), and wave-turbulence interaction (WT). The General Ocean Turbulence Model (GOTM) is employed to investigate the effects of the three mechanisms concerning wave-induced mixing. The numerical investigation is carried out for three turbulence closure schemes, say, k-epsilon, k-omega and Mellor-Yamada (1982), with the observational data from OSC Papa station and wave data from ECMWF. The mixing enhancement by various waved-induced mixing mechanisms is investigated and verified.

Global Mapping of Oceanic and Continental Shelf Crustal Thickness and Ocean-Continent Transition Structure

NASA Astrophysics Data System (ADS)

Kusznir, Nick; Alvey, Andy; Roberts, Alan

2017-04-01

The 3D mapping of crustal thickness for continental shelves and oceanic crust, and the determination of ocean-continent transition (OCT) structure and continent-ocean boundary (COB) location, represents a substantial challenge. Geophysical inversion of satellite derived free-air gravity anomaly data incorporating a lithosphere thermal anomaly correction (Chappell & Kusznir, 2008) now provides a useful and reliable methodology for mapping crustal thickness in the marine domain. Using this we have produced the first comprehensive maps of global crustal thickness for oceanic and continental shelf regions. Maps of crustal thickness and continental lithosphere thinning factor from gravity inversion may be used to determine the distribution of oceanic lithosphere, micro-continents and oceanic plateaux including for the inaccessible polar regions (e.g. Arctic Ocean, Alvey et al.,2008). The gravity inversion method provides a prediction of continent-ocean boundary location which is independent of ocean magnetic anomaly and isochron interpretation. Using crustal thickness and continental lithosphere thinning factor maps with superimposed shaded-relief free-air gravity anomaly, we can improve the determination of pre-breakup rifted margin conjugacy and sea-floor spreading trajectory during ocean basin formation. By restoring crustal thickness & continental lithosphere thinning to their initial post-breakup configuration we show the geometry and segmentation of the rifted continental margins at their time of breakup, together with the location of highly-stretched failed breakup basins and rifted micro-continents. For detailed analysis to constrain OCT structure, margin type (i.e. magma poor, "normal" or magma rich) and COB location, a suite of quantitative analytical methods may be used which include: (i) Crustal cross-sections showing Moho depth and crustal basement thickness from gravity inversion. (ii) Residual depth anomaly (RDA) analysis which is used to investigate OCT bathymetric anomalies with respect to expected oceanic values. This includes flexural backstripping to produce bathymetry corrected for sediment loading. (iii) Subsidence analysis which is used to determine the distribution of continental lithosphere thinning. (iv) Joint inversion of time-domain deep seismic reflection and gravity anomaly data which is used to determine lateral variations in crustal basement density and velocity across the OCT, and to validate deep seismic reflection interpretations of Moho depth. The combined interpretation of these independent quantitative measurements is used to determine crustal thickness and composition across the ocean-continent-transition. This integrated approach has been validated on the Iberian margin where ODP drilling provides ground-truth of ocean-continent-transition crustal structure, continent-ocean-boundary location and magmatic type.

ISEA (International geodetic project in SouthEastern Alaska) for rapid uplifting caused by glacial retreat: (4) Gravity tide observation

NASA Astrophysics Data System (ADS)

Sato, T.; Miura, S.; Sun, W.; Kaufman, A. M.; Cross, R.; Freymueller, J. T.; Heavner, M.

2006-12-01

The southeastern Alaska shows a large uplift rate as 30 mm/yr at most, which is considered to be closely related to the glacial isostatic adjustment (GIA) including two effects of the past and present-day ice melting (Larsen et al., 2004). So, this area is important to improve our knowledge of the viscoelastic property of the earth and to consider the global changes. Combing the displacement and gravity observations is useful to constrain the model computation results for GIA (Sato et al., 2006). In order to progress the previous work by the group of Univ. Alaska, Fairbanks (UAF), an observation project by Japan and USA groups was started in 2005 (Miura et al., this meeting). Under this project, June 2006, the continuous GPS measurements started (M. Kufman et al., this meeting) and the absolute gravity (AG) measurements were conducted (W. Sun et al., this meeting). Precise correction for the effect of ocean tide loading is one of the key to increase the observation accuracy of the GPS and gravity observations, especially for the AG measurement. Thanks for the satellite sea surface altimeters such as TOPEX/Poseidon and Jason-1, the accuracy of global ocean tide models based on these data has been much improved, and its accuracy is estimated at a level better than 1.3 cm as a RMS error of the vector differences of the 8 main tidal waves (Matsumoto et al., 2006). However, on the other hand, it is known that the southeastern Alaska is a place that shows a large discrepancy among the proposed global ocean tide models mainly due to a complex topography and bathymetry of the fjord area. In order to improve the accuracy of the ocean tide correction, we started the gravity tide observation at Juneau from June 2006. Two kinds of gravimeters are used for the observation. Sampling interval of the data is at every 1 min. We analyzed the 1 month data from the beginning of the observation and compared the tidal analysis results with the model tide including both effects of the solid and ocean tides. For this computation, we used the Love numbers and the loading Green function for the PREM earth model (Dziewonski & Anderson, 1981) and a global ocean tide model by Schwiderski (1980). Our comparison clearly indicates that a possibility to improve the accuracy of the model prediction by taking into account the actual tidal harmonics observed in the southeastern Alaska.

Global Distribution of Seamounts as Inferred from Ship Depth Soundings and Satellite Altimetry

NASA Astrophysics Data System (ADS)

Wessel, P.; Kim, S.; Sandwell, D. T.

2006-12-01

Traditionally, seamounts are active or extinct undersea volcanoes rising more than 1 km above the abyssal plain, but scientists now regularly apply the seamount label to features of just a few tens of meters in height. As constructional features they represent a small but significant fraction of the total volcanic extrusive budget for oceanic seafloor and their distribution provides key information on the variations in intraplate volcanic activity through space and time. Furthermore, they sustain significant ecological communities, determine habitats for fish, and act as obstacles to ocean currents, thus enhancing tidal energy dissipation and ocean mixing. Consequently, it is of some importance to locate and characterize seamounts. Two approaches are used to map the global distribution of seamounts. Depth soundings from single- and multi-beam echo sounders can provide the most detailed maps with up to 100--200 m horizontal resolution. However, soundings from the 5600 publicly available cruises sample only a small fraction of the ocean floor. Direct radar measurements of the ocean surface by satellite-borne altimeters have been used to infer the marine gravity field. By examining such gravity data one can characterize seamounts taller than ~2 km and such studies have produced seamount catalogues holding almost 15,000 seamounts. Recent retracking of the original radar altimeter waveforms to improve the accuracy of the gravity field has resulted in a two-fold increase in resolution. By extrapolating the inferred power-law that relates seamount size to frequency we estimate that 45,000 smaller seamounts taller than 1.5 km still remain uncharted. Future altimetry missions could improve on resolution and decrease noise levels even further, allowing for an even larger number of small (1--1.5 km) seamounts to be separated from the background abyssal hill fabric. Mapping the complete global distribution of seamounts will help constrain competing models of seamount formation as well as facilitate the understanding of marine habitats and deep ocean circulation.

Rigorous covariance propagation of geoid errors to geodetic MDT estimates

NASA Astrophysics Data System (ADS)

Pail, R.; Albertella, A.; Fecher, T.; Savcenko, R.

2012-04-01

The mean dynamic topography (MDT) is defined as the difference between the mean sea surface (MSS) derived from satellite altimetry, averaged over several years, and the static geoid. Assuming geostrophic conditions, from the MDT the ocean surface velocities as important component of global ocean circulation can be derived from it. Due to the availability of GOCE gravity field models, for the very first time MDT can now be derived solely from satellite observations (altimetry and gravity) down to spatial length-scales of 100 km and even below. Global gravity field models, parameterized in terms of spherical harmonic coefficients, are complemented by the full variance-covariance matrix (VCM). Therefore, for the geoid component a realistic statistical error estimate is available, while the error description of the altimetric component is still an open issue and is, if at all, attacked empirically. In this study we make the attempt to perform, based on the full gravity VCM, rigorous error propagation to derived geostrophic surface velocities, thus also considering all correlations. For the definition of the static geoid we use the third release of the time-wise GOCE model, as well as the satellite-only combination model GOCO03S. In detail, we will investigate the velocity errors resulting from the geoid component in dependence of the harmonic degree, and the impact of using/no using covariances on the MDT errors and its correlations. When deriving an MDT, it is spectrally filtered to a certain maximum degree, which is usually driven by the signal content of the geoid model, by applying isotropic or non-isotropic filters. Since this filtering is acting also on the geoid component, the consistent integration of this filter process into the covariance propagation shall be performed, and its impact shall be quantified. The study will be performed for MDT estimates in specific test areas of particular oceanographic interest.

Ocean Wave Studies with Applications to Ocean Modeling and Improvement of Satellite Altimeter Measurements

NASA Technical Reports Server (NTRS)

Glazman, Roman E.

1999-01-01

Combining analysis of satellite data (altimeter, scatterometer, high-resolution visible and infrared images, etc.) with mathematical modeling of non-linear wave processes, we investigate various ocean wave fields (on scales from capillary to planetary), their role in ocean dynamics and turbulent transport (of heat and biogeochemical quantities), and their effects on satellite altimeter measuring accuracy. In 1998 my attention was focused on long internal gravity waves (10 to 1000 km), known also as baroclinic inertia-gravity (BIG) waves. We found these waves to be a major factor of altimeter measurements "noise," resulting in a greater uncertainty [up to 10 cm in terms of sea surface height (SSH) amplitude] in the measured SSH signal than that caused by the sea state bias variations (up to 5 cm or so). This effect still remains largely overlooked by the satellite altimeter community. Our studies of BIG waves address not only their influence on altimeter measurements but also their role in global ocean dynamics and in transport and turbulent diffusion of biogeochemical quantities. In particular, in collaboration with Prof Peter Weichman, Caltech, we developed a theory of turbulent diffusion caused by wave motions of most general nature. Applied to the problem of horizontal turbulent diffusion in the ocean, the theory yielded the effective diffusion coefficient as a function of BIG wave parameters obtainable from satellite altimeter data. This effort, begun in 1997, has been successfully completed in 1998. We also developed a theory that relates spatial fluctuations of scalar fields (such as sea surface temperature, chlorophyll concentration, drifting ice concentration, etc.) to statistical characteristics of BIG waves obtainable from altimeter measurements. A manuscript is in the final stages of preparation. In order to verify the theoretical predictions and apply them to observations, we are now analyzing Sea-viewing Wide Field of view Sensor (SeaWiFS) and Field of view Sensor (SeaWiFS) and Advanced Very High-Resolution Radiometer (AVHRR) data on sea surface temperature (SST) and chlorophyll concentration jointly with TOPEX/POSEIDON data on SSH variations.

Tide Gauge Records Reveal Improved Processing of Gravity Recovery and Climate Experiment Time-Variable Mass Solutions over the Coastal Ocean

NASA Astrophysics Data System (ADS)

Piecuch, Christopher G.; Landerer, Felix W.; Ponte, Rui M.

2018-05-01

Monthly ocean bottom pressure solutions from the Gravity Recovery and Climate Experiment (GRACE), derived using surface spherical cap mass concentration (MC) blocks and spherical harmonics (SH) basis functions, are compared to tide gauge (TG) monthly averaged sea level data over 2003-2015 to evaluate improved gravimetric data processing methods near the coast. MC solutions can explain ≳ 42% of the monthly variance in TG time series over broad shelf regions and in semi-enclosed marginal seas. MC solutions also generally explain ˜5-32 % more TG data variance than SH estimates. Applying a coastline resolution improvement algorithm in the GRACE data processing leads to ˜ 31% more variance in TG records explained by the MC solution on average compared to not using this algorithm. Synthetic observations sampled from an ocean general circulation model exhibit similar patterns of correspondence between modeled TG and MC time series and differences between MC and SH time series in terms of their relationship with TG time series, suggesting that observational results here are generally consistent with expectations from ocean dynamics. This work demonstrates the improved quality of recent MC solutions compared to earlier SH estimates over the coastal ocean, and suggests that the MC solutions could be a useful tool for understanding contemporary coastal sea level variability and change.

Validation of Mean Absolute Sea Level of the North Atlantic obtained from Drifter, Altimetry and Wind Data

NASA Technical Reports Server (NTRS)

Maximenko, Nikolai A.

2003-01-01

Mean absolute sea level reflects the deviation of the Ocean surface from geoid due to the ocean currents and is an important characteristic of the dynamical state of the ocean. Values of its spatial variations (order of 1 m) are generally much smaller than deviations of the geoid shape from ellipsoid (order of 100 m) that makes the derivation of the absolute mean sea level a difficult task for gravity and satellite altimetry observations. Technique used by Niiler et al. for computation of the absolute mean sea level in the Kuroshio Extension was then developed into more general method and applied by Niiler et al. (2003b) to the global Ocean. The method is based on the consideration of balance of horizontal momentum.

Interaction of the Climate System and the Solid Earth: Analysis of Observations and Models

NASA Technical Reports Server (NTRS)

Bryan, Frank

2001-01-01

Under SENH funding we have carried out a number of diverse analyses of interactions of the climate system (atmosphere, ocean, land surface hydrology) with the solid Earth. While the original work plan emphasized analysis of excitation of variations in Earth rotation, with a lesser emphasis on time variable gravity, opportunities that developed during the proposal period in connection with preparations for the GRACE mission led us to a more balanced effort between these two topics. The results of our research are outlined in several topical sections: (1) oceanic excitation of variations in Earth rotation; (2) short period atmosphere-ocean excitation of variations in Earth rotation; (3) analysis of coupled climate system simulation; (4) observing system simulation studies for GRACE mission design; and (5) oceanic response to atmospheric pressure loading.

Assessing modern rates of river sediment discharge to the ocean using satellite gravimetry

NASA Astrophysics Data System (ADS)

Mouyen, Maxime; Longuevergne, Laurent; Steer, Philippe; Crave, Alain; Lemoine, Jean-Michel; Save, Himanshu; Robin, Cécile

2017-04-01

Worldwide rivers annually export about 19 Gigatons of sediments to the ocean that mostly accumulate in the coastal zones and on the continental shelves. This sediment discharge testifies of the intensity of continental erosion and records changes in climate, tectonics and human activity. However, natural and instrumental uncertainties inherent to the in-situ measurements of sediment discharge prevent from conclusive estimates to better understand these linkages. Here we develop a new method, using the Gravity Recovery and Climate Experiment (GRACE) satellite data, to infer mass-integrative estimates of sediment discharge of large rivers to the ocean. GRACE satellite provides global gravity time series that have proven useful for quantifying mass transport, including continental water redistribution at the Earth surface (ice sheets and glaciers melting, groundwater storage variations) but has been seldom used for monitoring sediment mass transfers so far. Here we pair the analysis of regularized GRACE solutions at high spatial resolution corrected from all known contributions (hydrology, ocean, atmosphere) to a particle tracking model that predicts the location of the sediment sinks for 13 rivers with the highest sediments loads in the world. We find that the resulting GRACE-derived sediment discharges off the mouth of the Amazon, Ganges-Brahmaputra, Changjiang (Yangtze), Indus, Magdalena, Godavari and Mekong rivers are consistent with in-situ measurements. Our results suggest that the lack of time continuity and of global coverage in terrestrial sediment discharge measurements could be reduced by using GRACE, which provides global and continuous data since 2002. GRACE solutions are regularly improved and new satellite gravity missions are being prepared hence making our approach even more relevant in a near future. The accumulation of sediments over time will keep increasing the signal to noise ratio of the gravity time series, which will improve the precision of the GRACE-derived sediment discharges values.

Acoustic Gravity Waves Generated by an Oscillating Ice Sheet in Arctic Zone

NASA Astrophysics Data System (ADS)

Abdolali, A.; Kadri, U.; Kirby, J. T., Jr.

2016-12-01

We investigate the formation of acoustic-gravity waves due to oscillations of large ice blocks, possibly triggered by atmospheric and ocean currents, ice block shrinkage or storms and ice-quakes.For the idealized case of a homogeneous weakly compressible water bounded at the surface by ice sheet and a rigid bed, the description of the infinite family of acoustic modes is characterized by the water depth h and angular frequency of oscillating ice sheet ω ; The acoustic wave field is governed by the leading mode given by: Nmax=\\floor {(ω h)/(π c)} where c is the sound speed in water and the special brackets represent the floor function (Fig1). Unlike the free-surface setting, the higher acoustic modes might exhibit a larger contribution and therefore all progressive acoustic modes have to be considered.This study focuses on the characteristics of acoustic-gravity waves generated by an oscillating elastic ice sheet in a weakly compressible fluid coupled with a free surface model [Abdolali et al. 2015] representing shrinking ice blocks in realistic sea state, where the randomly oriented ice sheets cause inter modal transition and multidirectional reflections. A theoretical solution and a 3D numerical model have been developed for the study purposes. The model is first validated against the theoretical solution [Kadri, 2016]. To overcome the computational difficulties of 3D models, we derive a depth-integrated equation valid for spatially varying ice sheet thickness and water depth. We show that the generated acoustic-gravity waves contribute significantly to deep ocean currents compared to other mechanisms. In addition, these waves travel at the sound speed in water carrying information on ice sheet motion, providing various implications for ocean monitoring and detection of ice-quakes. Fig1:Snapshots of dynamic pressure given by an oscillating ice sheet; h=4500m, c=1500m/s, semi-length b=10km, ζ =1m, omega=π rad/s. Abdolali, A., Kirby, J. T. and Bellotti, G., 2015, Depth-integrated equation for hydro-acoustic waves with bottom damping, Journal of Fluid Mechanics, 766, R1 doi:10.1017/jfm.2015.37 Kadri, U., 2016, Generation of Hydroacoustic Waves by an Oscillating Ice Block in Arctic Zones, Advances in Acoustics and Vibration. 2016. doi:10.1155/2016/8076108

Gravimetric investigations on the North American Datum (1972 - 1973)

NASA Technical Reports Server (NTRS)

Mather, R. S.

1975-01-01

All the available unclassified gravity data on the North American Datum (NAD) and in the surrounding oceans was assembled late in 1972 for the investigation of the gravity field in North America and its relation to North American Datum 1927 (NAD 27). The gravity data in Canada and the United States was compiled on a common datum compatible with the International Gravity Standardization Network 1971 (IGSN 71). The variation in the error of representation in the region is studied along with the correlation characteristics of gravity anomalies with elevation. A free air geoid (FAG 73) was computed from a combination of surface gravity data and Goddard Earth Model (GEM) 4 and this was used as the basis for the computation of the non-Stokesian contributions to the height anomaly. The geocentric orientation parameters obtained by this astrogravimetric method are compared with those obtained by satellite techniques. The differences are found to be no greater than those between individual satellite solutions. The differences between the astrogravimetric solution and satellite solutions GSFC 73 and GEM 6 are studied in detail with a view to obtaining a better understanding of these discrepancies.

Damping of short gravity-capillary waves due to oil derivatives film on the water surface

NASA Astrophysics Data System (ADS)

Sergievskaya, Irina; Ermakov, Stanislav; Lazareva, Tatyana

2016-10-01

In this paper new results of laboratory studies of damping of gravity-capillary waves on the water surface covered by kerosene are presented and compared with our previous analysis of characteristics of crude oil and diesel fuel films. Investigations of kerosene films were carried out in a wide range values of film thicknesses (from some hundreds millimetres to a few millimetres) and in a wide range of surface wave frequencies (from 10 to 27 Hz). The selected frequency range corresponds to the operating wavelengths of microwave, X- to Ka-band radars typically used for the ocean remote sensing. The studied range of film thickness covers typical thicknesses of routine spills in the ocean. It is obtained that characteristics of waves, measured in the presence of oil derivatives films differ from those for crude oil films, in particular, because the volume viscosity of oil derivatives and crude oil is strongly different. To retrieve parameters of kerosene films from the experimental data the surface wave damping was analyzed theoretically in the frame of a model of two-layer fluid. The films are assumed to be soluble, so the elasticity on the upper and lower boundaries is considered as a function of wave frequency. Physical parameters of oil derivative films were estimated when tuning the film parameters to fit theory and experiment. Comparison between wave damping due to crude oil, kerosene and diesel fuel films have shown some capabilities of distinguishing of oil films from remote sensing of short surface waves.

A geopotential model from satellite tracking, altimeter, and surface gravity data: GEM-T3

NASA Technical Reports Server (NTRS)

Lerch, F. J.; Nerem, R. S.; Putney, B. H.; Felsentreger, T. L.; Sanchez, B. V.; Marshall, J. A.; Klosko, S. M.; Patel, G. B.; Williamson, R. G.; Chinn, D. S.

1994-01-01

An improved model of Earth's gravitational field, Goddard Earth Model T-3 (GEM-T3), has been developed from a combination of satellite tracking, satellite altimeter, and surface gravimetric data. GEM-T3 provides a significant improvement in the modeling of the gravity field at half wavelengths of 400 km and longer. This model, complete to degree and order 50, yields more accurate satellite orbits and an improved geoid representation than previous Goddard Earth Models. GEM-T3 uses altimeter data from GEOS 3 (1975-1976), Seasat (1978) and Geosat (1986-1987). Tracking information used in the solution includes more than 1300 arcs of data encompassing 31 different satellites. The recovery of the long-wavelength components of the solution relies mostly on highly precise satellite laser ranging (SLR) data, but also includes Tracking Network (TRANET) Doppler, optical, and satellite-to-satellite tracking acquired between the ATS 6 and GEOS 3 satellites. The main advances over GEM-T2 (beyond the inclusion of altimeter and surface gravity information which is essential for the resolution of the shorter wavelength geoid) are some improved tracking data analysis approaches and additional SLR data. Although the use of altimeter data has greatly enhanced the modeling of the ocean geoid between 65 deg N and 60 deg S latitudes in GEM-T3, the lack of accurate detailed surface gravimetry leaves poor geoid resolution over many continental regions of great tectonic interest (e.g., Himalayas, Andes). Estimates of polar motion, tracking station coordinates, and long-wavelength ocean tidal terms were also made (accounting for 6330 parameters). GEM-T3 has undergone error calibration using a technique based on subset solutions to produce reliable error estimates. The calibration is based on the condition that the expected mean square deviation of a subset gravity solution from the full set values is predicted by the solutions' error covariances. Data weights are iteratively adjusted until this condition for the error calibration is satisfied. In addition, gravity field tests were performed on strong satellite data sets withheld from the solution (thereby ensuring their independence). In these tests, the performance of the subset models on the withheld observations is compared to error projections based on their calibrated error covariances. These results demonstrate that orbit accuracy projections are reliable for new satellites which were not included in GEM-T3.

Tsunami normal modes with solid earth and atmospheric coupling and inversion of the TEC data to estimate tsunami water height in the case of the Queen Charlotte tsunami.

NASA Astrophysics Data System (ADS)

Rakoto, V.; Lognonne, P. H.; Rolland, L.

2016-12-01

Large underwater earthquakes (Mw > 7) can transmit part of their energy to the surrounding ocean through large sea-floor motions, generating tsunamis that propagate over long distances. The forcing effect of long period ocean surface vibrations due to tsunami waves on the atmosphere trigger atmospheric internal gravity waves (IGWs) that induce ionospheric disturbances when they reach the upper atmosphere. In this poster, we study the IGWs associated to tsunamis using a normal modes 1D modeling approach. Our model is first applied to the case of the October 2012 Haida Gwaii tsunami observed offshore Hawaii. We found three resonances between tsunami modes and the atmospheric gravity modes occurring around 1.5 mHz, 2 mHz and 2.5 mHz, with a large fraction of the energy of the tsunami modes transferred from the ocean to the atmosphere. At theses frequencies, the gravity branches are interacting with the tsunami one and have large amplitude in the ocean. As opposed to the tsunami, a fraction of their energy is therefore transferred from the atmosphere to the ocean. We also show that the fundamental of the gravity waves should arrive before the tsunami due to higher group velocity below 1.6 mHz. We demonstrate that only the 1.5 mHz resonance of the tsunami mode can trigger observable ionospheric perturbations, most often monitored using GPS dual-frequency measurements. Indeed, we show that the modes at 2 mHz and 2.5 mHz are already evanescent at the height of the F2 peak and have little energy in the ionosphere. This normal modes modeling offers a novel and comprehensive study of the transfer function from a propagating tsunami to the upper atmosphere. In particular, we can invert the perturbed TEC data induced by a tsunami in order to estimate the amplitude of the tsunami waveform using a least square method. This method has been performed in the case of the Haida Gwaii tsunami. The results showed a good agreement with the measurement of the dart buoy.

Gravity Maps of Antarctic Lithospheric Structure from Remote-Sensing and Seismic Data

NASA Astrophysics Data System (ADS)

Tenzer, Robert; Chen, Wenjin; Baranov, Alexey; Bagherbandi, Mohammad

2018-02-01

Remote-sensing data from altimetry and gravity satellite missions combined with seismic information have been used to investigate the Earth's interior, particularly focusing on the lithospheric structure. In this study, we use the subglacial bedrock relief BEDMAP2, the global gravitational model GOCO05S, and the ETOPO1 topographic/bathymetric data, together with a newly developed (continental-scale) seismic crustal model for Antarctica to compile the free-air, Bouguer, and mantle gravity maps over this continent and surrounding oceanic areas. We then use these gravity maps to interpret the Antarctic crustal and uppermost mantle structure. We demonstrate that most of the gravity features seen in gravity maps could be explained by known lithospheric structures. The Bouguer gravity map reveals a contrast between the oceanic and continental crust which marks the extension of the Antarctic continental margins. The isostatic signature in this gravity map confirms deep and compact orogenic roots under the Gamburtsev Subglacial Mountains and more complex orogenic structures under Dronning Maud Land in East Antarctica. Whereas the Bouguer gravity map exhibits features which are closely spatially correlated with the crustal thickness, the mantle gravity map reveals mainly the gravitational signature of the uppermost mantle, which is superposed over a weaker (long-wavelength) signature of density heterogeneities distributed deeper in the mantle. In contrast to a relatively complex and segmented uppermost mantle structure of West Antarctica, the mantle gravity map confirmed a more uniform structure of the East Antarctic Craton. The most pronounced features in this gravity map are divergent tectonic margins along mid-oceanic ridges and continental rifts. Gravity lows at these locations indicate that a broad region of the West Antarctic Rift System continuously extends between the Atlantic-Indian and Pacific-Antarctic mid-oceanic ridges and it is possibly formed by two major fault segments. Gravity lows over the Transantarctic Mountains confirms their non-collisional origin. Additionally, more localized gravity lows closely coincide with known locations of hotspots and volcanic regions (Marie Byrd Land, Balleny Islands, Mt. Erebus). Gravity lows also suggest a possible hotspot under the South Orkney Islands. However, this finding has to be further verified.

Gravity model for the North Atlantic ocean mantle: results, uncertainties and links to regional geodynamics

NASA Astrophysics Data System (ADS)

Barantsrva, O.; Artemieva, I. M.; Thybo, H.

2015-12-01

We present the results of gravity modeling for the North Atlantic region based on interpretation of GOCE gravity satellite data. First, to separate the gravity signal caused by density anomalies within the crust and the upper mantle, we subtract the lower harmonics in the gravity field, which are presumably caused by deep density structure of the Earth (the core and the lower mantle). Next, the gravity effect of the upper mantle is calculated by subtracting the gravity effect of the crustal model. Our "basic model" is constrained by a recent regional seismic model EUNAseis for the crustal structure (Artemieva and Thybo, 2013); for bathymetry and topography we use a global ETOPO1 model by NOAA. We test sensitivity of the results to different input parameters, such as bathymetry, crustal structure, and gravity field. For bathymetry, we additionally use GEBCO data; for crustal correction - a global model CRUST 1.0 (Laske, 2013); for gravity - EGM2008 (Pavlis, 2012). Sensitivity analysis shows that uncertainty in the crustal structure produces the largest deviation from "the basic model". Use of different bathymetry data has little effect on the final results, comparable to the interpolation error. The difference in mantle residual gravity models based on GOCE and EMG2008 gravity data is 5-10 mGal. The results based on two crustal models have a similar pattern, but differ significantly in amplitude (ca. 250 mGal) for the Greenland-Faroe Ridge. The results demonstrate the presence of a strong gravity and density heterogeneity in the upper mantle in the North Atlantic region. A number of mantle residual gravity anomalies are robust features, independent of the choice of model parameters. This include (i) a sharp contrast at the continent-ocean transition, (ii) positive mantle gravity anomalies associated with continental fragments (microcontinents) in the North Atlantic ocean; (iii) negative mantle gravity anomalies which mark regions with anomalous oceanic mantle and the Mid-Atlantic Ridge. To understand better a complex geodynamics mosaic in the region, we compare our results with regional geochemical data (Korenaga and Klemen, 2000), and find that residual mantle gravity anomalies are well correlated with anomalies in epsilon-Nd and iron-depletion.

Spacebased Observation of Water Balance Over Global Oceans

NASA Astrophysics Data System (ADS)

Liu, W.; Xie, X.

2008-12-01

We demonstrated that ocean surface fresh water flux less the water discharge into the ocean from river and ice melt balances the mass loss in the ocean both in magnitude and in the phase of annual variation. The surface water flux was computed from the divergence of the water transport integrated over the depth of the atmosphere. The atmospheric water transport is estimated from the precipitable water measured by Special Sensor Microwave Imager, the surface wind vector by QuikSCAT, and the NOAA cloud drift wind through a statistical model. The transport has been extensively validated using global radiosonde and data and operational numerical weather prediction results. Its divergence has been shown to agree with the difference between evaporation estimated from the Advanced Microwave Scanning Radiometer data and the precipitation measured by Tropical Rain Measuring Mission over the global tropical and subtropical oceans both in magnitude and geographical distribution for temporal scales ranging from intraseasonal to interannual. The water loss rate in the ocean is estimated by two methods, one is from Gravity Recovery and Climate Experiment and the other is by subtracting the climatological steric change from the sea level change measured by radar altimeter on Jason. Only climatological river discharge and ice melt from in situ measurements are available and the lack of temporal variation may contribute to discrepancies in the balance. We have successfully used the spacebased surface fluxes to estimate to climatological mean heat transport in the Atlantic ocean and is attempting to estimate the meridional fresh water (or salt) transport from the surface flux. The approximate closure of the water balance gives a powerful indirect validation of the spacebased products.

Worldwide complete spherical Bouguer and isostatic anomaly maps

NASA Astrophysics Data System (ADS)

Bonvalot, S.; Balmino, G.; Briais, A.; Peyrefitte, A.; Vales, N.; Biancale, R.; Gabalda, G.; Reinquin, F.

2011-12-01

We present here a set of digital maps of the Earth's gravity anomalies (surface "free air", Bouguer and isostatic), computed at Bureau Gravimetric International (BGI) as a contribution to the Global Geodetic Observing Systems (GGOS) and to the global geophysical maps published by the Commission for the Geological Map of the World (CGMW). The free air and Bouguer anomaly concept is extensively used in geophysical interpretation to investigate the density distributions in the Earth's interior. Complete Bouguer anomalies (including terrain effects) are usually computed at regional scales by integrating the gravity attraction of topography elements over and beyond a given area (under planar or spherical approximations). Here, we developed and applied a worldwide spherical approach aimed to provide a set of homogeneous and high resolution gravity anomaly maps and grids computed at the Earth's surface, taking into account a realistic Earth model and reconciling geophysical and geodetic definitions of gravity anomalies. This first version (1.0) has been computed by spherical harmonics analysis / synthesis of the Earth's topography-bathymetry up to degree 10800. The detailed theory of the spherical harmonics approach is given in Balmino et al., (Journal of Geodesy, submitted). The Bouguer and terrain corrections have thus been computed in spherical geometry at 1'x1' resolution using the ETOPO1 topography/bathymetry, ice surface and bedrock models from the NOAA (National Oceanic and Atmospheric Administration) and taking into account precise characteristics (boundaries and densities) of major lakes, inner seas, polar caps and of land areas below sea level. Isostatic corrections have been computed according to the Airy Heiskanen model in spherical geometry for a constant depth of compensation of 30km. The gravity information given here is provided by the Earth Geopotential Model (EGM2008), developed at degree 2160 by the National Geospatial Intelligence Agency (NGA) (Pavlis et al., 2008), which represents the best up-to-date global gravity model (including surface gravity measurements from land, marine and airborne surveys as well as gravity and altimetry satellite measurements). The surface gravity anomaly (free air) is computed at the Earth's surface in the context of Molodensky theory and includes corrections from the mass of the atmosphere. The way gravity anomalies are computed on a worldwide basis slightly differs from the classical usage, but meets modern concerns which tend to take the real Earth into account. The resulting anomaly maps and grids will be distributed for scientific and education purposes by the Commission for the Geological Map of the World (CGMW) with support of UNESCO and other institutions. Upgraded versions might be done as soon as new global gravity model is available (including satellite GOCE and new surface measurements: ground, airborne). Visit / contact BGI (http://bgi.omp.obs-mip.fr) and CCMW (http://ccgm.free.fr) for more information.

Ocean-Wave Coupled Modeling in COAMPS-TC: A Study of Hurricane Ivan (2004)

DTIC Science & Technology

2013-08-15

America, Stennis Space Center, MS 39529, United States a r t i c l e i n f oArticle history : Received 12 December 2012 Received in revised form 16 May 2013...Information for Society. Venice , Italy, 21– 25 September 2009. ESA Publication. Sullivan, P.P., McWilliams, J.C., Melville, W.K., 2007. Surface gravity

Optimal Geoid Modelling to determine the Mean Ocean Circulation - Project Overview and early Results

NASA Astrophysics Data System (ADS)

Fecher, Thomas; Knudsen, Per; Bettadpur, Srinivas; Gruber, Thomas; Maximenko, Nikolai; Pie, Nadege; Siegismund, Frank; Stammer, Detlef

2017-04-01

The ESA project GOCE-OGMOC (Optimal Geoid Modelling based on GOCE and GRACE third-party mission data and merging with altimetric sea surface data to optimally determine Ocean Circulation) examines the influence of the satellite missions GRACE and in particular GOCE in ocean modelling applications. The project goal is an improved processing of satellite and ground data for the preparation and combination of gravity and altimetry data on the way to an optimal MDT solution. Explicitly, the two main objectives are (i) to enhance the GRACE error modelling and optimally combine GOCE and GRACE [and optionally terrestrial/altimetric data] and (ii) to integrate the optimal Earth gravity field model with MSS and drifter information to derive a state-of-the art MDT including an error assessment. The main work packages referring to (i) are the characterization of geoid model errors, the identification of GRACE error sources, the revision of GRACE error models, the optimization of weighting schemes for the participating data sets and finally the estimation of an optimally combined gravity field model. In this context, also the leakage of terrestrial data into coastal regions shall be investigated, as leakage is not only a problem for the gravity field model itself, but is also mirrored in a derived MDT solution. Related to (ii) the tasks are the revision of MSS error covariances, the assessment of the mean circulation using drifter data sets and the computation of an optimal geodetic MDT as well as a so called state-of-the-art MDT, which combines the geodetic MDT with drifter mean circulation data. This paper presents an overview over the project results with focus on the geodetic results part.

Spectral evolution and extreme value analysis of non-linear numerical simulations of narrow band random surface gravity waves.

NASA Astrophysics Data System (ADS)

Socquet-Juglard, H.; Dysthe, K. B.; Trulsen, K.; Liu, J.; Krogstad, H. E.

2003-04-01

Numerical simulations of a narrow band gaussian spectrum of random surface gravity waves have been carried out in two and three spatial dimensions [7]. Different types of non-linear Schr&{uml;o}dinger equations, [1] and [4], have been used in these simulations. Simulations have now been carried with a JONSWAP spectrum associated with a spreading function of the type cosine-squared [5]. The evolution of the spectrum, skewness, kurtosis, ... will be presented. In addition, some results about stochastic properties of the surface will be shown. Based on the approach found in [2], [3] and [6], the results are presented in terms of deviations from linear Gaussian theory and the standard second order small slope perturbation theory. begin{thebibliography}{9} bibitem{kk96} Trulsen, K. &Dysthe, K. B. (1996). A modified nonlinear Schr&{uml;o}dinger equation for broader bandwidth gravity waves on deep water. Wave Motion, 24, pp. 281-289. bibitem{BK2000} Krogstad, H.E. and S.F. Barstow (2000). A uniform approach to extreme value analysis of ocean waves, Proc. ISOPE'2000, Seattle, USA, 3, pp. 103-108. bibitem{PRK} Prevosto, M., H. E. Krogstad and A. Robin (2000). Probability distributions for maximum wave and crest heights, Coast. Eng., 40, 329-360. bibitem{ketal} Trulsen, K., Kliakhandler, I., Dysthe, K. B. &Velarde, M. G. (2000) On weakly nonlinear modulation of waves on deep water, Phys. Fluids, 12, pp. L25-L28. bibitem{onorato} Onorato, M., Osborne, A.R. and Serio, M. (2002) Extreme wave events in directional, random oceanic sea states, Phys. Fluids, 14, pp. 2432-2437. bibitem{BK2002} Krogstad, H.E. and S.F. Barstow (2002). Analysis and Applications of Second Order Models for the Maximum Crest height, % Proc. 21nd Int. Conf. Offshore Mechanics and Arctic Engineering, Oslo. Paper no. OMAE2002-28479. bibitem{JFMP} Dysthe, K. B., Trulsen, K., Krogstad, H. E. and Socquet-Juglard, H. (2002, in press) Evolution of a narrow band spectrum of random surface gravity waves, J. Fluid Mech.

GRAVITY DATA OBTAINED DURING CHAIN CRUISE NO. 70.

DTIC Science & Technology

Profiles, tabulations, and charts are presented of principal gravity facts, free-air gravity anomalies, and simple Bouguer gravity anomalies obtained in the western North Atlantic Ocean during R/V CHAIN cruise 70. (Author)

Sea surface determination from space: The GSFC geoid

NASA Technical Reports Server (NTRS)

Vonbun, F. O.; Mcgoogan, J.; Marsh, J.; Lerch, F. J.

1975-01-01

The determination of the sea surface/geoid and its relative variation were investigated and results of the altimeter experiment on Skylab to test the geoid are discussed. The spaceborne altimeter on Skylab revealed that the sea surface of the world's oceans can be measured with an accuracy in the meter range. Surface variations are discussed as they relate to those computed from satellite orbital dynamics and ground based gravity data. The GSFC geoid was constructed from about 400,000 satellite tracking data (range, range rate, angles) and about 20,000 ground gravity observations. One of the last experiments on Skylab was to measure and/or test this geoid over almost one orbit. It was found that the computed water surface deviates between 5 to 20 m from the measured one. Further outlined are the influence of orbital errors on the sea surface, and numerical examples are given based upon real tracking data. Orbital height error estimates were computed for geodetic type satellites and are found to be in the order of 0.2 to 5 meters.

Estimation and Validation of Oceanic Mass Circulation from the GRACE Mission

NASA Technical Reports Server (NTRS)

Boy, J.-P.; Rowlands, D. D.; Sabaka, T. J.; Luthcke, S. B.; Lemoine, F. G.

2011-01-01

Since the launch of the Gravity Recovery And Climate Experiment (GRACE) in March 2002, the Earth's surface mass variations have been monitored with unprecedented accuracy and resolution. Compared to the classical spherical harmonic solutions, global high-resolution mascon solutions allows the retrieval of mass variations with higher spatial and temporal sampling (2 degrees and 10 days). We present here the validation of the GRACE global mascon solutions by comparing mass estimates to a set of about 100 ocean bottom pressure (OSP) records, and show that the forward modelling of continental hydrology prior to the inversion of the K-band range rate data allows better estimates of ocean mass variations. We also validate our GRACE results to OSP variations modelled by different state-of-the-art ocean general circulation models, including ECCO (Estimating the Circulation and Climate of the Ocean) and operational and reanalysis from the MERCATOR project.

Spreading rate dependence of gravity anomalies along oceanic transform faults.

PubMed

Gregg, Patricia M; Lin, Jian; Behn, Mark D; Montési, Laurent G J

2007-07-12

Mid-ocean ridge morphology and crustal accretion are known to depend on the spreading rate of the ridge. Slow-spreading mid-ocean-ridge segments exhibit significant crustal thinning towards transform and non-transform offsets, which is thought to arise from a three-dimensional process of buoyant mantle upwelling and melt migration focused beneath the centres of ridge segments. In contrast, fast-spreading mid-ocean ridges are characterized by smaller, segment-scale variations in crustal thickness, which reflect more uniform mantle upwelling beneath the ridge axis. Here we present a systematic study of the residual mantle Bouguer gravity anomaly of 19 oceanic transform faults that reveals a strong correlation between gravity signature and spreading rate. Previous studies have shown that slow-slipping transform faults are marked by more positive gravity anomalies than their adjacent ridge segments, but our analysis reveals that intermediate and fast-slipping transform faults exhibit more negative gravity anomalies than their adjacent ridge segments. This finding indicates that there is a mass deficit at intermediate- and fast-slipping transform faults, which could reflect increased rock porosity, serpentinization of mantle peridotite, and/or crustal thickening. The most negative anomalies correspond to topographic highs flanking the transform faults, rather than to transform troughs (where deformation is probably focused and porosity and alteration are expected to be greatest), indicating that crustal thickening could be an important contributor to the negative gravity anomalies observed. This finding in turn suggests that three-dimensional magma accretion may occur near intermediate- and fast-slipping transform faults.

Ocean Thermal and Color Evolution During the 1997/1998 ENSO Event

NASA Technical Reports Server (NTRS)

Rienecker, Michele

1998-01-01

A reduced gravity primitive equation modeling and assimilation system is used to study the evolution of the tropical Pacific during the 1997/1998 ENSO cycle. The modeling/assimilation scheme ingests satellite altimeter data and TAO temperature profiles and uses SSM/I satellite derived winds as surface boundary forcing. The four-dimensional structure of the upper ocean circulation structure will be compared against available in situ observations across the Pacific basin. In particular, variability near the Galapagos Islands will be highlighted during the spring of 1998 when phytoplankton concentrations were observed to increase a hundred-fold over a two week period.

Three-wave and four-wave interactions in gravity wave turbulence

NASA Astrophysics Data System (ADS)

Aubourg, Quentin; Campagne, Antoine; Peureux, Charles; Ardhuin, Fabrice; Sommeria, Joel; Viboud, Samuel; Mordant, Nicolas

2017-11-01

Weak-turbulence theory is a statistical framework to describe a large ensemble of nonlinearly interacting waves. The archetypal example of such system is the ocean surface that is made of interacting surface gravity waves. Here we describe a laboratory experiment dedicated to probe the statistical properties of turbulent gravity waves. We set up an isotropic state of interacting gravity waves in the Coriolis facility (13-m-diam circular wave tank) by exciting waves at 1 Hz by wedge wave makers. We implement a stereoscopic technique to obtain a measurement of the surface elevation that is resolved in both space and time. Fourier analysis shows that the laboratory spectra are systematically steeper than the theoretical predictions and the field observations in the Black Sea by Leckler et al. [F. Leckler et al., J. Phys. Oceanogr. 45, 2484 (2015), 10.1175/JPO-D-14-0237.1]. We identify a strong impact of surface dissipation on the scaling of the Fourier spectrum at the scales that are accessible in the experiments. We use bicoherence and tricoherence statistical tools in frequency and/or wave-vector space to identify the active nonlinear coupling. These analyses are also performed on the field data by Leckler et al. for comparison with the laboratory data. Three-wave coupling is characterized by and shown to involve mostly quasiresonances of waves with second- or higher-order harmonics. Four-wave coupling is not observed in the laboratory but is evidenced in the field data. We discuss temporal scale separation to explain our observations.

Anomalous Structure of Oceanic Lithosphere in the North Atlantic and Arctic Oceans: A Preliminary Analysis Based on Bathymetry, Gravity and Crustal Structure

NASA Astrophysics Data System (ADS)

Barantsrva, O.

2014-12-01

We present a preliminary analysis of the crustal and upper mantle structure for off-shore regions in the North Atlantic and Arctic oceans. These regions have anomalous oceanic lithosphere: the upper mantle of the North Atlantic ocean is affected by the Iceland plume, while the Arctic ocean has some of the slowest spreading rates. Our specific goal is to constrain the density structure of the upper mantle in order to understand the links between the deep lithosphere dynamics, ocean spreading, ocean floor bathymetry, heat flow and structure of the oceanic lithosphere in the regions where classical models of evolution of the oceanic lithosphere may not be valid. The major focus is on the oceanic lithosphere, but the Arctic shelves with a sufficient data coverage are also included into the analysis. Out major interest is the density structure of the upper mantle, and the analysis is based on the interpretation of GOCE satellite gravity data. To separate gravity anomalies caused by subcrustal anomalous masses, the gravitational effect of water, crust and the deep mantle is removed from the observed gravity field. For bathymetry we use the global NOAA database ETOPO1. The crustal correction to gravity is based on two crustal models: (1) global model CRUST1.0 (Laske, 2013) and, for a comparison, (2) a regional seismic model EUNAseis (Artemieva and Thybo, 2013). The crustal density structure required for the crustal correction is constrained from Vp data. Previous studies have shown that a large range of density values corresponds to any Vp value. To overcome this problem and to reduce uncertainty associated with the velocity-density conversion, we account for regional tectonic variations in the Northern Atlantics as constrained by numerous published seismic profiles and potential-field models across the Norwegian off-shore crust (e.g. Breivik et al., 2005, 2007), and apply different Vp-density conversions for different parts of the region. We present preliminary results, which we use to examine factors that control variations in bathymetry, sedimentary and crustal thicknesses in these anomalous oceanic domains.

Gravity Field Changes due to Long-Term Sea Level Changes

NASA Astrophysics Data System (ADS)

Makarynskyy, O.; Kuhn, M.; Featherstone, W. E.

2004-12-01

Long-term sea level changes caused by climatic changes (e.g. global warming) will alter the system Earth. This includes the redistribution of ocean water masses due to the migration of cold fresh water from formerly ice-covered regions to the open oceans mainly caused by the deglaciation of polar ice caps. Consequently also a change in global ocean circulation patterns will occur. Over a longer timescale, such mass redistributions will be followed by isostatic rebound/depression due to the changed surface un/loading, resulting in variable sea level change around the world. These, in turn, will affect the gravity field, location of the geocentre, and the Earth's rotation vector. This presentation focuses mainly on gravity field changes induced by long-term (hundredths to many thousand years) sea level changes using an Earth System Climate Model (ESCM) of intermediate complexity. In this study, the coupled University of Victoria (Victoria, Canada) Earth System Climate Model (Uvic ESCM) was used, which embraces the primary thermodynamic and hydrological components of the climate system including sea and land-ice information. The model was implemented to estimate changes in global precipitation, ocean mass redistribution, seawater temperature and salinity on timescales from hundreds to thousands years under different greenhouse warming scenarios. The sea level change output of the model has been converted into real mass changes by removing the steric effect, computed from seawater temperature and salinity information at different layers also provided by Uvic ESCM. Finally the obtained mass changes have been converted into changes of the gravitational potential and subsequently of the geoid height using a spherical harmonic representation of the different data. Preliminary numerical results are provided for sea level change as well as change in geoid height.

Treatment of ocean tide aliasing in the context of a next generation gravity field mission

NASA Astrophysics Data System (ADS)

Hauk, Markus; Pail, Roland

2018-07-01

Current temporal gravity field solutions from Gravity Recovery and Climate Experiment (GRACE) suffer from temporal aliasing errors due to undersampling of signal to be recovered (e.g. hydrology), uncertainties in the de-aliasing models (usually atmosphere and ocean) and imperfect ocean tide models. Especially the latter will be one of the most limiting factors in determining high-resolution temporal gravity fields from future gravity missions such as GRACE Follow-On and Next-Generation Gravity Missions (NGGM). In this paper a method to co-parametrize ocean tide parameters of the eight main tidal constituents over time spans of several years is analysed and assessed. Numerical closed-loop simulations of low-low satellite-to-satellite-tracking missions for a single polar pair and a double pair Bender-type formation are performed, using time variable geophysical background models and noise assumptions for new generation instrument technology. Compared to the single pair mission, results show a reduction of tide model errors up to 70 per cent for dedicated tidal constituents due to an enhanced spatial and temporal sampling and error isotropy for the double pair constellation. Extending the observation period from 1 to 3 yr leads to a further reduction of tidal errors up to 60 per cent for certain constituents, and considering non-tidal mass changes during the estimation process leads to reductions of tidal errors between 20 and 80 per cent. As part of a two-step approach, the estimated tide model is used for de-aliasing during gravity field retrieval in a second iteration, resulting in more than 50 per cent reduction of ocean tide aliasing errors for a NGGM Bender-type formation.

Determining the Ocean's Role on the Variable Gravity Field on Earth Rotation

NASA Technical Reports Server (NTRS)

Ponte, Rui M.

1999-01-01

A number of ocean models of different complexity have been used to study changes in the oceanic mass field and angular momentum and their relation to the variable Earth rotation and gravity field. Time scales examined range from seasonal to a few days. Results point to the importance of oceanic signals in driving polar motion, in particular the Chandler and annual wobbles. Results also show that oceanic signals have a measurable impact on length-of-day variations. Various circulation features and associated mass signals, including the North Pacific subtropical gyre, the equatorial currents, and the Antarctic Circumpolar Current play a significant role in oceanic angular momentum variability.

Transport of inertial anisotropic particles under surface gravity waves

NASA Astrophysics Data System (ADS)

Dibenedetto, Michelle; Koseff, Jeffrey; Ouellette, Nicholas

2016-11-01

The motion of neutrally and almost-neutrally buoyant particles under surface gravity waves is relevant to the transport of microplastic debris and other small particulates in the ocean. Consequently, a number of studies have looked at the transport of spherical particles or mobile plankton in these conditions. However, the effects of particle-shape anisotropy on the trajectories and behavior of irregularly shaped particles in this type of oscillatory flow are still relatively unknown. To better understand these issues, we created an idealized numerical model which simulates the three-dimensional behavior of anisotropic spheroids in flow described by Airy wave theory. The particle's response is calculated using a simplified Maxey-Riley equation coupled with Jeffery's equation for particle rotation. We show that the particle dynamics are strongly dependent on their initial conditions and shape, with some some additional dependence on Stokes number.

Reorientation of Sputnik Planitia implies a subsurface ocean on Pluto.

PubMed

Nimmo, F; Hamilton, D P; McKinnon, W B; Schenk, P M; Binzel, R P; Bierson, C J; Beyer, R A; Moore, J M; Stern, S A; Weaver, H A; Olkin, C B; Young, L A; Smith, K E

2016-12-01

The deep nitrogen-covered basin on Pluto, informally named Sputnik Planitia, is located very close to the longitude of Pluto's tidal axis and may be an impact feature, by analogy with other large basins in the Solar System. Reorientation of Sputnik Planitia arising from tidal and rotational torques can explain the basin's present-day location, but requires the feature to be a positive gravity anomaly, despite its negative topography. Here we argue that if Sputnik Planitia did indeed form as a result of an impact and if Pluto possesses a subsurface ocean, the required positive gravity anomaly would naturally result because of shell thinning and ocean uplift, followed by later modest nitrogen deposition. Without a subsurface ocean, a positive gravity anomaly requires an implausibly thick nitrogen layer (exceeding 40 kilometres). To prolong the lifetime of such a subsurface ocean to the present day and to maintain ocean uplift, a rigid, conductive water-ice shell is required. Because nitrogen deposition is latitude-dependent, nitrogen loading and reorientation may have exhibited complex feedbacks.

Reorientation of Sputnik Planitia implies a subsurface ocean on Pluto

NASA Astrophysics Data System (ADS)

Nimmo, F.; Hamilton, D. P.; McKinnon, W. B.; Schenk, P. M.; Binzel, R. P.; Bierson, C. J.; Beyer, R. A.; Moore, J. M.; Stern, S. A.; Weaver, H. A.; Olkin, C. B.; Young, L. A.; Smith, K. E.; Moore, J. M.; McKinnon, W. B.; Spencer, J. R.; Beyer, R.; Binzel, R. P.; Buie, M.; Buratti, B.; Cheng, A.; Cruikshank, D.; Ore, C. Dalle; Earle, A.; Gladstone, R.; Grundy, W.; Howard, A. D.; Lauer, T.; Linscott, I.; Nimmo, F.; Parker, J.; Porter, S.; Reitsema, H.; Reuter, D.; Roberts, J. H.; Robbins, S.; Schenk, P. M.; Showalter, M.; Singer, K.; Strobel, D.; Summers, M.; Tyler, L.; White, O. L.; Umurhan, O. M.; Banks, M.; Barnouin, O.; Bray, V.; Carcich, B.; Chaikin, A.; Chavez, C.; Conrad, C.; Hamilton, D. P.; Howett, C.; Hofgartner, J.; Kammer, J.; Lisse, C.; Marcotte, A.; Parker, A.; Retherford, K.; Saina, M.; Runyon, K.; Schindhelm, E.; Stansberry, J.; Steffl, A.; Stryk, T.; Throop, H.; Tsang, C.; Verbiscer, A.; Winters, H.; Zangari, A.; Stern, S. A.; Weaver, H. A.; Olkin, C. B.; Young, L. A.; Smith, K. E.

2016-12-01

The deep nitrogen-covered basin on Pluto, informally named Sputnik Planitia, is located very close to the longitude of Pluto’s tidal axis and may be an impact feature, by analogy with other large basins in the Solar System. Reorientation of Sputnik Planitia arising from tidal and rotational torques can explain the basin’s present-day location, but requires the feature to be a positive gravity anomaly, despite its negative topography. Here we argue that if Sputnik Planitia did indeed form as a result of an impact and if Pluto possesses a subsurface ocean, the required positive gravity anomaly would naturally result because of shell thinning and ocean uplift, followed by later modest nitrogen deposition. Without a subsurface ocean, a positive gravity anomaly requires an implausibly thick nitrogen layer (exceeding 40 kilometres). To prolong the lifetime of such a subsurface ocean to the present day and to maintain ocean uplift, a rigid, conductive water-ice shell is required. Because nitrogen deposition is latitude-dependent, nitrogen loading and reorientation may have exhibited complex feedbacks.

Low degree Earth's gravity coefficients determined from different space geodetic observations and climate models

NASA Astrophysics Data System (ADS)

Wińska, Małgorzata; Nastula, Jolanta

2017-04-01

Large scale mass redistribution and its transport within the Earth system causes changes in the Earth's rotation in space, gravity field and Earth's ellipsoid shape. These changes are observed in the ΔC21, ΔS21, and ΔC20 spherical harmonics gravity coefficients, which are proportional to the mass load-induced Earth rotational excitations. In this study, linear trend, decadal, inter-annual, and seasonal variations of low degree spherical harmonics coefficients of Earth's gravity field, determined from different space geodetic techniques, Gravity Recovery and Climate Experiment (GRACE), satellite laser ranging (SLR), Global Navigation Satellite System (GNSS), Earth rotation, and climate models, are examined. In this way, the contribution of each measurement technique to interpreting the low degree surface mass density of the Earth is shown. Especially, we evaluate an usefulness of several climate models from the Coupled Model Intercomparison Project phase 5 (CMIP5) to determine the low degree Earth's gravity coefficients using GRACE satellite observations. To do that, Terrestrial Water Storage (TWS) changes from several CMIP5 climate models are determined and then these simulated data are compared with the GRACE observations. Spherical harmonics ΔC21, ΔS21, and ΔC20 changes are calculated as the sum of atmosphere and ocean mass effect (GAC values) taken from GRACE and a land surface hydrological estimate from the selected CMIP5 climate models. Low degree Stokes coefficients of the surface mass density determined from GRACE, SLR, GNSS, Earth rotation measurements and climate models are compared to each other in order to assess their consistency. The comparison is done by using different types of statistical and signal processing methods.

Global Bathymetric Prediction For Ocean Modeling and Marine Geophysics

NASA Technical Reports Server (NTRS)

Sandwell, David T.; Smith, Walter H. F.; Sichoix, Lydie; Frey, Herbert V. (Technical Monitor)

2001-01-01

We proposed to construct a complete bathymetric map of the oceans at a 3-10 km resolution by combining all of the available depth soundings collected over the past 30 years with high resolution marine gravity information provided by the Geosat, ERS-1/2, and Topex/Poseidon altimeters. Detailed bathymetry is essential for understanding physical oceanography and marine geophysics. Currents and tides are controlled by the overall shapes of the ocean basins as well as the smaller sharp ocean ridges and seamounts. Because erosion rates are low in the deep oceans, detailed bathymetry reveals the mantle convection patterns, the plate boundaries, the cooling/subsidence of the oceanic lithosphere, the oceanic plateaus, and the distribution of off-ridge volcanoes. We proposed to: (1) Accumulate all available depth soundings collected over the past 30 years; (2) Use the short wavelength (< 160 km) satellite gravity information to interpolate between sparse ship soundings; (3) Improve the resolution of the marine gravity field using enhanced estimates along repeat altimeter profiles together with the dense altimeter measurements; (4) Refine/improve bathymetric predictions using the improved resolution gravity field and also by investigating computer-intensive methods for bathymetric prediction such as inverse theory; and (5) Produce a 'Globe of the Earth' similar to the globe of Venus prepared by the NASA Magellan investigation. This will also include the best available digital land data.

Evaluation and Comparison of the Processing Methods of Airborne Gravimetry Concerning the Errors Effects on Downward Continuation Results: Case Studies in Louisiana (USA) and the Tibetan Plateau (China).

PubMed

Zhao, Qilong; Strykowski, Gabriel; Li, Jiancheng; Pan, Xiong; Xu, Xinyu

2017-05-25

Gravity data gaps in mountainous areas are nowadays often filled in with the data from airborne gravity surveys. Because of the errors caused by the airborne gravimeter sensors, and because of rough flight conditions, such errors cannot be completely eliminated. The precision of the gravity disturbances generated by the airborne gravimetry is around 3-5 mgal. A major obstacle in using airborne gravimetry are the errors caused by the downward continuation. In order to improve the results the external high-accuracy gravity information e.g., from the surface data can be used for high frequency correction, while satellite information can be applying for low frequency correction. Surface data may be used to reduce the systematic errors, while regularization methods can reduce the random errors in downward continuation. Airborne gravity surveys are sometimes conducted in mountainous areas and the most extreme area of the world for this type of survey is the Tibetan Plateau. Since there are no high-accuracy surface gravity data available for this area, the above error minimization method involving the external gravity data cannot be used. We propose a semi-parametric downward continuation method in combination with regularization to suppress the systematic error effect and the random error effect in the Tibetan Plateau; i.e., without the use of the external high-accuracy gravity data. We use a Louisiana airborne gravity dataset from the USA National Oceanic and Atmospheric Administration (NOAA) to demonstrate that the new method works effectively. Furthermore, and for the Tibetan Plateau we show that the numerical experiment is also successfully conducted using the synthetic Earth Gravitational Model 2008 (EGM08)-derived gravity data contaminated with the synthetic errors. The estimated systematic errors generated by the method are close to the simulated values. In addition, we study the relationship between the downward continuation altitudes and the error effect. The analysis results show that the proposed semi-parametric method combined with regularization is efficient to address such modelling problems.

Evaluation and Comparison of the Processing Methods of Airborne Gravimetry Concerning the Errors Effects on Downward Continuation Results: Case Studies in Louisiana (USA) and the Tibetan Plateau (China)

PubMed Central

Zhao, Qilong; Strykowski, Gabriel; Li, Jiancheng; Pan, Xiong; Xu, Xinyu

2017-01-01

Gravity data gaps in mountainous areas are nowadays often filled in with the data from airborne gravity surveys. Because of the errors caused by the airborne gravimeter sensors, and because of rough flight conditions, such errors cannot be completely eliminated. The precision of the gravity disturbances generated by the airborne gravimetry is around 3–5 mgal. A major obstacle in using airborne gravimetry are the errors caused by the downward continuation. In order to improve the results the external high-accuracy gravity information e.g., from the surface data can be used for high frequency correction, while satellite information can be applying for low frequency correction. Surface data may be used to reduce the systematic errors, while regularization methods can reduce the random errors in downward continuation. Airborne gravity surveys are sometimes conducted in mountainous areas and the most extreme area of the world for this type of survey is the Tibetan Plateau. Since there are no high-accuracy surface gravity data available for this area, the above error minimization method involving the external gravity data cannot be used. We propose a semi-parametric downward continuation method in combination with regularization to suppress the systematic error effect and the random error effect in the Tibetan Plateau; i.e., without the use of the external high-accuracy gravity data. We use a Louisiana airborne gravity dataset from the USA National Oceanic and Atmospheric Administration (NOAA) to demonstrate that the new method works effectively. Furthermore, and for the Tibetan Plateau we show that the numerical experiment is also successfully conducted using the synthetic Earth Gravitational Model 2008 (EGM08)-derived gravity data contaminated with the synthetic errors. The estimated systematic errors generated by the method are close to the simulated values. In addition, we study the relationship between the downward continuation altitudes and the error effect. The analysis results show that the proposed semi-parametric method combined with regularization is efficient to address such modelling problems. PMID:28587086

Evaluation and Comparison of the Processing Methods of Airborne Gravimetry Concerning the Errors Effects on Downward Continuation Results: Case Studies in Louisiana (USA) and the Tibetan Plateau (China)

NASA Astrophysics Data System (ADS)

Zhao, Q.

2017-12-01

Gravity data gaps in mountainous areas are nowadays often filled in with the data from airborne gravity surveys. Because of the errors caused by the airborne gravimeter sensors, and because of rough flight conditions, such errors cannot be completely eliminated. The precision of the gravity disturbances generated by the airborne gravimetry is around 3-5 mgal. A major obstacle in using airborne gravimetry are the errors caused by the downward continuation. In order to improve the results the external high-accuracy gravity information e.g., from the surface data can be used for high frequency correction, while satellite information can be applying for low frequency correction. Surface data may be used to reduce the systematic errors, while regularization methods can reduce the random errors in downward continuation. Airborne gravity surveys are sometimes conducted in mountainous areas and the most extreme area of the world for this type of survey is the Tibetan Plateau. Since there are no high-accuracy surface gravity data available for this area, the above error minimization method involving the external gravity data cannot be used. We propose a semi-parametric downward continuation method in combination with regularization to suppress the systematic error effect and the random error effect in the Tibetan Plateau; i.e., without the use of the external high-accuracy gravity data. We use a Louisiana airborne gravity dataset from the USA National Oceanic and Atmospheric Administration (NOAA) to demonstrate that the new method works effectively. Furthermore, and for the Tibetan Plateau we show that the numerical experiment is also successfully conducted using the synthetic Earth Gravitational Model 2008 (EGM08)-derived gravity data contaminated with the synthetic errors. The estimated systematic errors generated by the method are close to the simulated values. In addition, we study the relationship between the downward continuation altitudes and the error effect. The analysis results show that the proposed semi-parametric method combined with regularization is efficient to address such modelling problems.

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.